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ORGANISATION 
GENERALE DU COURS 



Ce cours du tronc commun "Dynamique des Fluides" du Master Recherche 
"Dynamique des Fluides, Energetique et Transferts" (DET), intitule " ON- 
DES ET INSTABILITES" , repond aux objectifs suivants : 

• Manipuler des modeles de base de la mecanique des fluides a travers 
l'etude des leurs ondes lineaires pour les systemes conservatifs et de 
leurs instabilites pour les systemes dissipatifs forces. Les principaux 
phenomenes physiques pris en compte sont les suivants : compressibilite, 
gradient vertical de densite, cisaillement vertical de vitesse ou encore 
presence d'une surface libre. 

• Assimiler plusieurs concepts de base souvent rencontres en mecanique 
des fluides : linearisation autour d'etats stationnaires, recherche de 
modes propres en tenant compte des symetries du probleme, calcul 
des relations de dispersions generalisees et interpretation physique des 
champs decomposes sur une base de modes propres. 

• Assimiler les concepts de base de la dynamique spatio-temporelles des 
trains d'ondes : relation de dispersion, vitesse de groupe d'un paquet 
d'onde, reponse impulsionnelle d'un milieu dispersif. 

• Assimiler les concepts de base de la dynamique instabilites temporelles : 
relations de dispersions generalisees, bifurcations de l'equilibre, insta- 
bilites de cisaillement ou thermiques dans les fluides. 

• Susciter un interet pour certains approfondissements : train d'ondes 
emis par un obstacle mobile ou oscillant, expression du flux d'energie 
ondulatoire et refraction des ondes, equations d'amplitudes, ... 

Le support ecrit de ce cours est structure sous la forme de plusieurs "Articles 
Pedagogiques Multimedia" independants. 
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Un "Article Pedagogique Multimedia" (APM) est un document electronique 
qui regroupe des textes de type polycopie mais aussi des animations, des 
programmes de demonstrations ou des illustrations difficiles a transcrire sous 
forme papier. 

Tous ces APM sont accessibles a l'adresse suivante : 

http : / / thual . perso . enseeiht . f r/ otapm 

Ce cours comporte deux parties decrites ci-dessous. 

PARTIE I : "Des ondes et des fluides" 

Cette partie s'appuie sur les cinq "Articles Pedagogiques Multimedia" (APM) 
suivants : 

• reladi : Relations de dispersion des ondes dans les fluides 

• paqund : Dispersion d'un paquet d'ondes ID 

• emiobs : Emission d'ondes par un obstacle oscillant ou mobile 

• fluxen : Flux d'energie des ondes dans les fluides 

• refrac : Refraction des ondes et trace de rayons 

Ces articles peuvent etre lus et travailles independamment les uns des autres. 

La partie ecrite de ces cinq articles font l'objet de l'ouvrage suivant : "Des on- 
des et des fluides" , O. Thual, Cepadues-Editions 2004. La version electronique, 
qui complete l'ouvrage par des animations ou la photocopie de transparents, 
est disponible a l'adresse : 

http : / / thual . perso . enseeiht . f r/ otapm 

Les pre-requis necessaires pour la comprehension de ce cours peuvent etre 
limites a une connaissance de base des equations de la mecanique des flu- 
ides. Parmi les tres nombreux livres permettant d'acquerir cette connais- 
sance, j'en cite un qui presente l'avantage de partager la plupart des nota- 
tions utilisees dans cet ouvrage : "Introduction a la Mecanique des Milieux 
Continus Deformables" , O. Thual, Cepadues-Editions 1997. 

PARTIE II : "Des instability et des fluides" 

Cette partie s'appuie sur les quatre "Articles Pedagogiques Multimedia" (APM) 
suivants : 



ONDES ET INSTABILITIES 



7 



• cisail : Instabilities de cisaillement 

• rayben : Convection de Rayleigh-Benard 

• bifgen : Bifurcations generiques de l'equilibre 

• sysdyn : Systemes dynamiques simples 

Ces articles peuvent etre lus et travailles independamment les uns des autres. 

La partie ecrite de ces quatre articles est contenue dans le present polycopie. 
La version electronique, qui complete ce polycopie par des animations ou la 
photocopie de transparents, est disponible a l'adresse : 

http : / / thual . perso . enseeiht . f r/ otapm 

SYLLABUS 

ONDES ET INSTABILITIES 

Semestre 1 Cours : 17 h TD : Pro jet : Credits : 

Mots-Cles : ondes sonores, ondes de gravite internes, ondes de surface, 
relations de dispersion, refraction, sillage d'onde. 

Bibliographie : [1] LIGHTILL, "Waves in fluids", Cambridge (1990), [2] 
DRAZIN et REID, "Hydrodynamics stability", Cambridge (1984), [3] Thual, 
"Des ondes et des fluides", Cepadues (2004), [4] Charru, "Instabilites Hydro- 
dynamiques", EDP Sciences (2007). 

Objectif : Assimiler les concepts generaux sur les ondes et les instabilites 
tout en manipulant des modeles de base de la mecanique des fluides. 

Programme : Le cours est repose sur une dizaine "dArticles Pedagogiques 
Multimedia" (APM) independants regroupes en deux parties. 

La partie intitulee "Des ondes et des fluides" comprend cinq APM : (1) 
derivation des relations de dispersion des ondes sonores, de gravite interne 
et de surface, (2) notion de vitesse de groupe a travers l'etude de paquets 
d'ondes, (3) train d'onde emis par un obstacle oscillant ou mobile, (4) calcul 
de Penergie d'une onde monochromatique et de son flux, (5) refraction des 
ondes et trace de rayons. 

La partie intitulee "Des instabilites et des fluides" comprend quatre APM : (1) 
derivation des relations de dispersion de plusiseurs instabilites de cisaillement, 
(2) calcul du seuil de la convection de Rayleigh-Benard, (3) dynamique des 
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bifurcations generiques de l'equilibre, (4) etude d'instabilites pour quelques 
systemes dynamiques simples. 

Le cours se concentre sur la moitie de ces APM tout en motivant la lecture 
des autres. Plusieurs exemples d'applications des notions introduites sont 
abordes sous forme d'exercices. 

O. THUAL 

DEROULEMENT PRATIQUE 

• Cours magistral : le cours oral proprement dit dure entre lh30 sans 
pause. 

• Evaluations : revaluation est effectuee a l'aide d'un examen a la fin 
du cours. Ce controle ecrit est individuel. Tous les documents sont 
autorises. 

• Documentation du cours : la documentation du cours est constitute 
de neuf "Articles Pedagogiques Multimedia" (APM). Chaque APM peut 
etre lu et travaille independamment des autres. Certains exercices et 
problemes poses lors d'examens precedents y sont inclus. Ces APM, 
avec leurs dernieres mises a jour, sont disponibles a l'adresse : 

http : //thual .perso . enseeiht . f r/otapm 

Ces pages contiennent aussi des animations illustrant le cours et des 
programmes Matlab a telecharger. 

La partie "Des ondes et des fluides" fait l'objet d'un ouvrage publie aux 
editions Cepadues (2004). Le pret de cet ouvrage peut-etre envisage 
par l'enseignant. La partie "Des instabilites et des fluides" pourra etre 
distribute sous forme de polycopie. 

• Travaux des eleves : certaines notions abordees dans ce cours sont 
utilisee lors de projets realise par les eleves de PENSEEIHT et publiees 
a l'adresse suivante : 

http : //www. enseeiht . f r/travaux 

• Bibliographie : pret dans les bibliotheques de plusieurs etablissement 
des ouvrages de la bibliographie du livre. 

• Fiche d'evalution du cours : les etudiants sont invites a remplir une 
fiche d'evaluation du cours (voir ci-dessous) pour la remettre a Tissue 
de l'examen final. 
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PROGRAMME DETAILLE DU COURS 

Pour l'annee 2009-2010, les "Articles Pedagogiques Multimedia" suivants 
seront exposes lors du cours magistral : 

1. bifgen : Bifurcations generiques de Pequilibre 

2. reladi : Relation de dispersion des ondes dans les fiuides 

3. rayben : Convection de Rayleigh-Benard 

4. paqund : Dispersion d'un paquet d'ondes ID 

5. cisail : Instabilites de cisaillement 

6. refrac : Refraction des ondes et trace de rayons 
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FICHE D 'EVALUATION DU COURS 
ONDES ET INSTABILITIES 

Afin d'etablir un bilan du cours et eventuellement, en fonction de vos reponses, 
d'envisager des modifications de l'enseignement pour les annees suivantes, je 
vous serais reconnaissant de bien vouloir remplir ce questionnaire et me le 
transmettre. 

Olivier THUAL 



NOM (facultatif : ) : 

Pour chacune des rubriques, portez un croix dans la colonne retenue et com- 
mentez votre reponse 





Tros Bien 


Bien 


Moyen 


Passable 


Mauvais 


Commcntaires 


Definition des 
objectifs du cours 














Documentation 
ecrite du cours 














Intervention 
de l'enseignant 














Controle 

des connaissances 














Atteinte des 
objectifs du cours 















Commentaires supplement aires : 
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La partie intitulee "DES ONDES ET DES FLUIDES" s'appuie sur les cinq 
"Articles Pedagogiques Multimedia" (APM) suivants : 

• disper : Relations de dispersion des ondes dans les fluides 

• paqund : Dispersion d'un paquet d'ondes ID 

• emiobs : Emission d'ondes par un obstacle oscillant ou mobile 

• fluxen : Flux d'energie des ondes dans les fluides 

• refrac : Refraction des ondes et trace de rayons 

Ces articles peuvent etre lus et travailles independamment les uns des autres. 

La partie ecrite de ces cinq articles font l'objet de l'ouvrage suivant : "Des on- 
des et des fluides" , O. Thual, Cepadues-Editions 2004. La version electronique, 
qui complete l'ouvrage par des animations ou la photocopie de transparents, 
est disponible a l'adresse : 

http : //thual .perso . enseeiht . f r/otapm 

Les pre-requis necessaires pour la comprehension de ce cours peuvent etre 
limites a une connaissance de base des equations de la mecanique des flu- 
ides. Parmi les tres nombreux livres permettant d'acquerir cette connais- 
sance, j'en cite un qui presente l'avantage de partager la plupart des nota- 
tions utilisees dans cet ouvrage : "Introduction a la Mecanique des Milieux 
Continus Deformables" , O. Thual, Cepadues-Editions 1997. 
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DES INSTABILITIES ET 
DES FLUIDES 
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La partie intitulee "DES INSTABILITIES ET DES FLUIDES" s'appuie sur 
les quatre "Articles Pedagogiques Multimedia" (APM) suivants : 

• cisail : Instabilites de cisaillement 

• rayben : Convection de Rayleigh-Benard 

• bifgen : Bifurcations generiques de Pequilibre 

• sysdyn : Systemes dynamiques simples 

Ces articles peuvent etre lus et travailles independamment les uns des autres. 

La partie ecrite de ces quatre articles est contenue dans le present polycopie. 
La version electronique, qui complete ce polycopie par des animations ou la 
photocopie de transparents, est disponible a l'adresse : 

http : / /thual . perso . enseeiht . f r/otapm 
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APM 1 

Instabilities de cisaillement 



FICHE SIGNALETIQUE 



Article Pedagogique Multimedia (APM) 


Titre 




Instabilites de cisaillement 


Auteur 




0. THUAL, INPT /ENSEEIHT 


Reference 




APM-INPT thu-cisail (2004) 


Mois de creation 




septembre 2004 


Mise a jour 




November 14, 2009 


Comite de lecture 




neant 


URL multimedia 




http : //www . enseeiht . f r/hmf /enseignants/thual/ 


Niveau 




Master I ou II 


Cours oral 




lh30 a 2h 


Animations 




neant 


Mots cles 




instabilites, croissance temporelle 




relations de dispersion generalisees 



OBJECTIFS PEDAGOGIQUES 

L'objectif de cette article pedagogique est de presenter la notion d'instabilite 
sur les exemples des exemple simples issus de la mecanique des fluides en se 
limitant aux ecoulement cisailles. Plusieurs objectifs de formation sont vises : 

• Assimiler la notion de linearisation autour d'un equilibre. 

• Assimiler la notion de relation de dispersion. 

• Saisir le sens physique du declenchement de l'instabilite "roll waves" 
d'un ecoulement a surface libre sur un plan incline. 

• Saisir le sens physique du declenchement de l'instabilite de Kelvin Hel- 
moltz. 
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NOTATIONS 



• Assimiler les notions de bases sur la stabilite des ecoulements paralleles. 

Les competence a acquerir lors de l'etude de cet article pedagogique sont les 
suivantes : 

• Linearisation d'un systeme d'equations aux derivees partielles autour 
d'un equilibre. 

• Recherche de la forme des modes propres en tenant compte des invari- 
ances par translations dans certaines directions d'espace. 

• Calcul de la relation de dispersion et des modes propres. 

• Discussion de la stabilite en faisant varier un parametre. 

PRE-REQUIS 

Le niveau requis pour la lecture de cette article pedagogique se situe autour 
de celui d'une Licence scientifique. 

NOTATIONS 



a Amplitude reelle 

A Amplitude complexe 

A Vecteur complexe 

A_i Vecteur propre de L(ik±, ik2,ikz) 

b= (bi, 62 , ■•■) Vecteur de composantes bi 

B Amplitude complexe 

Cf Coefficient de frottement 

W Plan complexe 

Derivee partielle par rapport a t 

-J^ Derivee partielle par rapport a x 

div Operateur divergence 

D Operateur 4- 

V Sous-domaine de l'espace complexe (U 

e x ,e y ,e z Vecteurs unitaires de la base canonique 

f(x, y, z, t) Fonction quelconque 

F Nombre de Froude 

F(x, t) Equation de la surface libre F(x, t) = 

F c Nombre de Froude critique 

2l(U) Operateur vectoriel, differentiel non lineaire 

g Gravite 

g' Notation g' = g cos a 
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grad Operateur gradient 

h(x, t) Hauteur d'une couche fluide 

h* Grandeur h adimensionnee 

h(x, t) Perturbation de h 

h* Perturbation h adimensionnee 

h n Const ante 

I Pente / = sin a 

ki "Vecteur d'onde" ID ou composante de k 

k = (ki,k2,ks) Vecteur d'onde 3D 

k Module \k\ ou norme ||fe|| 

k* Grandeur k\ adimensionnee 

K Nombre d'onde k adimensionne 

K Notation K = fa/Cf 

C Operateur differentiel P {j^j 

Ll(£i,&, 6) Matrice n x n des polynomes (^1,^2,6 

Lij Composante de 

0(b) Du meme ordre de grandeur que b 

P(£) Polynome 

p Constante ou parametre de controle 

p(x, t) Champ de pression 

pi , P2 Champs de pression dans les couches 1 ou 2 

po(z) Profil de pression 

p r Pression constante 

R Nombre sans dimension 

i?* Valeur particuliere de R 

R c Valeur particuliere de R 

M Droite des reels 

s Valeur propre complexe s = o — iu 

S Notation S = s/C f 

Si Valeurs propres complexes 

t Temps 

t* Temps adimensionne 

tanh Tangente hyperbolique 

u(x, t) Champ ID pour un modele lineaire 

uq Valeur d'equilibre de u 

u m (z) Profil complexe 

u Perturbation de u autour d'un equilibre 

U(x, t) Vitesse moyenne d'une couche fluide 

U* Grandeur U adimensionnee 

Hi , U.2 Champ de vitesse dans les couches 1 ou 2 

U±,U2 Vitesses constantes 

Uq(z) Profil de vitesse 
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NOTATIONS 



Uq(z) Derivee premiere de Uq(z) 

Uq(z) Derivee seconde de Uq(z) 

U(x, t) Perturbation de U 

U* Perturbation U adimensionnee 

U n Constante 

U(x, t) Champ 3D de dimension n 

U_q Solution stationnaire 

(U\, U2, •••) Composantes de U 

U(x, y, z) Champ de vitesse 3D 

(u, v, w) Composantes de U(x, y, z) 

U(x, t) Perturbation de U autour d'un equilibrc 

v Vitesse dans la direction e y 

v m (z) Profil complexe 

w Vitesse dans la direction e z 

w m {z) Profil complexe 

w(x, t) Vitesse verticale 

w\ , W2 Vitesses w dans les couches 1 et 2 

x Coordonnee sur l'axe e x 

x Variable d'un systeme dynamique 

x* Grandeur x adimensionnalisee 

x Partie reelle de z 

x e Equilibre d'un systeme dynamiques 

£+ Equilibre particulier 

X- Equilibre particulier 

x(t) Trajectoire solution d'un systeme dynamique 

X_ Variables d'un systeme dynamique 

X e Equilibre d'un systeme dynamiques 

X + Equilibre particulier 

X__ Equilibre particulier 

X_{t) Trajectoire solution d'un systeme dynamique 

y Coordonnee sur l'axe e y 

y Partie imaginaire de z 

y{t) Notation pour —x(t) 

z Coordonnee sur l'axe e z 

z = x + iy Variable complexe d'un systeme dynamique 

z(t) Trajectoire complexe solution d'un systeme dynamiques 

\z\ Module de z 

a Angle d'un plan incline avec l'horizontale 

a±,ct3 Constantes 

Constante ou parametre de controle 

7 Constante 

S Constante ou parametre de controle 
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A Operateur Laplacien 

A Discriminant d'une equation du second degre 

r](x,y,t) Elevation de la surface libre 

A(z, ki, fo) Operateur de derivation 

A(z, k\ Operateur A(z, ki, 0) 

ix Const ante 

A*o j ^2 , M4 Const antes 

£ Variable du polynome P(£) 

Ci > 6 , 6 Variables des polynomes L,j (£i ,6,6) 

p Module de z 

i) Masse volumique 

p\ , p2 Masses volumiques dans les couches 1 et 2 

u Partie reelle de s 

o\ , ci2 Valeur particulieres de a 

S(/ci) Partie reelle de la relation de dispersion generalisee 

Sj(fe) Partie reelle d'une relation de dispersion generalisee 

S + , S_ Parties reelles d'une relation de dispersion generalisee 

v Viscosite cinematique 

v Viscosite turbulente 

(f>(x, t) Potentiel des vitesses 

4>\ , 4>2 Potentiel des vitesses dans les couches 1 et 2 
$1(2), $2(2) Profils dans les couches 1 et 2 

uj Partie imaginaire de — s 

uji , UJ2 Valeur particulieres de a 

Cl(ki) Partie imaginaire de la relation de dispersion generalisee 

Q,i(k) Partie imaginaire d'une relation de dispersion generalisee 

Q- Parties imaginaire d'une relation de dispersion generalisee 



COURS ECRIT 
Introduction 

On considere ici les modeles d'ecoulement les plus simples permettant de 
mettre en evidence la notion d'instabilite de cisaillement. Dans un premier 
temps, on presente le concept de linearisation autour d'une solution station- 
naire et de croissance exponentielle d'une petite perturbation sur des exem- 
ples d'equations aux derivees partielles ID tres simples. Le premier exemple 
d'ecoulement physique presente est alors le modele des equations de Saint 
Venant qui decrit la dynamique des ondes de surface d'une couche fluide 
peu profonde coulant sur un plan incline. Le deuxieme exemple traitable de 
maniere analytique est celui de Pinstabilite de Kelvin Helmotz qui decrit le ci- 
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saillement de deux couches fluides irrotationelles, la vorticite etant concentree 
a l'interface. Enfin, le probleme au valeurs propres permettant de determiner 
la stabilite des ecoulements paralleles est aborde a travers la derivation de 
l'equation d'Orr-Sommerfeld. 

1 Modeles prototype ID 

Pour introduire les notions d'instabilite, de relations de dispersion ou de 
modes propres de la maniere la plus simple possible, on presentent ici des 
modeles bases sur des equations aux derivees partielles unidimensionnelles 
(ID). II est fait ici abstraction de l'eventuelle signification physique des solu- 
tions, bien qu'il serait possible d'expliciter des systemes dont la dynamique 
se laisse decrire par ces modeles. 

1.1 Modele lineaire 

Pour introduire le concept d'instabilite, on considere le modele lineaire ID 
suivant : 

du du d 3 u d 2 u <9 4 u 

En se plagant dans la cadre d'un domaine infini, on cherche des solutions 
complexes de la forme 

u(x,t) = Ae iklX+st = Ae iklX e at - iujt (1.2) 

avec A G (T, k\ G M et s = a — i u> 6 W. La forme de ces solutions est justifiee 
par le fait que l'equation aux derivees partielles consideree est invariante par 
translations d'espace et de temps. On pourrait considerer des solutions plus 
generale de la forme u(x, t) = a exp (7 x + i ki x) exp (at — icut), pour pren- 
dre en compte une croissance spatiale, mais on se restreint ici au cas ou les 
solutions sont bornees a l'infmi. Dans la mesure ou le modele est lineaire, 
les solutions generales sont obtenues en sommant les solutions elementaires 
considerees. En particulier, les solutions reelles peuvent s'obtenir en prenant 
la partie relle des solutions complexes. Pour pouvoir etendre les notations aux 
cas 2D ou 3D, on dira que k\ G M est le "vecteur d'onde" et que sa norme 
(valeur absolue) k = \k±\ est le "nombre d'onde". 

En reportant la forme des solutions recherchees dans l'equation, on voit que 
l'on obtient des solutions non triviales (^4 7^ 0) a condition que les relations 
de dispersions generalisees suivantes soient verifiees : 

a = T,(k\) = fiQ — fj,2 k\ + fj,4 k l et uj = Vl(k\) = a\k\ — 03 k\ . (1.3) 



1. MODELES PROTOTYPE ID 
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Un autre point de vue consiste a ecrire le modele sous la forme 2jr = C u, 

ou C = P (jj^j est l'operateur differentiel qui contient les derivee spatiales et 

P est le polynome defini par P(£) = no — "i £ + A*2 £ 2 — 03 £ 3 + A*4 £ 4 - Les 
fonctions Ae lkx sont les "modes propres" de l'operateur C. On peut ecrire 

£e iklX = P(ih)e iklX = [E(Jfei)-ifi(fci)] e ifclX . (1.4) 




Figure 1.1: Dispersion et stabilite du modele prototype dans le cas ol\ = 1, 
«3 = 1, fi2 = —2 et |U4 = —1. aj Taux de croissance a = S(/ci) pour 
/io 6 { — 1.5,-1.0,-0.5,0.0,0.5}. b) Relation de dispersion to = U(ki). c) 
Seuil d'instabilite R c (k) = — 1 dans la plan (k,Ho). 

La deuxieme relation de l'equation (1.3) (figure 2.2.b) est la relation de dis- 
persion de l'equation de Korteweg de Vries lineaire que l'on obtient dans le 
cas ou les coefficients /ij du modele sont nuls. 

La premiere relation de l'equation (1.3) (figure 2. 2. a) permet de determiner si 
les solutions du modeles sont amplifiees (<r > 0) ou amorties (<r < 0). On voit 
que l'on doit supposer \i± < pour que les modes de grand nombre d'onde k 
(petites periodes spatiales) soient amortis. Dans le cas contraire, le modele est 

non seulement "non physique" , mais aussi mal pose sur le plan mathematique. 

2 

Si H4 < 0, le maximum de S(fei) est atteint pour k\ = et vaut no — 

On voit alors que lorsque le nombre sans dimension R = 4 110 f" 1 est inferieur 

fj, 2 

a la valeur critique R c = —1, tous les modes sont stables. Lorsque R est 
superieur a cette valeur critique, les modes sont instables pour un intervalle 
finie de nombre d'onde fei, le mode le plus instable correspond a la valeur 
critique k c \J-^- 

1.2 Linearisation autour d'un equilibre 

On considere le modele non lineaire ID suivant : 
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avec u(x, t) G M, 7 > 0, U2 > et U4 < 0. Pour u > 0, on remarque 
l'existence des deux solutions stationnaires symetriques n+ = no et n_ = —no 
avec no = ^/n/7. On s'interesse a la stabilite temporelle de equilibre n+ = no 
au regard des petites pertubations. Pour cela, on pose 

u(x, t) = Uq + u(x, t) (1-6) 

et l'on s'interesse a revolution des perturbations n(x, t) que l'on peut choisir 
petites a un instant initial donne. En reportant dans l'equation du modele 
on obtient : 

du du ^du ( 2 \ ~ ^ 2 n ^ u 

-7 3n n 2 + n 3 . (1.7) 

Tant que la perturbation u(x, t) reste petite, on peut negliger les terme "non 
lineaires" u || et [3nQU + n 3 ] devant les termes lineaires. En posant uo = 
u — 3 7^ = — 2 u, le systeme "linearise autour de l'equilibre no" s'ecrit done 

du du „ d 2 u d 4 u 

-+n - = u n + a 2 ^ + u 4 ^. (1.8) 

On retrouve le modele lineaire etudie dans le paragraphe precedent a condition 
de poser a\ = no et 03 = 0. 

1.3 Calcul de stabilite dans le cas general 

Pour generaliser le concept de stabilite temporelle d'un equilibre au cas 3D, 
on considere le modele general 

f = < 1J » 

oil U(x, t) = (Ui, U2, U n ) G M n est un vecteur dans un espace de dimension 
n, x = (x±, X2, X3) G M 3 est le vecteur des coordonnees spatiales et T_ est un 
operateur non lineaire faisant intervenir les derivees spatiales (gf^-, gf^, g§^) 
des composantes de U. 

On suppose que l'etat constant U est un equilibre du modele, e'est-a-dire 
qu'iljverifie T(U ) = Q. On pose alors U(x, t) = U + U(x, t) et l'on suppose 
que U_ est une perturbation initiallement petite. En reportant dans le modele 
on obtient _ 

f (■*')• 

ou L(£i, ^2 j £3) est une matrice nxn dont les composantes (£1, £2, £3) sont 
des polynomes et O (j|ZZ|| 2 ) est un terme negligeable tant que la perburbation 
est petite. 
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On cherche alors des solutions du systeme linearisee j^U = L U sous la forme 

U(x,t) = Ae i ^ +8t (1.11) 

ou k = (ki, k2, ks) € M 3 est appele le "vecteur d'ondes" et A est une amplitude 
complexe constante. On obtient alors 

s A = LXik 1 ,ik 2 ,ik 3 ) A . (1-12) 

Les solutions non nulles sont obtenues en choissant s, k et A de telle sorte 
que .A soit un vecteur propre de la matrice L(iki,ik2,iks). La valeur propre 
est alors s et on dit que la solution ainsi trouvee est un "mode propre" . 
Dans le cas ou L(iki,ik2,iks) est une matrice (n x n) diagonalisable, cette 
recherche de valeurs propres et vecteurs propres conduit a n relations de 
dispersion generalisees de la forme 

<7j = Sj(fe) et u>i = £li(k) pour i = 1, n . (1-13) 

Les vecteurs propres Ai associes , solutions de L A { = Sj A i sont les amplitudes 
complexes des modes propres A i exp(i k-x) qui engendrent des solutions dont 
le comportement temporel est donne par la fonction e Sit = e a i t ~ tuJ i t . 
On en deduit que cij = Sj(fc) est negatif pour les i = l,...,n et tous les 
vecteurs d'ondes k consideres, le systeme lineaire est stable et, par consequent, 
Pequilibre ?7 autour duquel le modele a ete linearise. En revanche, il suffit 
qu'il existe une valeur propre Sj dont la partie reelle U{ est positive pour une 
bande, meme petite, de vecteurs d'ondes k, pour que l'equilbre considere soit 
instable. 



2 L'instabilite "roll waves" d'un ecoulement incline 

L'ecoulement d'une lame d'eau dont l'epaisseur est faible devant l'echelle hor- 
izontale des phenomenes que Ton observe est modelisee par les equations de 
Saint Venant (shallow water). Dans le cas d'un ecoulement gravitaire sur un 
plan incline, on observe une instabilite de la surface libre lorsque la pente 
depasse un seuil critique. Cette instabilite porte le nom de "roll waves" du 
nom des structure en forme de rouleaux qui se forment a la surface libre. Nous 
limitons ici l'etude au calcul de seuil de declenchement de l'instabilite dans le 
cas le plus simple. 

2.1 Modele des equations de Saint Venant 

On considere un ecoulement a surface libre sur un plan incline faisant un angle 
a avec Phorizontale. On note g la gravite. 
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Figure 1.2: Ecoulement a surface libre d'une lame d'eau d'epaisseur h et de 
vitesse moyenne U sur un plan incline. 

Dans le cadre de "P approximation des eaux peu profondes" , le modele decrivant 
l'ecoulement est celui des equations de Saint Venant qui s'ecrivent 



dU u dU 

dt dx 
9h d /TT1 . 

m + d- x {Uh) 



-g' ^+gl 
Ox 



1 U\U\ 



o, 



(1.14) 



ou U(x,t) est la vitesse moyenne de la couche fluide, h(x,t) sa hauteur, g' = 
g cos q la projection de la gravite sur la normale au fond et / = sin a le 
sinus de l'angle d'inclinaison, proche de la pente lorsque celle-ci est petite. Le 
nombre sans dimension Cf modelise le frottement du fond sur la couche. On 
suppose ici que Cf est constant. 

Les solutions stationnaires (U n ,h n ) de ce modele sont solutions de l'equation 



U n 



(1.15) 



II existe done une famille a un parametre de solutions stationnaires que l'on 

appelle "ecoulements normaux" en hydraulique. 

On cherche a determiner la stabilite de ces equilibres. On pose alors 



U(x,t) = U n + U(x,t) 



et 



h(x, t) = h n + h(x, t) . 



(1.16) 



En reportant dans le modele est en negligeant les termes d'ordre 2, on obtient 
le systeme linearise suivant 



dU dU 

~dt +Un dx~ 

dh dh dU 
at ox ox 



~ 9 dx +C f 



Un 

K 



u 




2 hi 



(1.17) 



On cherche alors des solutions sous la forme 



(U,h) = (A, B) e ik ' x+st 



(1.18) 
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oil A et B sont de amplitudes complexes, et s = a — i to £ W. En reportant 
dans les equations, on obtient le systeme lineaire suivant : 



ihK s + ikiU n / \BJ \0 



(1.19) 



On obtient des solutions non triviales a condition que le determinant de la 
matrice 2x2 soit nul, ce qui conduit a la relation de dispersion generalised 
suivante : 

(s + ihUn + Cf^ (s + ihUrJ-ihg'K (d'ih-Cf^j = 0. 

(1-20) 

La stabilite des equilibres (U n , h n ) s'obtient en trouvant les racines complexes 
s\ et S2 de ce polynome de degre deux pour chaque valeur de k±, dont on 
deduit les relations de dispersion generalisees 

£7i = Si(fei), LJi = fti(fci) et a 2 = S 2 (fci), ^2 = ^2(h) . (1.21) 

A priori, la stabilite depend des quatre parametres (h n ,g',I,Cf), la vitesse 
U n y etant reliee. L'analyse dimensionnelle permet de reduire la dependance 
a deux parametres sans dimensions, dans la mesure oil les unites de longueur 
et de temps retirent deux parametres. Nous allons retrouver ce resultat dans 
le paragraphe suivant en ecrivant le modele sous une forme adimensionnee. 
Nous allons ensuite montrer que la structure particuliere des equations fait 
qu'il n'y a en fait qu'un seul nombre sans dimension a considerer. 

2.2 Equations adimensionnees 

Pour etudier la stabilite d'un ecoulement stationnaire (U n , h n ) on choisit 
d'adimensionner les grandeurs de la maniere suivante 

- h + - h^-t 

X t^Tl X% 5 ^ T T 

U(x,t) = U n U*(x*,t*), h(x,t) = h n h*(x*,U) , (1-22) 

oil (x*, t*,U*,h*) sont des grandeurs sans dimensions. En reportant ce change- 
ment de variables dans les equations on obtient le systeme adimensionne suiv- 
ant : 

dU. dU. dK i c ( _ uf\ 

du * dx* F 2 dx* 2 f { K) 

?]± + JL{U.h.) = 0, (1.23) 
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oil F = 



— — ./M- est le "nombre de Froude". Pour ces equations adi- 

mensionnees, l'equilibre considers? est (U*,h#) = (1,1). 
On etudie sa stabilite en posant 

j ) — 1 ~\~ XJ* (yX^f. j ij^f. ) ct h% 5 ) — 1 ~~ \~ (yXjf. ^ ) (^1.24) 

et en cherchant revolution des perturbations sous la forme 

(U*X) = (A*,B*)e ik * x * +s * u . (1.25) 

La stabilite de Pequilibre ([/*, h*) = (1, 1) s'obtient en annulant le determinant 
de l'operateur linearise, ce qui s'ecrit 



~t~ i -\- C j p2 ^ k* 2 ^ f 
i k* ~\~ % k* 

(s* + i k* + Cj) (s* + i ) i /c* ( i k* — Cf 



. (1.26) 



On pose S = s*/C/ et if = k*/Cf. Avec ces notations, la relation de disper- 
sion generalised s'ecrit 



(S + iK + l)(S + iK) + ^K 2 + ^iK = . 



(1.27) 



On voit que le parametre Cf ne joue aucun role dans la determination du 
seuil d'instabilite. 



2.3 Calcul de stabilite 

La stabilite de l'equilibre ne depend done plus que du nombre de Froude 
F. En posant Y = S + i K pour simplifier la recherche des racines de cette 
equation du second degre, on se ramene a 

Y 2 + Y + ^K 2 + \iK = . (1.28) 

Les solutions s'ecrivent 

Y l {K)= l 2 (-l + ^) 
avec A 

la convention utilisee pour prendre la racine du discriminant complexe Aeff 
etant indifferente dans la mesure ou Ton s'interesse aux deux solutions. 
En posant S = T,(K) — i 0,(K) = —i K + Y(K), on peut tracer sur la figure 
1.5 les relations de dispersions generalisees Y,-(K) et T, + (K) pour la partie 



et Y 2 (K) = 
4 , 

l-2iK - —,K' 
F 2 



(- 



1 - A2 



(1.29) 



2. L'INSTABILITE "ROLL WAVES" D'UN ECOULEMENT INCLINE 31 




Figure 1.3: Relations de dispersion generalisees pour les nombres de Froude 
F = .2 : .2 : 3. a) E+(K) et E-(K). bj fi + (ff) et fl-(K). 



reelle et fl-(K) et 0+(if) pour la partie imaginaire, en adoptant la convention 
E_(tf) < E+(#). ^ 

En effectuant des developpements limites, on obtient les comportements asymp- 
totiques suivants : 




Sur le trace des courbes, on remarque que la courbe E_(iiT) est toujours 
negative tandis que la courbe T, + (K) est negative pour les nombres petits et 
positive, au moins pour certaines valeurs de K, pour F assez grand. Pour le 
Froude critique F = F c = 2, on remarque que les racines Y sont solutions de 
l'equation 

Y 2 + Y + -K 2 + -iK = ( Y + i - K ) 
4 2 V 2 y 

Pour F = 2, on a done 

E + (ir) = 0, = et E_(K) = -1, 0_(K) = ^ET. 

(1.33) 

On montre alors que E+(if) < pour F < 2 et E+(i^) < pour F > 2. Le 
seuil d'instabilite est done F = F c = 2. Tous les modes deviennent instables 



1 



i-K 



(1.32) 
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au seuil. Au seuil, la vitesse de phase des ondes est C + (K) = VL + (K)/K = | 
et C-(K) = Q-(K)/K = \. En revenant aux grandeurs dimensionnelles, 
ces vitesses sont done | U n et \ U n . On remarque que ces ondes sont non 
dispersives (vitesse de phase uj/k constante) au seuil F = F c = 2, alors que 
les modes stables ou instables sont dispersifs (quoique seulement legerement 
des que k est sumsamment grand). 




Figure 1.4: Mise en evidence experimentale de roll waves. 

Les modes instables que l'on voit apparaitre pour F > F c = 2 sont des struc- 
tures periodiques que Ton appelle "roll- waves" . A priori, tous les modes de 
toutes les longueurs d'ondes deviennent simultanement instables pour F = 2, 
les modes de petites longueurs d'ondes (k grand) etant les plus instables. 
Cette situation n'est pas physique et il faut prendre en compte, dans un 
modele plus complet, les effets de dissipation d'energie a petite echelle. Un 
exemple de modelisation des petites echelles, souvent turbulentes, consiste a 
ajouter un terme en v e dans le second membre de l'equation de conser- 
vation de la quantite de mouvement, v e etant une viscosite turbulente que 
l'on peut, par exemple, supposer constante. La saturation de l'instabilite 
s'effectue par la formation de ressauts (chocs) dont Pepaisseur est controlee 
par la valeur de la viscosite. Ces considerations ne sont pas developpees ici. 

3 Instability de Kelvin Helmoltz 

Lorsque deux couches fluides de vitesses et / ou de masses volumiques differentes 
sont mises en contact, on observe, pour certaines valeurs des parametres, 
le developpement d'une instabilite. Ce phenomene est appele, de maniere 
generique "instabilite de Kelvin Helmoltz". Plus specifiquement, cette deno- 
mination se rapporte a la situation idealisee de l'interface de deux fluides 
parfaits individuellement homogenes. 

3.1 Equations d'Euler incompressible avec interface 

On considere deux couches fluides semi-infinies separees par une surface libre. 
On suppose que les fluides sont parfaits et incompressibles et on note p\ et p2 
les masses volumiques des fluides respectivement situes en bas et en haut. On 
suppose que les fluides sont parfaits, si bien que leur dynamique est modelisee 
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-Z A 



-► u 2 



P2 



pi 



Figure 1.5: Couches fluides superposees de vitesses U et de masses volumiques 
p differentes. 



par les equations d'Euler incompressibles 

d 



-grad Pi- Pi g e z 



div U 1 = et pi — C/i 
div f7 2 = et p 2 -jU_ 2 = -grad p 2 - P2 9 e z 



(1.34) 



ou Ui(x, t) = (ui,vi,w\) et C/ 2 (x, £) = (u 2 , ^2,^2) sont les champs de vitesse 
et pi(x, t) et P2(x, t) les champs de pression. Le vecteur unitaire e z est vertical 
et g est l'intensite de la gravite. La notation A designe la derivee particulaire 



A. 

dt 

du haut. 



dt 

§i+LLi- giad pour le fluide du bas et ^ = + U_ 2 ■ grad pour le fluide 



On suppose que Ton a les conditions aux limites 



lim ~U_ X = U\ e. 



et 



lim Un 



U 2 e. 



(1.35) 



ou e x est un vecteur unitaire horizontal. 

Les conditions aux limites a l'interface des deux fluides sont 



dr] 
dt 



+ Hi ■ g?ad 7] = wi, Pl =p 2 , 



dr] 
dt 



+ U 2 ■ grad r] = w 2 (1.36) 



sur la surface libre d'equation F(x,t) = z — i](x,y,t) = 0. L'elevation de 
la surface libre est done r](x,y,t). Les deux conditions cinematiques ^ = 
w\ — = w 2 — ^j=0 s'obtiennent en prenant la derivee particulaire de la 
fonction F(x, t) = 0. La condition dynamique p\ = p 2 s'obtient en assurant 
la continuity des efforts de contact et en negligeant done l'effet de tension 
superficielle. 

On suppose que Pecoulement est irrotationnel dans chacune des couches flu- 
ides (rot Hi = rot H 2 = 2)) ce permet d'ecrire les champ de vitesses sous 
la forme 



U_ 1 = U\X + grad </>i et 
Le systeme d'equations s'ecrit alors 

A^i = et grad 



U_ 2 = U 2 x + grad fa ■ 



(1.37) 



tyl . TT 94>i 1 2 P! 

QT + U 1 — + -(grMl ( t>i) +-+9Z 
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A02 = et grad 
avec les conditions aux limites 



)<h . TT d(/> 2 1 . . 2 . P2 . 

df + U2 ^ + 2 (grJLd ^ 2) + J 2 +9Z 



lim grad 0i = et lim grad 02 = Q 



+ ) ?? + grad 02 • grad 77 



dz 



dz 



et 



Pi = P2 



en 2 



Comme 0i et <p2 sont respectivement dermis a une "constante" (en espace) 
Ci(t) et C*2(t) pres, on peut choisir une pression de reference p r arbitraire 
permettant d'eliminer la pression en ecrivant 



Pi = Pr ~ Pi 

P2=Pr~ P2 



<90i d0i 1 ,,s2. 
^ + f/i^ + 2(giad0i) +gz 

d(p2 dfa 1 2 , 



(1.38) 



3.2 Linearisation autour de l'equilibre 

On s'interesse a l'etat de base \J_ X = U\ e x , U 2 = U2 e x et rj = 0. La pression 
est alors po(z) = Pr — Pig z pour z < et po(^) = Pr — P2g z pour z > ou 
p r est une pression de reference arbitraire qui peut etre choisie sufnsamment 
grande pour eviter des pressions negatives dans une grande partie de la couche 
du haut qui est, en pratique, d'extension finie (mais grande). 
On linearise alors autour de cet etat de base en posant p\ = po(z) + p~\. La 
linearisation conduit a negliger les termes non lineaires dans les equations mais 
aussi a remplacer les conditions aux limites sur la surface mobile d 'equation 
z = r](x,y,t) par des conditions aux limites sur la surface fixe d'equation 
z = 0. En effet, pour un champ quelconque f(x, y, z, t) on peut ecrire 



f[x, y, rj(x, y, t),t] = f(x, y, 0, t) [1 + 0( V )} . 

Le modele est alors constitue de l'equation de Laplace A(f>i 
le fluide avec les conditions aux limites : 



(1.39) 

A(f>2 = dans 



lim grad <pi 

z^—00 

d d \ d<h 

dt + Ul dx: r = ^z- 



et lim grad 02 = Q 

( 9 TT 9 \ 

Pl=P2 ' ldt + U2 dx: r 



dz 



en z = . 



En remplagant la pression par sa valeur et en linearisant on remplace la con- 
dition p\ = P2 par la condition aux limites 



Pi 



— + U — 
dt 1 dx 



h + gv 



P2 



— + U — 
dt 2 dx 



^>2+gr] 



en z = . 
(1.40) 
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Dans la mesure oil le probleme est invariant par translations en temps et en 
espace dans les directions horizontals, on cherche des solutions sous la forme 



(01,77,02) = [* 1 (z),A,* 2 (z)]e iklX+ih *y + ' t . 

On note ici k = \Jk\ + k 2 . Le probleme a resoudre est done 

$'/ - k 2 $i = et $2 - A; 2 <J> 2 = avec 

(s + ik 1 U 1 )A = <f>' 1 (0) et (8 + ik 1 U 2 )A = & 2 (0) 

pi [(s + i ki Ui) <J>i(0) + g A) = p 2 [(s + ik 1 U 2 )<f> 2 (0) + gA] 

lim $[(z) = et lim & 2 (z) = 

On en deduit ®i(z) = Ai e kz et ® 2 (z) = A 2 e~ kz avec 



(1.41) 



(1.42) 



(s + iki Ui)A = A-y k et (s + i k x U 2 )A = -A 2 k 
pi [(s + i ki Ui) Ai + g A] = p 2 [(a + %k x U 2 ) A 2 + g A] (1.43) 

En reportant les valeurs de A\ et A 2 en fonction de A et en simplifiant par A 
on obtient la relation de dispersion generalisee 



Pi 



gk + is + ikiUyf 



P2 



gk-{s + ikiU 2 f 



(1.44) 



En resolvant cette equation de second degre, on obtient deux branches de 
solutions que Ton peut noter 



[Si(fei, k 2 ),tti(ki, k 2 )] et [S 2 (fei, k 2 ),n 2 (ki, k 2 )] . 



(1.45) 



3.3 Calcul de stabilite dans deux cas particuliers 

On montre qu'une condition necessaire et sumsante pour que Pequilibre soit 
instable est 

g y/kf + k\ [p\ - pf) < k\ pi p 2 (Ui - XJ 2 f . (1.46) 

II suffit, par exemple, de poser X = is dans (1.44) et de calculer le dis- 
criminant de l'equation du second degre en X ainsi obtenue. On voit done 
Pecoulement est instable pour que U\ / U 2 et que les modes de petites de 
longueur d'ondes (grands k) sont les plus instables. Une modelisation plus 
physique consiste alors a prendre en compte la viscosite qui a pour effet de 
dissiper l'energie et aux petites echelles. 

Plutot que d'etudier de maniere systematique tous les cas possibles, nous 
presentons ici quelques cas particuliers. Pour simplifier, on s'interesse aux 
ecoulements 2D, ce qui revient a poser ^2 = 0. On a done k = \ki\. En 
utilisant la symetrie du au fait que le probleme est reel, on peut alors n'etudier 
que le cas des vecteurs d'ondes k\ > 0. 
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Un premier cas particulier consiste a supposer que p\ = p 2 . Le fluide est 
homogene et l'instabilite ne provient que du cisaillement des vitesses. La 
relation de dispersion generalisee s'ecrit alors 

(a + ih C/i) 2 + (s + ih U 2 ) 2 = . (1.47) 

On en deduit les deux branches de modes suivants 

Ei(fci) = y (U! - U 2 ), ^(h) = | (Ui + *7 2 ) 

S 2 (fci) = y (172 - l/i), n 2 (Ari) = y (C/i + tf 2 ) (1-48) 

On voit que ces deux modes se progagent a la vitesse de phase c = (Ui + U 2 ) /2 
qui est la vitesse moyenne des deux couches fluides. L'un des deux modes est 
instable des que U\ ^ U 2 tandis que l'autre est stable. 

Un deuxieme cas particulier consiste a supposer que U\ = U 2 - En se plagant 
dans le repere qui se deplace a cette vitesse commune, on peut se ramener au 
cas U\ = U 2 = 0. La relation de dispersion generalisee s'ecrit alors 

Pi (gk + s 2 ) =p 2 (gk-s 2 ) s 2=Pl^£± g k. (1.49) 

V ) V J pi+ p 2 

Si p\ est plus grand que p 2 les deux racines s\ et s 2 sont imaginaires pure 
et le sytemes est marginal. Les ondes qui se developpent a l'interface des 
deux fluides sont des ondes de gravite internes. Si p 2 est plus grand que pi, 
c'est-a-dire si le fluide du haut est le plus lourd, il existe une families des 
modes instables dont le taux de croissance croit avec k. II faut alors enrichir 
le modele avec des termes de dissipation qui vont stabiliser les tres petites 
echelle est selectionner l'echelle des ondes instables. 



4 Stabilite des ecoulements paralleles 

On s'interesse ici a la stabilite des ecoulements plans paralleles dont les 
vecteurs vitesse pointent dans une meme direction et varient dans une di- 
rection orthogonale. Les instabilites de cisaillement ainsi obtenues permet- 
tent d'expliquer de nombreux phenomenes ce transition de la mecanique des 
fluides. 

4.1 Equations d'Orr-Sommerfeld 

On se place dans le cadre des ecoulements de fluide visqueux incompressibles 
qui sont decrits par les equations de Navier-Stokes 

d 1 

divC7; = et — U + U ■ grad U = grad p + v AU (1.50) 

dt po 
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ou U(x, t) = U(x, y, z, t) = (u, v, w) est le champ de vitesse, po la masse 
volumique constante, p la pression et v la viscosite cinematique. On suppose 
que l'ecoulement est compris entre les plans z = z\ et z = Z2, les bornes z\ 
et Z2 etant eventuellement infinies dans le cas de domaines non bornes. Dans 
le cas visqueux (y ^ 0), on impose les conditions aux limites U = U\e x en 
z = z\ et U = U2e x en z = z-i- Dans le cas inviscide (y = 0), seules les 
conditions aux limites cinematique w = en z = z\ et z = Z2 subsistent. 
On s'interesse aux solutions particulieres de la forme U = Uo(z) e x , ou Uq(z) 
est un profil de vitesse ne dependant que de z et e x un vecteur perpendiculaire 
a e z . On appelle "ecoulement parallele" une telle solution stationnaire dans 
la mesure ou les trajectoires et les lignes de champ de la vitesse (confondues 
dans le cas stationnaire) sont des droites paralleles (a e x ). 
En reportant cette forme de solution dans les equations de Navier-Stokes, la 
loi de conservation de la masse div U = est automatiquement satisfaite et 
la loi de conservation de la quantite de mouvement entraine ^ 
jpjj = et g| = 0, ou Uq designe la derivee seconde de Uq(z). 
Lorsque la viscosite moleculaire v est negligeable, les equations s'ecrivent 
M = % = If' ce ^ implique que p = p re f est une constante. Aucune 
contrainte n'est imposee pour la vitesse si bien que tous les profils Uo(z), de 
formes quelconques, sont solutions. 

Dans le cas ou l'ecoulement est visqueux, les equations entrainent ^ = G, 
|2 = ^2=0etG' = pqvUq{z), ou G est une constante arbitraire. On 
en deduit que p = po(x) avec po(x) = p TC { + G x, ou p re f est une constante 
arbitraire, et Uq(z) = n G z 2 + b z + c ou b et c sont des constantes arbitraires 
ou bien determinees par d'eventuelles conditions aux limites. 
Pour etudier la stabilite de ces ecoulements paralleles, on pose 



U(x,t) = U (z) e x + U(x,t) et p(x, t) = p (x) + p(x, t) (1-51) 
et Ton suppose que U = (u,v,w) et p sont des petites perturbations. En 
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Figure 1.6: Ecoulement parallele. 
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reportant dans les equations et en ne retenant que les termes lineaires, on 
obtient : 



dwU = et -j= + U (z) e x +wU' Q {z)e a 



1 

Po' 



gradp+z^ AU . (1.52) 



Compte-tenu de l'invariance de ce systeme lineaire vis-a-vis des translations 
en temps et dans les direction x et y, on cherche des solutions sous la forme 



[U(x,t),p(x,t)] = [U m (z),p m (z)] e 



i ki x+i ki y+st 



(1.53) 



avec U m = (u mjV m ,w m ) et s — (7 - iu £ ff. En reportant cette forme de 
solutions dans le modeles linearise, on obtient le systeme 

Ufa \ 

>o 1 \ 



(s + A(z,fa,k 2 ) o U' Q (z) 

s + A(z,ki,k 2 ) 
s + A(z,k 1 ,k 2 ) 

V i k\ ik 2 D 



po 

Mk 2 

PO 





\ 


/u\ 






V 







w 





/ 


\pj 





(1.54) 



avec A(z, fa, k 2 ) = i fa U (z) + u(D 2 - fa\), k 2 H = k\ + k\ et D = £. 
Les modes instables, s'ils existent, sont ceux pour lesquels la valeur propre s 
admet une partie reelle positive. Le theoreme de Squire indique que pour tout 
mode instable tridimensionnel (3D), il existe un mode bidimensionnel (2D) 
dans le plan (x, z) qui est encore plus instable. Ce theoreme se demontre en 
effectuant la transformation de Squire qui consiste a remplacer kjj par fa, 
fau + k 2 v par fa v\ et knp par fa p pour passer du cas 3D au cas 2D. Nous 
l'admettrons ici pour nous concentrer sur les perturbations 2D de la forme 



[U(x,z,t),p(x,z,t)} = [u m (z),v m (z),p m (z)} e iklX+st 
qui verifient done les equations 

1 

= 



(1.55) 



[s + A(z, fa)]u m + Uq(z) w m H i fa Pm 

1 po 

[s + A(z, fa)]w m H Dp m = et i fa u m + D w m = 0(1.56) 

Po 

avec A(z, fa) = ifa Uq(z) + v(D 2 — kf). En eliminant p m et u m de ce systeme, 
on obtient nnallement l'equation d'Orr-Sommerfeld qui s'ecrit 



s + ifa Uo(z) - v (D 2 — k 2 )} (D 2 - kj) w m = Ud(z) w m . 



(1.57) 



Les quatres conditions aux limites u m (Zl) = U m (z 2 ) = W m (zi) = W m (z 2 ) = 
sur les deux parois planes ou a l'infini conduisent aux quatres conditions 
w m (zi) = w m (z 2 ) = et w' m (zi) = w' m (z 2 ) = pour cette equation differen- 
tielle ordinaire de degre quatre en z. 



4. STABILITE DES ECOULEMENTS PARALLELES 
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Pour un vecteur d'onde donne fei, ce systeme n'admet de solutions triviales 
que pour un ensemble discret de valeurs propres s n {k\) si le domaine est borne 
ou pour un ensemble discret de valeurs propres s n {k\) combinee a un spectre 
continu s £ T>{k\) ou T>{k\) est une courbe ou une surface dans le plan com- 
plexe (D. La determination des valeurs propres et vecteurs propres associes 
a un ecoulement parallele Uq(z) est un probleme difficile et Ton a souvent 
recours a l'outil numerique pour le resoudre. En principe, seuls les profils 
Uq{z) de forme paraboliques peuvent etre considered dans le cas visqueux. En 
pratique, on etudie l'equation d'Orr-Sommerfeld des profils Uq(z) de formes 
quelconques c'est-a-dire pour la classe des solutions des equations de Navier- 
Stokes inviscides (equations d'Euler). Cette approximation revient a negliger 
le terme de dissipation visqueuse v Uq(z) dans la mesure ou il modifie la 
forme de Uo(z) sur des echelles de temps longues par rapport au temps car- 
acteristique de croissance des instabilites. 

4.2 Theoreme du point d'inflexion dans le cas inviscide 

Dans le cas inviscide v = 0, les equations de Navier-Stokes sont appelees 
"equations d'Euler" et l'equation d'Orr-Sommerfeld est appelee "l'equation 
de Rayleigh" qui s'ecrit 



Le quatres conditions aux limites sur le profil w m {z) de l'equation d'Orr- 
Sommerfeld se reduisent aux deux conditions w m {z\) = w m {z2) = 0, dans la 
mesure ou le nombre de conditions aux limites est reduit lorsque Ton passe des 
equations de Navier-Stokes aux equations d'Euler. Par exemple, une condition 
d'adherence a la paroi U_ = pour le cas des equations de Navier-Stokes se 
reduit a une condition cinematique w = lorsque l'on prend la limite v = 
des equations d'Euler. 
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Figure 1.7: Spectre de l'operateur dans le plan complexe. 
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Comme on l'a deja indique, tous les ecoulements paralelles Uq{z) e x sont 
solutions dans le cas inviscide, de maniere exacte. Un theoreme, demontre par 
Rayleigh en 1880, indique qu'une condition necessaire pour qu'un ecoulement 
parallele soit instable est l'existence d'un point d'inflexion dans le profil Uq(z). 
Pour obtenir ce result at, il suffit de multiplier l'equation de Rayleigh par 
Wm(z), le complexe conjuge de w m (z). En integrant P equations ainsi obtenu 
sur le domaine compris entre les plans z = z\ et z = z-i, on obtient : 



f 



Z2 w* m (D 2 - k\) w m dz = i h r U [ (Z 1 ^ m|2 dz . (1.59) 

J Zl S + IKiUq\Z) 



En integrant par parite le terme w* m D w m et en utilisant les conditions aux 
limites w m (zi) = w m (z2) = 0, on obtient 

f Z2 (\Dw m \ 2 + kl\w m \ 2 ) dz + l k, r U [ {Z 1 ^ dz = (1.60) 
hi v ' Jzi s + ikiU (z) 

dont la partie imaginaire est 



ok ^\ nrmr-rr 7^2 dz ( 1 - 61 ) 



Zl 



\s + ik l U {z)\ 2 



avec s = a — ioj. On voit que si Uq(z) ne change pas de signe, c'est-a-dire si 
Uq(z) n'admet pas de point d'inflexion, l'integrale de cette equation est non 
nulle. On a done a = 0, ce qui interdit toute instabilite. 

L'existence d'un point d'inflexion dans le profil Uq(z) est une condition necessaire 
mais non suffisante pour l'existence d'une instabilite. En effet, il faut non 
seulement que a soit non nul, mais aussi positif pour certaines valeurs de k\ 
si Ton doit observer des modes instables. 



4.3 Instabilities visqueuses 

Dans le cas general v ^ 0, l'equation d'Orr-Sommerfeld permet de trouver 
les valeurs propres s G W et, eventuellement, le spectre continu s G V U W 
pour un vecteur d'onde k\ donnee. On rappelle que l'on considere la stabilite 
de profils Uq{z) de formes quelconques, meme si, en principe, seules les pro- 
fils paraboliques sont des solutions exactes des equations de Navier-Stokes. 
Cette approche est rendue possible en negligeant le terme de viscosite pour 
revolution du mouvement moyen cette viscosite etant neanmoins presente 
dans l'equation de Orr-Sommerfeld. 

Nous nous contentons ici d'indiquer schematiquement l'allure des frontieres de 
stabilite dans un plan (^,k±) pour differentes formes de profils Uq(z). L'etude 
approfondie des solutions de l'equation d'Orr-Sommerfeld est un probleme 
assez difficile qui fait l'objet de plusieurs monographies citees en reference. 



STABILITE DES ECOULEMENTS PARALLELES 
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Figure 1.8: Seuil d'instabilite pour un profil avec point d'inflexion 
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Figure 1.9: Seuil d'instabilite pour un profil sans point d'inEexion 
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La figure 1.8 presente Failure schematique de la frontiere entre les regimes 
stables et instable dans le plan (^, k\) dans le cas ou le profil Uq(z) admet un 
point d'inflexion. On voit que dans la limite v — ► zero il existe un intervalle 
de vecteurs d'ondes k\ pour lesquels certains modes sont instables. 
La figure 1.9 represente schematiquement le cas ou Uq(z) n'adment pas de 
point d'inflexion. Comme predit par le theoreme du point d'inflexion, l'ecoulement 
est stable dans la limite v — > 0. C'est done la viscosite qui est responsable 
de Pinstabilite observee pour des v non negligeables. Le mecanisme physique 
responsable de cette instabilite visqueuse est assez subtil et n'est pas explicite 
ici. En gros, la presence de viscosite induit un dephasage dans la direction 
x des composantes periodiques du champ de vitesse, ce qui conduit a des 
mecanismes amplificateurs de la perturbation. 

Pour conclure, on peut evoquer le cas de l'ecoulement de Couette qui est 
caracterise par un profil lineaire, e'est-a-dire a la frontiere entres les profils 
avec ou sans points d'inflexions. L'etude de stabilite lineaire montre qu'il 
n'existe aucun mode instable pour cet ecoulement. Les instabilites observees 
dans la realite sont due a des perturbations d'amplitudes finies, et done a des 
effets non-lineaires, ce qui sort du cadre de la presente etude. 

Conclusion 

Les quatre exemples qui ont ete aborde ont permis de detailler le calcul de 
la relation de dispersion generalised de modeles lineaires ou linearise autour 
d'un equilibre. Les modeles prototypes ID ont permis de ramener la recherche 
des branches de valeurs propres, parametrees par le vecteur d'onde k±, a la 
recherche des racines d'un polynome. L 'instabilite des rolls- waves pour les 
equations de Saint Venant, qui se declenche pour un nombre de Froude egale 
a deux, a permis de voir les limitations du modele lorsque la dissipation 
visqueuse des petites echelles est negligee et que celles-ci sont les plus instables. 
L'etude de la stabilite des ecoulements paralleles a permis de montrer que 
l'existence d'un point d'inflexion dans le profil influence fortement la nature 
des instabilites. On peut voir Pinstabilite de Kevlin-Helmoltz comme un cas 
limite d'un ecoulement parallele admettant un point d'inflexion. 
Pour les quatre exemples, on remarque que les modeles sont invariants par 
translations en temps et en espace dans la direction x. Ces invariances ex- 
pliquent pourquoi Ton a recherche les modes propres a l'aide de fonction 
exp(i k\ x + st). Lorsque le systeme n'est pas invariant par translation en z, 
comme c'est le cas pour Pinstabilite des ecoulements paralleles, la forme des 
modes propres doit etre trouvee en appliquant des conditions aux limites. 
On aurait pu aussi chercher des solutions sous la forme exp(7x + i k\ x + 
iuot), ou encore exp(7x + ik\x + at + iut) pour prendre en compte des 
amplifications en espace. Cette approche de la theorie des instabilites n'a pas 
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ete abordee ici, mais fait partie des outils important pour aller plus loin dans 
la comprehension des phenomenes realistes en particulier pour les ecoulements 
ouverts, comme par exemple les ecoulements paralleles. 



FORMULAIRE 



MODELES ID 

Modele lineaire 

du du d 3 u d 2 u d 4 u 

a = S(fci) = hq — fj,2 k\ + ^4 kf et uj = Q,(k{) = a\k\ — «3 k\ 
Linearisation autour d'un equilibre 

du du o d 2 u d A u 



du du „ d 2 u d 4 u 



Calcul de stabite 



dU _ f d d d_ \ ~ 
dt dt = \dxi' dx 2 : dx^J ~ 

U(x,t) = Ae l ^ +St =► s A = Uiki,ik 2 ,ik 3 ) A 
=>■ Oi = Sj(fc) et uji = £li(k) pour i = 1, ...,n 
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FORMULAIRE 



INSTABILITY "ROLL WAVES" 




Modele de Saint Venant 




dU TT dU .dh T 1 „ U\U\ dh 

m +u a, = - s a x - s, -2 C < h et at + 




Linearisation autour de (h n , U n ) : 




dt n dx 3 dx f \h n 


*K ) 


dh TT dh dU 
dt dx dx 




Equations adimensionnees 




dU. dU* _ 1 dh. 1 ( V*\ dh 
dt +U * dx ~ dx + 2 Cf { 1 K) dt 4 


d 

W (TT U \ — ft 

- -q^ \<~>* At* J — U , 


Calcul de stabite 




(S + iK + l)(S + iK) + ^if 2 + \ iK = 





Instable pour tout K des que F > 2. 
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INSTABILITE DE KELVIN HELMOLTZ 
Equations d'Euler avec interface 

div U j = et pi ~^Li = —grad p% — pi g e z pour i = 1, 2 



Conditions aux limites a l'interface 
drj 



dt 



+ U_i ■ grad rj = w 1 , pi = p 2 , 



at 



+ J7 2 • grad rj = w 2 



Linearisation autour de l'equilibre 



Hi = Ui x + grad fa et A<^ = 
avec les conditions aux limites 



Pi 



8 TT d 



P2 



9 TT 9 , 



en z = 



Calcul de stabite 

pi gk + (s + i k\ Ui) 2 



P2 



gk-(s + ikiU 2 f 



Cas pi = p 2 : (s + « fei ^i) 2 + (s + i fci f/ 2 ) 2 = 
Cas U! = U 2 = 0: s 2 = p -^± g k 

1 Z P1+P2 3 



STABILITE DES ECOULEMENTS PARALLELES 
Equation d'Orr-Sommerfeld 

s + i h U {z) -u(D 2 - k 2 )] (D 2 - kj) w m = U^z) w m . 

Theoreme du point d'inflexion 

Une condition necessaire pour que Uq{z) soit instable est l'existence d'un point 
d'inflexion. 

Instabilites visqueuses 





1/V 




1/V 
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EXERCICES ET PROBLEMES 

EXERCICE 1.1 Obstacle immerge et roll-waves 




Figure 1.10: Ecoulement a surface libre sur un plan incline 

On considere un ecoulement peu profond a surface libre sur un plan incline 
faisant un angle a avec l'horizontale dans un champ de gravite d'intensite g. 
On suppose que cet ecoulement est unidimensionnel (ID) et on note U(x, t) la 
vitesse moyennee dans un plan tranverse a l'ecoulement et h(x, t) la hauteur 
de la surface libre dans ce plan. On note I = sin a et g' = g cos a. 

Modele d'onde cinematique 

On considere tout d'abord le modele dans lequel la vitesse U est reliee a la 
hauteur h par la relation U = K \fl ha ou K est une constante. 

1) Ecrire la loi de conservation de la masse sous la forme d'une equation aux 
derivees partielles ne faisant intervenir que la hauteur h(x, t). 

2) On s'interesse aux petites perturbations h(x, t) autour de la solution sta- 
tionnaire h n en posant h = h n + h. Ecrire l'equation linearisee regissant 
la dynamique de h. 

3) Justifier la recherche de solutions complexes sous la forme : 
h = B exp(i k\ x — ito t) ou B G W. 

4) Tracer la relation de dispersion en adoptant la convention co > 0. 

5) Comparer les vitesses de phase et de groupe des ondes a la vitesse U n de 
l'ecoulement. 

6) En x = 0, on fait osciller a la pulsation ljq un obstacle imerge dans 
l'ecoulement. Decrire et dessiner le train d'ondes resultant de cette oscil- 
lation dans le cadre du modele propose. 

7) A quelle vitesse se propagent les ondes emises par l'obstacle ? Comparer 
avec la vitesse des particules fluides. 
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Modele de Saint-Venant 

On considere maintenant le modele des equations de Saint-Venant dans l'appro- 
ximation des pentes faibles qui s'ecrit 

dU TT dU , dh dh d .„ , , . . 

8) On s'interesse aux petites perturbations U(x,t) et h(x,t) autour de la 
solution stationnaire (U n , h n ) G ]R+ en posant U = U n + U et h = h n + h. 
Ecrire les equations linearisees regissant la dynamique de U et h. 

9) Justifier la recherche de solutions complexes sous la forme : 
(U,h) = (A, B) exp(ik lX -iujt) ou (A, B) G d 2 . 

10) Tracer les relations de dispersion des ondes en adoptant la convention 
uj > 0. Discuter le trace en fonction du nombre de Froude F = U n / V 'g h n . 

11) En a; = 0, on fait osciller a la pulsation loq un obstacle imerge dans 
l'ecoulement. Decrire et dessiner les trains d'ondes resultant de cette 
oscillation dans le cadre du modele propose. 

Modele des "roll-waves" 

On considere maintenant le modele des equations de Saint-Venant avec frot- 
tement turbulent qui s'ecrit 

dU TT dU , dh r \U\ 

f t+ liUh) = 0, (1.63) 

et ou suppose que le coefficient de frottement est modelise par une loi de la 
forme C f {h) = D h~P avec D > et (3 > 0. 

12) On s'interesse aux petites perturbations U(x, t) et h(x, t) autour de la 
solution stationnaire (U n , h n ) G M+ en posant U = U n + U et h = h n + h. 
Ecrire les equations linearisees regissant la dynamique de U et h. 

13) Justifier la recherche de solutions complexes sous la forme : 

(U, h) = (A, B) exp(i k± x + s t) avec s = a - i u £ W et (A, B) G d 2 . 

14) Ecrire, par exemple a l'aide de grandeurs adimensionnees, la relation 
de dispersion generalisee regissant la stabilite de la solution stationnaire 
(U n , h n ). On pourra noter d = D h~@ et F = U n /V9 h n - 

15) On observe l'apparition du seuil d'instabilite lorsque F > 1.5. En deduire 
la valeur de l'exposant [3. 

16) On regie U n et h n pour se placer dans le cas marginal F = 1.5 et l'on 
fait osciller, en x = et a la pulsation ujq un obstacle imerge dans 
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l'ecoulement. Decrire et dessiner les trains d'ondes resultant de cette 
oscillation dans le cadre du modele propose. 

Corrige page 48 



CORRIGES DES EXERCICES ET PROBLEMES 



Corrige 1.1 Obstacle immerge et roll- waves 



Modele d'onde cinematique 

l)La loi de conservation de la masse s'ecrit ^ + ^U(h) || = oil U depend de 

h comme indique. 2)L'equation linearise autour de h n s'ecrit ^ + | U n || = 
0, avec U n = U(h n ). 3)L'equation lineaire etant invariante par transla- 
tion dans l'espace et le temps, on peut chercher des solutions sous forme 
d'exponentielles. On anticipe qu'il n'y aura ni amplification ni amortissement 
des ondes (le calcul le confirmerait). 4)La relation de dispersion uj = Q(ki) = 
\U n k\ n'est definie que pour k\ > avec la convention u> > 0. Son trace est 
done forme d'une demi-droite. 5) On a, = c g = | U n . Le milieu n'est pas 
dispersif. 6)Seule l'onde verifiant k\ = ko avec ko = ^u>o/U n est emise dans le 
sillage lointain. Ce sillage est uniquement observe dans la region x > dans 
la mesure ou la vitesse de groupe est positive. 7)Ces ondes, situees en aval 
de l'ecoulement, se propagent a une vitesse 5 U n /3 qui est done plus grande 
que la vitesse U n de l'ecoulement moyen. 



Modele de Saint- Venant 



8)Les equations linearisees s'ecrivent ^ + U n ^ + g' ' = et ^ + h n % + 

= 0. 9)Les equations sont invariantes par translation en espace et en 

, in x T ,v ., f-iu + Unifa g'iki \ ( A\ 

temps. lOlLe systeme s ecnt ... . TT . , ,-, = U. 

\ ik\n n —iuj + U n ik n J\BJ 

En annulant le determinant, on obtient la relation de dispersion (uj — kiU n ) 2 = 
g' h n k\. Avec la convention uj > et en posant c n = yV h n , on obtient les 
relations uj = (U n — c n )ki et uj = (U n + c n )k\ respectivement definies pour 
k\ < et k\ > lorsque F < 1 et pour k\ > seulement pour F > 1. Les 
traces de ces relations de dispersion sont des demi-droites. ll)Pour F < 1, 
on observe une onde de vecteur d'onde k\ = LJo/(U n + c n ) se propageant vers 
la droite a une vitesse c v = c g = U n + c n , et une onde de vecteur d'onde 
k\ = uJo/(U n — Cn) qui remonte le courant a la vitesse U n — c n < 0. Pour 
F > 1, ces deux ondes se propagent vers l'aval. 
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Modele des "roll-waves" 

12)Lej3ysteme s'ecrit jf + U n § + <?' f + D hn {l+P) U n - ±±£ D hn {2+/3) U 2 = 

® e * W + fi" + U n || = 0. 13)L'invariance par translations et le souhait de 
borner les solutions a l'infin conduit a cette forme de solutions element aires. 
14)En posant x = h n x*, t = jj^U, U(x,t) = U n U*(x*,t*) et h(x,t) = 

K h*(x*,U), on obtient le systeme g£ + £/* |g = §^ + \ d (l - 

et ^ + (U* h*) = 0. La relation de dispersion generalised s 'obtient en 
annulant le determinant du syteme linearise, ce qui s'ecrit 



5* ~\~ i k* ~\~ d 2 ^ 

i k* "I - ^ k% 



(s*+z&*+eZ)(s*+z £;*)-«/(;*( % - = . 



Cette relation de dispersion s'ecrit (S+iK+l)(S+iK) + j^K 2 + ^-iK = 
en posant S = s/d et K = k\jd. 15)Au seuil de l'instabilite, on a S = — ifi, 
ce qui entraine — i^ 2 ) 2 = K 2 /F 2 et£l — K = X en separant les parties 
reelles et imaginaires de l'equation. On en deduit que le seuil est obtenu pour 
F = j^j. Pour un seuil a F = 3/2, on obtient (3 = 1/3. On retrouve la 
parametrisation de Strickler. 16)Au seuil, la relation de dispersion s'ecrit 
= K, ce qui conduit, dans le cas (3 = 1/3, a!J = | K. En repassant 
aux unites dimensionnelles, on obtient la relation de dispersion uj = | U n k\ 
au seuil de l'instabilite. L'onde marginale emise par un obstacle est observee 
uniquement en aval de l'obstacle et a pour longueur d'onde k\ = ^u!o/U n . Sa 
vitesse et de phase et sa vitesse de groupe sont egale a | U n . L' autre onde 
emise par l'obstacle est amortie. 

QUESTIONNAIRES A CHOIX MULTIPLES 
CORRIGES DES QCM 
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OBJECTIFS PEDAGOGIQUES 

L'objectif de cet article pedagogique est de presenter la notion d'instabilite 
sur Pexemple de la convection de Rayleigh-Benard. Plusieurs objectifs de 
formation sont vises : 



• Se familiariser avec les equations de la mecanique des fluides dans le 
cadre de 1' approximation de Boussinesq. 

• Maitriser la formulation des differents jeux de conditions aux limites qui 
conditionnent le forgage thermique et l'effet des parois sur le champ de 
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NOTATIONS 



vitesse. 



• Maitriser et comparer plusieurs adimensionnalisation des equations d'un 
mo dele de mecanique des fluides. 

• Se familiariser avec la notion d'elimination de la pression dans le cas 
bidimensionnel avec l'introduction d'une fonction de courant. 

• Assimiler la notion de linearisation d'un modele autours d'une solution 
stationnaire. 

• Assimiler progressivement les etapes du calcul de la relation de disper- 
sion generalised autour du profil conductif. 

• Determiner les sueils critiques de l'instabilite convective a partir de 
l'expression des branches de valeurs propres. 

Les competence a acquerir lors de l'etude de cet article pedagogique sont les 
suivantes : 

• Etre capable de lineariser un systeme d'equations aux derivees partielles 
autour d'une solution stationnaire. 

• Etre capable d'adimensionner un systeme d'equations aux derivees par- 
tielles en effectuant un choix d'unites. 

• Maitriser le choix de la forme des vecteurs propres en fonction des in- 
variances par translation du probleme. 

• Maitriser le calcul de la relation de dispersion generaliser pour en deduire 
les seuils d'instabilite. 



Le niveau requis pour la lecture de cette article pedagogique se situe autour 
de celui d'une Licence scientifique. Quelques notions sur la formulation des 
equations de Navier-Stokes peuvent se reveler utiles. 



PRE-REQUIS 



NOTATIONS 



b= (bi, 62, •••) Vecteur de composantes bi 
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Derivee partielle par rapport a x 
Operateur ^ 

Tenseur des taux de deformation (s _1 ) 
Produit doublement contracte = tr (DD) ( s_2 ) 
Energie interne massique (J kg -1 ) 
Loi d'etat pour l'energie interne (J kg -1 ) 
Vecteurs unitaires de la base canonique (m) 
Densite volumique des force exterieures (N m~ 3 ) 
Gravite 

grad Operateur gradient m _1 

J(f,g) Operateur % g - % g (m^) 

k = (ki,k2,ks) Vecteur d'onde (m 1 ) 

k Norme du vecteur d'onde k (m _1 ) 

ko Nombre d'onde 2 n/d (m ~ 1 ) 

k c Nombre d'onde critique (m _1 ) 

L Periodicite horizontale (m) 

P Nombre de Prandtl 

p(x, t) Champ de pression (Pa) 

V Loi d'etat pour la pression (Pa) 

Po(z) Profil de pression (Pa) 

p r Pression constante (Pa) 

q Flux de chaleur (W mT 2 ) 

r Taux de production volumique de chaleur (J kg -1 

R Nombre de Rayleigh 

R n (ki) Nombre de Rayleigh critique 

R c {ki) Nombre de Rayleigh critique R\(ki) 

1Z Loi d'etat pour la masse volumique (kg m~ 3 ) 

i?* Valeur particuliere de R 

R c (ki,n) Rayleigh critique pour le mode n et pour k\ 

R c {k\) Rayleigh critique R c (k\, 1) 

Ri(k\) Rayleigh critique pour les modes impaires 

R p {k\ Rayleigh critique pour les modes paires 

M Droite des reels 

s Valeur propre complexe (s _1 ) 

T(x, t) Champ de temperature (°K) 

T\,T2 Temperatures constantes (°K) 

T c (z) Profil de temperature (°K) 

T r Temperature de reference (°K) 

t Temps (s) 

l[(x,t) Champ de vitesse (m s -1 ) 

U_ rn Vecteur amplitude complexe de vitesse (m s^ 1 ) 

(u, v, w) Composantes de U_(x, y, z) (m s _1 ) 
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V 


Vitesse dans la direction e y (m s _1 ) 


V m (z) 


Profil complexe (m s _1 ) 


W 


Vitesse dans la direction e z (m s _1 ) 


w m (z) 


Profil complexe (m s _1 ) 


X 


Coordonnee sur l'axe e x (m) 


x = (x,y,z) 


Coordonnees spatiales (m) 


y 


Coordonnee sur l'axe e y (m) 


z 


Coordonnee sur l'axe e z (m) 


a 


Coefficient d'expansion thermique (°K _1 ) 


P 


Rapport d/L 


r 


Gradient de temperature °K m _1 ) 


div 


Operateur divergence m _1 


A 


Operateur Laplacien m -2 


c 


Nombre sans dimension egal a 1 ou k/v 




Nombre sans dimension egal a 1 


0(x,t) 


Ecart de temperature (°K) 




Amplitude complexe de temperature (°K) 


[@] 


Unite de temperature (°K) 


@(z) 


Profil de temperature (°K) 


K 


Coefficient de conductivite thermique (m 2 s _1 ) 


A 


Constante (m _1 ) 


An 


Coefficient de viscosite de Lame (Pa s) 


Mn 


Coefficient de viscosite de Lame (Pa s) 


V 


Viscosite cinematique (m 2 s _1 ) 


n 


Pression generalisee (Pa) 


p(x, t) 


Champ de masse volumique (kg m~ 3 ) 


Pr 


Masse volumique de reference (kg m~ 3 ) 


ip(x,z) 


Fonction de courant (m 2 s _1 ) 




Amplitude comple de fonction de courant (m 2 s _1 ) 




Profil de fonction de courant (m 2 s _1 ) 



COURS ECRIT 
Introduction 

L'objectif est d'enoncer ici les equations du modele permettant de rendre 
compte des mouvements de convection rencontres dans un fluide soumis a un 
gradient de temperature puis de calculer les seuils de declenchement de cette 
instabilite. 

L' approximation de Boussinesq permet une formulation incompressible des 
equations de Navier-Stokes en prenant en compte des forces de flottabilite 
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Figure 2.1: Rouleaux de convection. 



(poussee d'Archimede) dues a la dilatation du fluide et induite par une vari- 
ation de la temperature. 

Le systeme d'equations de l'approximation de Boussinesq etant enonce, on 
s'interesse a la configuration idealisee d'une couche fluide horizontale forcee 
par un gradient thermique vertical, soit en imposant les temperatures aux 
frontieres, soit en imposant un flux thermique. On explore alors les differents 
choix d'adimensionnalisation permettant la formulation mathematique du 
probleme de Rayleigh-Benard. 

On envisage alors le cas ou les temperatures de plaques sont fixees ainsi que le 
cas ou le flux de chaleur est fixe. De meme, on compare le cas des conditions 
aux limites cinematiques libres et rigides. L'etat conductif est caracterise par 
une vitesse nulle et un profil lineaire de temperature a gradient vertical. On 
suppose que le bas de la couche est plus chaud que le haut et Ton augmente 
progressivement le forcage. Lorsque le forcage thermique est faible l'etat con- 
ductif est stable. Une perturbation de temperature ou de vitesse est dissipee 
par la diffusion thermique ou la viscosite avant que la force d'Archimede n'aie 
eu le temps d'agir. Lorsque le forcage thermique depasse un seuil critique, 
l'etat conductif devient instable et des mouvements de convection apparais- 
sent. 

Cette etude de stabilite s'effectue en linearisant les equations autour de l'etat 
conductif. On montre alors que la stabilite de l'etat conductif ne depend que 
de la valeur d'un nombre sans dimensions appele nombre de Rayleigh. On 
calcule ici le nombre de Rayleigh critique pour plusieurs configurations de 
conditions aux limites. 



1 Approximation de Boussinesq 

Pour decrire le mouvement de convection d'un fluide chauffe par le bas lorsque 
le champ de vitesse reste faible devant la vitesse du son, les equations de 
Navier-Stokes compressibles sont trop complexes et les equations de Navier- 
Stokes incompressibles ne rendent pas compte des forces de flotabilite. L'appro- 
ximation de Boussinesq permet de prendre en compte ces forces tout en filtrant 
les ondes sonores. 
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1.1 Fluides compressibles ou incompressibles 

On rappelle tout d'abord les equations de Navier-Stokes incompressibles en 
utilisant les notations de l'ouvrage "Introduction a la Mecanique des milieux 
continus deformables" (O. Thual, Cepadues 1997) : 



ou A n et p n sont les deux coefficients de viscosite de Lame (avec A„ = —2 // n /3 
selon l'hypothese de Stokes) et k le coefficient de conductivity thermique. Les 
notations suivantes sont utilisees pour les grandeurs physiques : p pour le 
champ de densite volumique, U pour le champ de vitesse, p pour le champ de 
pression, e pour le champ d'energie interne volumique, T pour le champ de 
temperature et D pour le champ de tenseur des taux de deformations obtenu 
en prenant la partie symetrique du gradient du champ de vitesse. La notation 
^j. = §t + U_- grad designe la derivee particulaire. 

Les deux lois d'etat £ et V dependent du fluide considere. Ici, la densite 
volumique des forces exterieures / = — p g e z est celle induite par le champ 
de gravite. On a suppose qu'il n'y avait pas de taux de production volumique 
de chaleur (r = 0). Pour fermer ce systeme d'equations, il faut specifier 
des conditions aux limites, et, dans le cas ou Ton s'interesse aux solutions 
instationnaires, des conditions initiales. 

L' approximation de fluide incompressible consiste a imposer p = p r partout 
et a neutraliser la loi d'etat p = V(p, e) qui n'a alors plus de sens. Ceci revient 
a se placer dans la limite ou la vitesse du son est infmiment grande devant la 
vitesse de l'ecoulement. La pression, qui etait une grandeur thermodynamique 
dans le cas compressible, devient une grandeur purement dynamique permet- 
tant de satisfaire la contrainte div [7 = (on dit que p est le multiplicateur de 
Lagrange associe a cette contrainte). On dit aussi que l'approximation de flu- 
ide incompressible permet de "filtrer" les ondes sonores. Ce filtrage permet 
d'augmenter considerablement le pas de temps des simulations numeriques 
qui devrait, sinon, etre plus petit que la periode des vibrations sonores. 

1.2 Equations de l'approximation de Boussinesq 

Pour certains ecoulements, l'approximation de fluide incompressible est trop 
drastique, car elle ne permet pas de prendre en compte les variations de 
densite dues aux variations de temperature. En presence de gravite, des 
inhomogeneftes de densite induisent un champ de force d'Archimede qui peut 



dp 



—p div U_ 




-grad p - p g e z + (A n + p n ) grad div U + p n AU 

k AT - p div U + A n (div U) 2 + 2 p n J2 : J2 
e = £{p,T) et p = P(p, e) 



(2.1) 
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mettre le fluide en mouvement. Ces mouvements ne sont pourtant pas de 
meme nature que les ondes sonores que l'on souhaite continuer a "filtrer" . 
L 'approximation de Boussinesq permet de prendre en compte ces variations 
de densite tout en supposant que la vitesse des ondes sonores est infinie devant 
la vitesse de Pecoulement. Elle consiste a ecrire les equations devolution du 
fluide sous la forme 

div U = 

-grad p- pge z + p n AU 
k AT + 2 fi n J2: J2 

£(p r ,T) et p = p r [l-a(T-T r )]. (2.2) 



Pr 
Pr 



dU 

~dt 
de 

e = 



On voit que Ton a suppose p egal a une densite constante p r partout, sauf 
dans le terme p g e z qui est responsable de la force d'Archimede. Cette 
hypothese d'incompressibilite a done entraine div U = 0, supprimant de ce fait 
plusieurs termes de l'equation de conservation de la quantite de mouvement 
et de l'equation de bilan de l'energie interne. 

D'autre part, cette hypothese d'incompressibilite interdit de prendre en compte 
une loi d'etat de la forme p = V(p, e) et la pression p devient un parametre 
dynamique decouple de la thermodynamique. Cette loi d'etat est alors rem- 
placee par une loi p = 1Z{T) que l'on suppose lineaire en ecrivant p = 
1Z(T) = p r [l — a(T — T r )]. Le coefficient a est appele "coefficient de dila- 
tion thermique" . II est positif dans la mesure ou un echauffement autour de 
la temperature T r induit une diminution de la densite. 



La loi d'etat e = £(p r ,T) permet d'ecrire = C v (p r ,T)^ avec C v = 




En reportant les lois d'etat dans le systeme d'equation, celui-ci se met sous 
la forme 



div U = 



dU 

~dt 
dT 



= -grad ( — + g z ) + a g (T - T r )e z + v AU 

\Pr / 



V 



k AT + 2 — J2_ : J2 



(2.3) 



ou v = p n /pr est la viscosite cinematique, n = k/(p r C v ) est le coefficient de 
diffusivite thermique a volume constant. 

Dans la plupart des applications, on neglige le terme 2 {v/C v ) D '■ D, qui 
correspond a l' echauffement du fluide par les forces visqueuses, devant le terme 
k AT, qui correspond a la diffusion de la temperature dans le fluide. 
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2 Probleme de Rayleigh-Benard 

L' approximation de Boussinesq ayant ete explicitee, il est maintenat possible 
de construire le modele rendant compte de l'ecoulement de convection observe 
lorsque Ton chauffe un fluide par le bas. Ce forcage est traduit dans la modele 
par la formulation de conditions aux limites. 

2.1 Conditions aux limites 

On considere un fluide compris entre deux plaques planes horizontales delimi- 
tant les frontieres d'equations z = et z = d. L'epaisseur de la couche fluide 
est done la distance d. 




Figure 2.2: Couche fluide chauSee par le bas. a) Schema des rouleaux de 
convection, b) Visualisation experimental du champ de temperature. 

Pour resoudre les equations de Navier-Stokes dans 1' approximation de Boussi- 
nesq, il est necessaire de specifier des conditions aux limites pour la tempe- 
rature et pour la vitesse. Pour ces deux categories de conditions aux limites, 
plusieurs variantes sont possibles suivant la nature du probleme physique con- 
sidere. Nous enumerons ici les principales configurations. 
On considere deux types de conditions aux limites en temperature : 

Temperatures fixees : T = T\ en z = 0, T = T<i en z = d. 

Flux fixe : =q en z = et z = d. (2.4) 

Les conditions aux limites en temperatures fixees correspondent au cas ou les 
parois planes sont infmiment conductrices et relient le fluide a des thermostats 
de temperatures fixees. Les conditions aux limites en flux fixe correspondent 
a un chauffage delivrant une quantite de chaleur determinee. Lorsque q < 
ou T\ < T2 le fluide est dit stratifie. Le probleme de Rayleigh-Benard, que 
Ton considere ici, s'interesse aux cas ou q > ou T% > T2 pour lesquels des 
mouvements de convection peuvent se developper. 
On considere deux types de conditions aux limites en vitesses : 

Libres : gj* = gj! = 0etu; = en z = et z = d. 
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Rigides : u = v = w = en z = et z = i (2.5) 

Les conditions aux limites libres correspondent au cas idealise oil la couche 
fluide considered est comprise entre deux autres fluides et oil la deformation 
des interfaces est negligee. Les conditions aux limites rigides correspondent 
au cas ou les frontieres sont des parois rigides. 

2.2 Etat conductif 

On appelle etat conductif la solution stationnaire telle que U = et T = 
T c (z) = Ti — Tz. Cette solution satisfait les conditions aux limites en tempe- 
ratures fixees en choisissant T = {T\ — T2)/d ou bien les conditions aux limites 
a flux fixe en choisissant T = q/n. 

On definit 9, l'ecart de temperature au profil conductif, a l'aide du changement 
de notation T(x,t) = T c (z) + 9(x,t). Les equations du mouvement s'ecrivent 
alors 



dU 

~dt 
d9 



div £/ = 

— grad II + ag9e z + vMJ_ 
T w + k A9 (2.6) 



dt 

avec LT = p/p r + g z+a gT r z—a gT\ z+\a gT z 2 +Ii r ou LT r est une constante 
arbitraire. Dans ces equations, on a neglige le terme 2 (y/C v ) D:D,l& vitesse 
verticale est notee w et Ton a utilise T c (z)e z = grad (T\ z — \ T z 2 ). Le terme 
T w represente l'advection de temperature par la vitesse due au gradient 
vertical de temperature T. 

2.3 Modele sous forme dimensionnelle 

En explicitant les derivees particulaires, les equations du mouvement s'ecrivent 

diviZ = 

dU 

— — h U ■ grad U_ = —grad H + ag9e z + v AU_ 
at 
r)9 

— + [/■ grad0 = Tw + kA9. (2.7) 

Les deux sortes de conditions aux limites pour la temperature considerees ici 
sont 



Temperatures fixees : 9 = en z = et z = d. 

.de 

"dz 

Les deux sortes de conditions aux limites en vitesses sont : 



Flux fixe : k3§ = en z = et z = d. (2.8) 



Libres : || = ||=0etu> = en z = et z = d. 

Rigides : u = v = w = en z = et z = d. (2.9) 
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3 Equations sous formes adimensionnees 

Le choix d'un systeme d'unite permet d'ecrire les equations du modele sous 
forme adimensionnee. Cette procedure est essentielle pour determiner le nom- 
bre de parametres sans dimension qui controlent le probleme. Toutes les adi- 
menssionnalisations sont equivalente tant que Ton ne neglige pas de termes. 
On presente ici deux adimensionnalisation qui different par le choix de l'unite 
de temps. 

3.1 Choix des unites 

Nous allons considerer plusieurs systemes d'adimensionnalisation en explo- 
rant plusieurs choix d'unites. L'unite de longueur sera toujours basee sur 
l'epaisseur d de la couche fiuide, ce que Ton note [L] = d. On note [r] = ( d 2 / n 
le choix des unites de temps, le facteur ( etant un nombre sans dimension que 
l'on pourra choisir ulterieurement. De meme, on note [0] = £ d T le choix des 
unites de temperature, le nombre sans dimension £ restant a choisir. On a sup- 
pose ici T > 0. On peut done resumer cette famille d'adimensionnalisations 
en ecrivant 

[L] =d [t] = ( d 2 / K et [6] = f d V . (2.10) 

Si Ton choisit, par exemple, £ = 1, l'unite de temps [r] = d 2 / k est le temps 
caracteristique de diffusion de la temperature. Si l'on choisit maintenant 
£ = k/u, l'unite de temps [r] = d 2 jv est le temps visqueux. Si l'on choisit 
£ = 1, l'unite de temperature [O] = T d = T\ — Ti est la difference de 
temperature entre le bas et le haut de la couche pour l'etat conductif considere. 
Quelque soit le choix de ( et £, l'unite de vitesse est le rapport entre l'unite 
de longueur et de temps ce que l'on ecrit 

[U] = [L\/[r] = \- d - (2.11) 

3.2 Equations adimensionnees 

En utilisant les memes notations pour les champs que dans le cas oil ces 
grandeurs sont dimensionnees, les equations adimensionnees s'ecrivent 



-grad U + ^RPde z + PAU 
lw + A8 (2.12) 



dU 

86 
di 



dwU = 
+ LL • g?ad U ) = 



+ U. ■ grad 8 
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ou apparaissent les deux nombres adimensionnes que sont P, le nombre de 
Prandtl, et R, le nombre de Rayleigh, dont les expressions sont les suivantes 

p= v - et fl = ^r = ^(r 1 -r 2 ) j 

K UK UK 

L'ecriture des conditions aux limites en z = et z = 1 (en coordonnees 
adimensionnees) est identique a la formulation dimensionnee. 

3.3 Adimensionnalisation visqueuse 

On choisit ici £ = 1/P et £ = 1 ce qui conduit au systeme d'unites [L] = d, 
[t] = d 2 /u, [U] = u/d et [©] = T d = T\ — Les equations adimensionnees 
s'ecrivent 

diviZ = 
^= + U ■ grad U = -grad ^ + R 6 e z + AU 

Q t — O — p — 

P^ + C/-grad^ = u; + A^ (2.14) 

3.4 Adimensionnalisation thermique 

On choisit ici Q = ^ = 1 ce qui conduit au systeme d'unites [L] = d, [r] = d 2 /K, 
[U] = n/d et [O] = T d = Ti — T 2 . Les equations adimensionnees s'ecrivent 

divt/ = 

-= + U ■ grad U = -grad U + RP6e z + PAU 

dt — ~ 

— grad0 = w + A^. (2.15) 

4 Probleme de Rayleigh-Benard bidimensionnel 

Par simplicity, on se restreint ici aux ecoulements bidimensionnels dans un 
plan vertical. Un raison physique permettant de justifier cette approche est le 
fait que Ton observe des rouleaux de convections 2D au seuil de Pinstabilite, 
du moins pour les conditions aux limites presentees ici. 



4.1 Fonction de courant 

On se place ici dans le cas 2D defini par le fait que U(x, z) et 9(x, z) ne 
dependent pas de la coordonnee y. Dans le cadre de Padimensionnalisation 
thermique, les equations du mouvement s'ecrivent alors 
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du du du dU „ . 

— + u — + w— = -— + PAu 
at ox oz ox 

dv dv dv dll _ . 

— + u — = -— + PAv 
at ox Oz dy 

dw dw dw dll „ „ „ „ . 

at ax az 
^ d9 d9 

dt +U dx +W dz = W + A °- ^ 

On voit que l'equation d'evolution de v est decouplee des equations d'evolution 
de u, w et 6. La composante v se comporte done comme un scalaire passif. 
En l'absence d'un gradient de pression dans la direction y, cette vitesse tend 
done vers zero. 

L'equation div U_ = |^ + ^ = entraine alors qu'il existe une fonction de 
courant tp(x,z) telle que u = et w = Les conditions aux limites en 
vitesses s'expriment alors de la maniere suivante 

Libres : = et g = en z = et z = 1. 

Rigides : §f = et J£ = en z = et z = 1. (2.17) 



4.2 Elimination de la pression 

Pour eliminer la pression generalised LT, il suffit de prendre le rotationnel 
de l'equation de quantite de mouvement. En remarquant que rot U_ = 
—A<p e^ 2 ) = — + )e^ 2 ^ et en projetant sur la direction e^ 2 ** orthog- 
onale au plan 2D, on obtient les equations du mouvement 

d df) 

— AV> + J(^,AV0 = RP-^ + PA 2 ^ 

oil la notation J est definie par J(f,g) = %. ff - |f 



5 Conditions aux limites periodiques 

Le calcul de stabilite est simple s'il on suppose que Pecoulement est periodique 
dans toutes les directions. En effet, ces symetries de translations permettent 
de considerer des modes propres ayant la forme d'exponentielles en espace 
et en temps, solution du modele linearise autour de la solution stationnaire 
conductive. Ce cas est presente ici pour permette une approche moins aride 
de la notation d'instabilite. Les calculs seront utiles pour des conditions aux 
limites plus realistes. 
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5.1 Linearisation 

On linearise les equations de Navier-Stokes dans 1' approximation de Boussi- 
nesq autour de l'etat conductif decrit par U = 0, T = T c (z) = T\ + T z 
et II = en choissisant la constante arbitraire II r convenab lenient. Avec la 
notation T = T c {z) + 9, cet etat est caracterise par 9 = 0. 
En choisissant les unites de longueur [L] = d, de temps [r] = d 2 / k et de 
temperature [0] = d T = T\ — T2 de Padimensionnalisation thermique, les 
equations linearisees s'ecrivent 

div U = 

-grad Il + RP6e z + PAU 
Tw + A6. (2.19) 

L 'expression des conditions aux limites abordees au chapitre precedent (tem- 
peratures fixees ou flux fixe, libres ou rigides) reste identique apres lineari- 
sation. Les deux sortes de conditions aux limites pour la temperature sont 
done 

Temperatures fixees : 9 = en z = et z = 1. 

Flux fixe : jjf = enz = 0etz = l. (2.20) 

Les deux sortes de conditions aux limites en vitesses sont done 

Libres : || = g = 0etu; = en z = et z = 1. 

Rigides : u = v = w = en z = et z = 1. (2.21) 

5.2 Modes de Fourier 

Dans un premier temps, nous allons considerer le cas academique des condi- 
tions aux limites periodiques dans la direction z et nous notons 2d la periode 
pour la variable z dimensionnelle ce qui conduit a la periodicite 2 pour la vari- 
able z adimensionnee. Meme si ces conditions aux limites ne correspondent 
pas, a priori, a une situation physique, elles presentent l'avantage d'introduire 
simplement des notions qui seront utiles pour l'examen de conditions aux lim- 
ites physiques. 

La linearite des equations, la periodicite dans la direction verticale et l'invari- 
ance par translation dans les directions horizontales conduisent a la recherche 
de solutions sous la forme 

[U(x, t),9(x, t)] = [U m , 9 m ] e* k -- x - +s * (2.22) 

ou k = (ki,k2,ks) avec k\ et k2 quelconques et ^3 = n tt avec n entier 
pour assurer la periodicite en z. Les amplitudes U_ m et 9 m sont des quantites 



dJZ 

~dt 
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Figure 2.3: Conditions aux limites periodiques dans la direction z. 

complexes. Comme le probleme est lineaire, les solutions reelles sont obtenues 
en prenant la partie reelle des solutions complexes. 

Pour etudier la stabilite temporelle de l'etat conductif, on suppose que k 
est reel et s est complexe. Si Ton trouve des solutions non nulles telles que 
la partie reelle de s est positive, on pourra afnrmer que l'etat conductif est 
instable dans la mesure ou une petite perturbation est exponentiellement 
amplifiee, tout du moins dans le cadre de l'analyse lineaire. Dans le cas 
contraire, l'etat conductif est stable. 

Pour un vecteur d'onde donne k, on peut choisir le systeme de coordonnees de 
telle sorte que ki = 0. On voit alors que l'on se ramene au cas bidimensionnel 
oil £/(x, z) et 9(x, z) ne dependent pas de la coordonnee y et oil la composante 
v(x, z) ne joue aucun role. On peut done introduire une fonction de courant 
ip(x, z) et l'on considere alors des modes de Fourier de la forme 

&(x,z,t),e(x,z,t)\ = 0m] e !(M+f:32)+si (2.23) 

oil ip m et 6 m sont des amplitudes complexes. 



5.3 Equation aux valeurs propres 

En reportant l'expression des modes de Fourier dans le systeme d'equations 
2D on en obtient la transformee de Fourier qui s'ecrit 



-s k ip, 



m 

s 6 m 



iki RP9 m + Pk A ip r , 



iki ip n 



k 



(2.24) 



oil k = J k\ + k\ designe le module du vecteur d'onde k. On peut ecrire ce 
systeme lineaire sous la forme 

'^m\_(-Pk 2 -RPik 1 /k 2 \ (ip m 
iki ~k 2 



(2.25) 
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Ce systeme admet des solutions (ip m , 9 m ) non triviales lorsque le determinant 
du systeme est nul c'est-a-dire 



-Pk 2 -s -RPih/k 2 



i k\ 



-k 2 







ce que Ton peut ecrire sous la forme 

s 2 + s(P + 1) k 2 + Pk 4 - RPkl/k 2 = 



(2.26) 



(2.27) 



5.4 Determination du signe des valeurs propres 

Le discriminant de ce polynome de degre deux en s est (P— l) 2 k A +4 R Pk\jk 2 . 
Son signe etant toujours positif, les racines sont reelles. Pour determiner leurs 
signes, on ecrit alors 1' equation sous la forme 

s 2 + s( p + 1)fc 2 + p| {M~ R ) =0 ' (2 " 28) 

Pour R < k 6 /k 2 leur produit est positif et les deux racines sont negatives 
dans la mesure ou leur somme Test. Les deux modes de vecteur d'onde k 
correspondant aux deux racines s\ et S2 sont alors amortis exponentiellement. 
Pour R > k 6 /k 2 un de ces modes croit exponentiellement dans la mesure ou 
le produit des racines du polynome en s est negatif. 



5.5 Determination du Rayleigh critique 

On rappelle ici que la periodicite en z requiert qu'un vecteur d'onde k = 
(ki,ks) soit tel que k^ = n ir avec n entier. Etant donne un tel vecteur 
d'onde, la valeur critique du nombre de Rayleigh a partir de laquelle un des 
modes correspondant croit exponentiellement est done 

flc(fci,n) = M±^2!. (2.29) 

Lorsque n varie, ce nombre est minimum pour n = 1 et l'on peut alors tracer 
la courbe R c (ki) = R c (k\, 1). 

Le minimum de cette fonction est atteint pour k c = Tr/y/2 et vaut R c = 
Rc{k c ) = 27 7r 4 /4. On peut alors affirmer que pour R < R c , tous les modes de 
vecteurs d'ondes k = (k±, n it) decroissent exponentiellement. L'etat conductif 
est alors stable. Pour R > R c , une famine de modes incluant le vecteur d'onde 
(kc,ir) croissent exponentiellement. L'etat conductif est alors instable. 
Lorsque que l'on considere des conditions aux limites periodiques dans la 
direction horizontale, la composante k\ du vecteur d'onde est de la forme 
k\ = m ko ou ko est le nombre d'onde fondamental et m est entier. Ce 
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Figure 2.4: a) Courbes marginales R c (ki,n) pour n = 1,2,3. b) Domaine 
stabilite dans le plan (k\,R). 

nombre d'onde peut s'ecrire &o = 2 n (3 si f3 = d/L est le rapport de la 
longueur d sur la periode horizontale L en coordonnees dimensionnees (1//3 
est la periode horizontale en coordonnees adimensionnees). Le nombre de 
Rayleigh critique associe a ces conditions aux limites est alors R c (ko). 
On remarque que le seuil de l'instabilite de Rayleigh-Benard est independant 
du nombre de Prandtl P. 

6 Conditions aux limites realistes 

On s'interesse maintenant au calcul de stabilite du probleme de Rayleigh- 
Benard dans le cas de conditions au limites realistes. Le cas des conditions 
aux limites libres se deduit des calculs effectues pour les conditions aux limites 
periodiques dans la direction verticale. Le cas des conditions aux limites 
rigides necessite un peu plus d'efforts. 




Figure 2.5: Visualisation numerique de rouleaux de convection, a) Champ de 
vitesse. b) Champ de temperature. 
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6.1 Probleme aux valeurs propres 



Pour etudier le cas des conditions aux limites physiques (temperatures fixees 
ou flux fixe, libres ou rigides) on cherche des solutions sous la forme 



{x,z,t),6(x,z,t)] = fi>{z),0{z)}e iklX+st 



(2.30) 



La dependance avec les variables horizontales est sinuso'idale dans la mesure ou 
les equations sont invariantes par les translations dans ces directions. On peut 
alors choisir les coordonnees de telle sorte que la composante k<i du vecteur 
d'onde horizontal soit nulle, ce qui ramene le probleme au cas bidimensionnel 
et done a une fonction de courant ip. 
En reportant dans les equations on obtient le systeme 



s {D 2 - k\) * 

s e 



-ik 1 RP@ + P{D 2 - ki) 2 * 
-ifci* + {D 2 -k\) e 



ou D = ^ est l'operateur de derivation par rapport a la variable z. 
Les deux sortes de conditions aux limites en temperatures s'ecrivent 



Temperatures fixees : 
Flux fixe : 



= en z = et z = 1. 
DQ = en z = et z = 1. 



Les deux sortes de conditions aux limites en vitesses sont : 

Libres : * = et D 2 ^f = en z = et z 

Rigides : * = et D^f = en z = et z 



1. 
d. 



(2.31) 



(2.32) 



(2.33) 



6.2 Elimination de 

Pour eliminer la fonction @(z) on ecrit le systeme sous la forme 
s-P{D 2 -k 2 )] {D 2 - k 2 ) * = -ikxRP® 



-iki ^ 



(2.34) 



s-{D 2 -kl) 6 = 

On applique l'operateur s — (D 2 — k 2 ) a la premiere equation ce qui conduit a 

s-{D 2 -kl)\ \s-P{D 2 -k 2 ] {D 2 - fcf)* = -kjRPy . (2.35) 

Comme ^ = 0enz = 0etz = lala fois pour les conditions aux limites libres 
et rigides, les deux sortes de conditions aux limites en temperatures s'ecrivent 

Temperatures fixees : sD 2 ^ - P{D 2 - k\) 2 ^ = en z = 0, 1. 

Flux fixe : sD 3 ^ - P(D 2 - kj) 2 D^ = en z = 0, 1. (2.36) 
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7 Conditions aux limites en temperatures fixees 



On se place ici dans le cas des conditions aux limites en temperatures fixees et 
on examine successivement le cas des conditions aux limites libres et rigides. 

7.1 Conditions aux limites libres 

Pour le cas des conditions aux limites libres, on doit resoudre P equation 



s - (D 2 - k 2 )] [s-P (D 2 - kl)] {D 2 - A; 2 )* 
avec les conditions aux limites 



-k 2 RPV 



* = 0, D = et DH = 



en z = et z = 1. 



(2.37) 



(2.38) 



La symetrie z — > — z de ce systeme entraine que les solutions sont toutes de 
la forme = ifj m sin(n tt z) ou n est un nombre entier. 

On voit alors que la discussion sur la stabilite de l'etat conductif est la meme 
que celle qui a ete menee pour le cas des conditions aux limites periodiques, 
ce qui conduit au nombre de Rayleigh critique et au nombre d'onde horizontal 
critique suivant : 



Rr. 



27 TT 4 



657, 5 et k c 



IT 
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2,2 



(2.39) 



7.2 Conditions aux limites rigides 

Pour le cas des conditions aux limites rigides, on doit resoudre Pequation 



s-P(D 2 -kf) 



(D 2 - k\)^ = -k\RP^ 



(2.40) 



s-(D 2 -k 2 ) 

avec les conditions aux limites 

* = 0, = et s D 2 ^ - P(D 2 - kf) 2 ^ = enz = 0,l. (2.41) 

Les solutions ^f(z) sont la superposition de fonction exp (±A z) ou A 2 est l'une 
des trois solutions complexes de Pequation 



s-(A 2 -A; 2 ) s-P(\ 2 -k\) (X 2 -kl) + RPk 2 = 



(2.42) 



La recherche du Rayleigh critique s'effectue en posant s = dans les equations. 
En ecrivant que le systeme avec ses conditions aux limites admet des solutions 
non triviales, on se ramene a une equation permettant de determiner deux 
courbes R p (k\) et Ri(k\) pour les modes pairs et impairs. En calculant le 
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minimum de ces fonctions, on obtient le nombre de Rayleigh critique des con- 
ditions aux limites rigides ainsi le nombre d'onde horizontal correspondant : 

R c ~ 1707.762 et k c = 3.117 . (2.43) 

Le detail de ces calculs est explicite dans le livre "Hydrodynamics and Hy- 
dromagnetic stability" de S. Chandrasekhar (Oxford University Press 1961, 
Dover 1981). 



Conclusion 

Nous avons done enonce le modele permettant de decrire le mouvement con- 
vectif d'une couche fiuide chauffee par le bas dans le cadre de 1' approximation 
de Boussinesq. Differents jeux de conditions aux limites ont ete detaille. La 
solution solution stationnaire observee lorsque le forcage est faible correspond 
a un transfert de chaleur par conduction. Lorsque le nombre de Rayleigh 
depasse une valeur critique, il apparait des rouleaux de convection. Le calcul 
du nombre de Rayleigh critique a ete effectue en linearisant les equations du 
modeles autour de la solution conductive, puis en resolvant le probleme aux 
valeurs propres associe. 



FORMULAIRE 

EXERCICES ET PROBLEMES 



PROBLEME 2.1 Rayleigh-Benard avec rotation 



Rotation seule 

On considere un fiuide pesant et compris entre deux plaques planes hori- 
zontales delimitant les frontieres z = et z = d. On note g = —g e*- 3 -* 
le champ de gravite. Le fiuide et les plaques sont soumis a une rotation 
uniforme d'axe vertical Qo = S^ 3 \ On suppose que la temperature du 
fiuide est constante et que le mouvement de ce fiuide est regi par le modele 
des equations de Navier-Stokes incompressibles : 

div U = 

Pr (^ + 2Q AUj = -grad p - p r g e {3) + fi n AU . (2.44) 

La notation ^ = + U ■ grad designe la derivee particulaire. On note 
v = Pnl Pr et U(x, y, z,t) = u + v + w e^K On modelise l'interaction 
avec les plaques en z = et z = d par les conditions aux limites libres : 

du dv / . 

— = — = et w = en z = et z = d. (2.45) 

oz oz 
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EXERCICES ET PROBLEMES 



On suppose de plus que le debit moyen de l'ecoulement est nul. On peut done 
modeliser le probleme a l'aide de conditions aux limites periodiques sur U et 
p dans les directions horizontales. 

1) Montrer que ces hypotheses entrainent que le repos U = est une solution 
du probleme. Expliciter le champ de pression p associe en supposant que 
p = Pi sur la plaque inferieure. 

2) Montrer que Ton peut remplacer le terme —grad p — p r ge^ par le terme 
— grad II, ou H(x, t) est un champ que Ton explicitera. 

3) Indiquer les dimensions de la quantite G r = d 3 g/u 2 (dans le systeme SI 
par exemple). 

4) En fabriquant des unites de longueur et de temps a partir des parametres 
d et v, ecrire les equations du modelc sous forme adimensionnee. On 
definit le nombre de Taylor T a par la relation T a = 4 d 4 /v 2 . 

5) On s'interesse aux mouvements quasi-2D (v / 0) definis par 

U_ = [u(x, z, t), v(x, z, t), w(x, z, t)] et p = p(x,z,t). (2.46) 

Ecrire les equations adimensionnees pour les composantes u, v et w. 

6) Montrer qu'il existe une fonction de courant ip(x, z, t) telle que 

— i et — E- (2 ' 47) 

7) Montrer qu'une partie des conditions aux limites peuvent s 'ecrire ip = ipi 
en z = et if) = 4 ! 2 en z = 1 ou tpi et V2 sont des constantes. 

8) En s'appuyant sur le fait que le debit horizontal moyen doit etre nul, 
montrer que l'on peut choisir tpi = ip2 = 0. 

9) Conclure en explicitant toutes les conditions aux limites en fonction de 
ip et de v. 

10) Eliminer la pression en prenant le rotationnel de l'equation de quantite 
de mouvement en projetant sur e^ 2 \ En deduire le systeme d'equations 
aux derivees partielles regissant revolution de ip et v. On pourra utiliser 
la notation J(/ )S ) = g|-J g. 

11) Ecrire le systeme d'equations linearisees autour de l'etat de repos U = 0. 

12) On cherche des solutions de la forme 

mx,z,t),v(x,z,t)] = Mz),V(z)]e^ x+st . (2.48) 

En notant D = ^, ecrire le systeme d'equations differentielles ordinaires 
que verifient les profils *&(z) et V(z) ainsi les conditions aux limites as- 
sociees. 
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13) Montrer que ces profils sont de la forme = \l/ m F(nir z) et V(z) = 
V m G(mr z) ou fy m £ (E et 6 (T sont des constantes, n £ JV un entier 
strictement positif et F(C) et G(C) des fonctions trigonometriques simples 
que l'on precisera. Ecrire le systeme lineaire que doit verifier le couple 
(* m , V m ) en notant k le nombre d'onde defini par k 2 = k\ + n 2 7r 2 . 

14) En deduire que le couple (* m , V m ) est different de zero si et seulement 
si les nombres s, k\, n et T a sont relies par une relation polynomiale 
P(s, ki, k ;T a ) = 0. En deduire l'expression des valeurs propres s en 
fonction de k±, n et T a . 

15) Discuter la stabilite du repos en fonction de la valeur du nombre de Taylor 
T a . 

16) Montrer que cette analyse quasi-2D rend compte de la stabilite du repos 
vis-a-vis de perturbations 3D. 

17) Indiquer comment varient qualitativement les solutions du probleme lorsque 
la gravite g ou la densite p r varient. 



Convection et rotation 

On suppose maintenant que la temperature des plaques est maintenue aux 
temperatures constantes T = T\ en z = et T = T2 en z = d et que le 
mouvement de ce fluide est regi par le modele des equations de Navier-Stokes 
dans 1' approximation de Boussinesq : 

div U = 

-grad ( — + g z\+ag{T - T )e (3) + u AU 
k AT (2.49) 



dU 
~dt 



+ 2 Oo A U = 
dT 

~dt 



ou a est le coefficient d'expansion thermique de la loi d' etat p = p r [l — a(T — 
To)] et k est la coefficient de diffusivite thermique. Dans ces equations, on 
a neglige la puissance des forces interieures visqueuses 2 ^- D : D. Pour le 
champ de vitesse, on reste dans le cadre des conditions aux limites libres sur 
les plaques avec un debit moyen nul. 

18) En se restreignant aux instabilites de type stationnaires (valeur propre s 
devenant positive), decrire la surface R = R c {k\,T a ) de stabilite marginale 
(s = 0) de l'etat conductif dans l'espace des parametres (k\,R, T a ). Mon- 
trer, en particulier, que la rotation a toujours un effet stabilisant. 

19) Indiquer une methode permettant de calculer la surface de stabilite mar- 
ginale des instabilites oscillatoires (on ne demande par de faire les calculs). 



Corrige page 74 
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PROBLEME 2.2 



Explorations numeriques 



On donne ici quelques indication pour des sujets d'etudes numeriques et an- 
alytiques pouvant conduire a des travaux d'approfondissements. 

1) Simuler les equations de Ginzburg-Landau complexe avec un /x(x) vari- 
able et un terme c g A x rendant compte d'une vitesse de groupe. On 
pourra alors separer les domaines d'instabilite absolue ou convective. 

2) Simuler le modele suivant de couplage d'oscillateurs : 

du T ,/, x d 2 u 

M=-V'(u) + v— 2 . (2.50) 

3) Simuler numeriquement le choc de P equation de Burgers : 

du du d 2 u . „„. 

« + " (2 ' 51) 

4) Simuler la dispersion des ondes dans 1' equation 

du du „d 3 u . Kri . 



CORRIGES DES EXERCICES ET PROBLEMES 



Corrige 2.1 Rayleigh-Benard avec rotation 



Rotation seule 

l)L'absence de gradient de debit moyen entraine que le repos est solution. 
La pression associee est alors la pression hydrostatique p = pi — p r g z. 2) On 
peut poser LT = p — p\ + p r g z. C'est ici Pecart a la pression hydrostatique. 
3)La quantite G r est un nombre sans dimension. 4) On choisit [L] = d comme 
unite de longueur et [r] = d 2 /v comme unite de temps. L'unite de vitesse 
est alors v/d et l'unite de pression est p r v 2 jd 2 . En notant les grandeurs 
adimensionnees de la meme fagon que les grandeurs physiques, les equations 
s'ecrivent 

dwU = 

^+r|e (3) Af/ = -gradn + A[/. (2.53) 
5)Les equations du mouvement s'ecrivent 
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du du du dli i 

— + u— + w— = -— + T a 2 v + Au 

dt OX OZ ox 

dv dv dv i 

— + n — + = -TJu + Av 

dt dx dz 

dw dw dw dll 

+ u 7 - + w— = -— +Aw. 2.54 

dt dx dz dz 

6) La relation d'incompressibilite |^ + ^ = et la nature 2D du mouvement 
entrainent qu'il existe une fonction de courant. 7)Les conditions aux limites 
w = s'ecrivent ^ = 0enz = 0etz = l (unite adimensionnees). Etant 
donnee la geometrie (plane) des plaques, on en deduit que la fonction de 
courant ip est constante sur chacune des plaques. 8)L'integrale de u sur une 
verticale est une constante, comme le montre l'egalite Jq u dz = — Jq 1 ^ dz = 
ip\ — ip2- Si ipi n'etait pas egal a %p2 le debit moyen ne serait pas nul. Comme 
la fonction de courant est definie a une constant pres, on peut done ecrire 
V'l = ^2 = 0. 9)Les conditions aux limites sont done ip = 0, %^ = et v = 
en z = et z = 1. 10)Le rotationnel de la vitesse U_ est parallele a e^ 2 \ Les 
equations du mouvement sont alors 

d I dv 

-AV + J(V,AV) = -Ti—+A 2 i, 
dt dz 

^ + Jty,v) = Ti^ + Av. (2.55) 

11) La linearisation des equations consiste a negliger les termes non-lineaires 
J(ip, Aip) et J(ip,v). 12)Les profils ^f(z) et V(z) verifient le systeme 

s(D 2 -kl)^ = -rj Dv + (D 2 - k{) 2 m 

sv = T} D$> + {D 2 - k\)v . (2.56) 

avec les conditions aux limites ^(0) = #(1) = 0, D 2 ^(0) = £> 2 *(1) = et 
DV(0) = DV{1) = 0. 13)Lc choix F(Q = sin(C) et G(Q = cos(C) permet de 
satisfaire les conditions aux limites. En reportant dans les equations, on est 
conduit a resoudre le systeme : 



m 




(2.57) 



14) On obtient des solutions non triviales en imposant au determinant du 
systeme lineaire d'etre nul : 



-k 2 -s -mrTJ/k 2 
i 

nvrT a 2 -k 2 -s 



t 2 + k 2 ) 2 + n 2 7r 2 T a /k 2 = . (2.58) 
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Les valeurs propres sont done s±(ki, n; T a ) = —k 2 ± irnrT^ /k. 15)Le repos 
est toujours stable car les valeurs propres s ont une partie reelle negative. 
16)L'analyse 3D amene a considerer des perturbations de la forme 

[^(z),V(z)]e iklX+ik2y+st . 

II suffit alors de choisir une nouvelle base dont le premier vecteur unitaire est 
dans la direction du vecteur k\ e^ 1 ) + k2§_^ pour se ramener au cas quasi-2D 
deja etudie. 17)Le champ de vitesse est independant de la gravite g ou de la 
densite p r ou du nombre sans dimension G r . Seule l'expression de la pression 
hydrostatique, qui n'influence pas le mouvement, depend du poids du fluide. 



Convection et rotation 

18) Pour discuter la stabilite de la convection de Rayleigh-Benard soumise 
a une rotation, il suffit de combiner l'analyse du cours et des questions 
precedentes. II est preferable de choisir l'adimensionnalisation visqueuse basee 
sur l'unite de temps [t] = d 2 /v pour rester dans le cadre developpe lors des 
questions precedentes. On note alors 6 l'ecart de temperature au profil lineaire 
conductif en unites adimensionnees. Les nombres sans dimension controlant 
le probleme sont T a , le nombre de Taylor, R le nombre de Rayleigh et P, le 
nombre de Prandtl. En posant 

[ip, 6, v] = O(z), V(z)]e ikl x+st , (2.59) 

avec ty(z) = ^> m sin(n7rz), <d(z) = <d m sin(n-7rz) et V(z) = V n cos(n7rz), les 
valeurs propres s sont donnees par l'equation 



-k 2 -s 

—i ki 

i 

mrT£ 



iRh/k 2 
-k 2 -P s 





-nirTa 2 /k 2 


-k 2 -s 



. 



(2.60) 



La surface de stabilite marginale des instabilites stationnaires pour n fixe est 
obtenue en posant s = dans cette equation, ce qui conduit a l'equation 



-k 2 iRk^k 2 
-ik\ 
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ce qui s ecnt R c {ki,n,T a ) = ^ ^ • Oe nombre est 



minimum 



pour n = 1, ce qui conduit a R c [k\,T a 



(k 2 +n 2 ) 3 +7T 2 T a 



. Comme R c (k\, Ta) > 



R c (ki,0), on voit que la rotation retarde l'apparition de la convection. Le 
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nombre d'onde critique k c (T a ) et le nombre de Rayleigh critique R c {T a ) peu- 
vent etre obtenus en minimisant la fonction R c (k\, T a ) pour T a fixe (equation 
du troisieme degre en k\ a resoudre). 19)La courbe marginale des insta- 
bilites oscillatoires est obtenue en cherchant des valeurs propres s = ico imag- 
inaire pure. En reportant cette valeur dans le determiant et en annulant 
les parties reelle et imaginaire de l'equation de dispersion ainsi obtenue, on 
obtient l'equation de la surface de stabilite marginale sous la forme R = 
R° c sal (k\,T a ,P) ainsi que la valeur u(k\, T a , P) de la pulsation des oscilla- 
tions au seuil. 



QUESTIONNAIRES A CHOIX MULTIPLES 



CORRIGES DES QCM 
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OBJECTIFS PEDAGOGIQUES 

Cet article pedagogique a pour but de presenter, de la maniere la plus simple 
possible, les bifurcations "generiques", c'est-a-dire les plus courantes, decrivant 
la destabilisation d'un systeme dynamique lorsque Ton varie un de ses para- 
metres. Plusieurs objectifs de formation sont vises : 

• Assimiler la notion de stabilite d'un equilibre. 

• Comprendre la notion de parametre de controle. 
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NOTATIONS 



• Se familiariser avec certains elements de vocabulaire des systemes dy- 
namiques. 

• Assimiler les trois diagrammes de bifurcation des trois exemples etudies. 

• Comprendre la portee des trois exemples simples etudie en assimilant 
la notion de bifurcation generique. 

• Assimiler la notion de brisure de symetrie dans le cas de la bifurcation 
fourche. 

Les competences a acquerir lors de l'etude de cet article pedagogique sont les 
suivantes : 

• Etre capable de calculer les proprietes de stabilite de l'equilibre pour 
des systemes dynamiques reels ou complexes a un degre de liberte. 

• Etre capable de tracer les trois diagrammes de bifurcations correspon- 
dant aux trois bifurcations generiques de l'equilibre. 

PRE-REQUIS 

Le niveau requis pour la lecture de cette article pedagogique se situe autour 
de celui d'une Licence scientifique. 

NOTATIONS 



a Notation pour \J n/a ou yS—fi/a 

b = (b\, 62, •••) Vecteur de composantes bi 

cosh Cosinus hyperbolique 

x(t) Derivee de la fonction x(t) 

f'(x) Derivee de la fonction f(x) 

f(x) Fonction quelconque 

KifJ-'-, 20 Fonction vectorielle quelconque 

0(b) Du meme ordre de grandeur que b 

p Constante ou parametre de controle 

q Constante ou parametre de controle 

sinh Sinus hyperbolique 

S Symetrie point x — > —x 

<S Matrice 3x3 des composantes d'une symetrie 

t Temps 

tanh Tangente hyperbolique 

u(t) Ecart a l'equilibre 

x Variable d'un systeme dynamique 
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X 


Partie reelle de z 


Xe 


Equilibre d'un systeme dynamiques 


X + 


Equilibre particulier 


X- 


Equilibre particulier 


x(t) 


Trajectoire solution d'un systeme dynamique 


X 


Variables d'un systeme dynamique 




Equilibre d'un systeme dynamiques 




Equilibre particulier 


X 


Equilibre particulier 




Trajectoire solution d'un systeme dynamique 


V 


Partie imaginaire de z 


y(t) 


Notation pour —x(t) 


z = x + iy 


Variable complexe d'un systeme dynamique 


z(t) 


Trajectoire complexe solution d'un systeme dynamiques 


\z\ 


Module de z 


a 


Constante ou parametre de controle 




Constante ou parametre de controle 


7 


Constante ou parametre de controle 


5 


Constante ou parametre de controle 


e 


Argument de z 


A 


Valeur propre 


A* 


Parametre de controle de la bifurcation 


ft 


Parametres de controle de la bifurcation 


a 

£-c 


Parametres critiques 


P 


Module de z 




Constante ou parametre de controle 




Pulsation a P equilibre 



COURS ECRIT 
Introduction 

Un grand nombre de modeles issus de la mecanique des milieux continus 
s'ecrivent sous la forme de systemes d'equations aux derivees partielles para- 
boliques ou de systemes d'equations differentielles ordinaires, le temps etant la 
variable decrivant revolution du systeme. On parle de "systemes dynamiques" 
pour les equations differentielles ordinaires et on appelle "nombre de degres de 
liberte" le nombre d'equations couplees du systeme. La resolution numerique 
des equations aux derivees partielles se ramene toujours a un systeme dy- 
namique, en discretisant l'espace ou en projetant sur un ensemble fini de 
fonctions. Le nombre de degres de liberte est alors grand si la resolution 
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spatiale est grande. 

Une premiere approche pour l'etude de ces equations d'evolution consiste a 
rechercher les equilibres, c'est-a-dire les solutions stationnaires ne presentant 
pas d'evolution temporelle. On dit qu'un equilibre est stable si une petite 
perturbation de cet equilibre evolue dans le temps en restant petite et en 
convergeant vers zero. Le systeme revient alors dans sa position d'equilibre. 
Ces equilibres ont done de grandes chances d'etre observes sur le systeme 
physique etudie ou sur la simulation numerique du modele propose. 
L'etape suivante consiste a faire varier un ensemble de parametres controlant 
le systeme, comme par exemple l'intensite du forgage ou les dimensions geo- 
metriques du domaine. On regarde alors ce que deviennent les equilibres du 
modele, en particulier ceux qui etaient stables avant de modifier les parametres 
de controle. Lorsqu'en variant les valeurs des parametres de controle, un 
equilibre stable devient instable ou disparait, on dit que l'on est en presence 
d'une "bifurcation". L'etude systematique des bifurcations fait appel a un 
arsenal mathematique developpe (formes normales) qui a permis de classifier 
un grand nombre de comportements observes dans les systemes physiques. 
Ici nous ne nous interessons qu'aux bifurcations les plus courantes pouvant 
affecter un equilibre lorsque l'on varie des parametres de controle. Faire varier 
continuement n parametres de controle peut etre vu comme un chemin continu 
dans un espace vectoriel de dimension n. Nous verrons que la destabilisation 
d'un equilibre correspond au changement de signe de la partie reelle d'une 
des valeurs propres de l'operateur lineaire gouvernant les petites perturba- 
tions du voisinage de cet equilibre. On voit done que ce changement de signe 
definit une surface de codimension 1 dans l'espace de controle. En suivant un 
chemin continu dans l'espace de controle, on a de bonnes chances de traverser 
cette surface et d'assister alors a la destabilisation de l'equilibre. Une telle 
bifurcation de codimension 1 est dite "generique" . Pour observer des bifurca- 
tions non generiques, c'est-a-dire de codimension superieure a 1, il faut faire 
varier plusieurs parametres de controle a la fois. II est done impossible de les 
rencontrer en variant au hasard les parametres de controle. 
La theorie des formes normales appliquees aux systemes dynamiques montre 
qu'il n'y a que trois bifurcations generiques de l'equilibre. Lorsque la valeur 
propre responsable de la destabilisation est reelle, on parle de bifurcation 
stationnaire. On est alors en presence d'une bifurcation "noeud-col" dans le 
cas general, et "fourche" lorsque qu'il y a "brisure d'une symetrie". Lorsque 
qu'une paire de valeurs propres complexes conjuguees traversent l'axe des 
imaginaires, on est en presence d'une bifurcation "de Hopf". 
Nous allons etudier ces trois bifurcations (noeud-col, fourche et Hopf) sur les 
exemples les plus simples, qui sont, par ailleurs, les "formes normales" de tous 
les systemes dynamiques a plus grand nombre de degres de liberte dont les 
equilibres subissent l'une de ces bifurcations. 
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1 Bifurcation noeud-col 

On considere l'equation differentielle ordinaire suivante : 

x = fi + a x 2 . (3-1) 
On s'interesse a l'ensemble de toutes les solutions de ce systeme dynamique. 

1.1 Elements de vocabulaire sur les systemes dynamiques 

On appelle "espace des phases" l'ensemble des etats x S M. On appelle 
espace de controle l'ensemble des parametres (fi, a) G M 2 . Pour ce systeme 
dynamique, l'espace des phases est done de dimension 1. On dit qu'il s'agit 
d'un systeme dynamique a un degre de liberte. L'espace de controle est de 
dimension 2. 

Les systemes dynamiques autonomes (independants du temps) a un degre de 
liberte s'ecrivent sous la forme : 

x = f(x)=F( E ,x) (3.2) 

oil \x = (fix, fj,2, Mm) appartient a l'espace de controle. On suppose que 
f(x) est une fonction continuement derivable. Une solution quelconque x(t) 
d'un tel systeme dynamique decrit une trajectoire dans l'espace des phases. 
Elle est issue d'une condition initiale x(0) = xq a t = 0. Les trajectoires 
dans un espace des phases de dimension 1 ont pour support quatre types de 
courbes : un point, un segment ouvert, une demi-droite ouverte ou la droite 
tout entiere. Un point d'equilibre x e , qui forme a lui seul une trajectoire, est 
solution de l'equation f(x e ) = F(fi, x e ) = 0. 




Figure 3.1: Systeme dynamique a trois degres de liberies. 

Les systemes dynamiques autonomes (independants du temps) a n degres de 
liberte s'ecrivent sous la forme 

X = F(fi;X) (3.3) 
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avec 2L & Si n appartenant a l'espace des phases. II s'agit d'un systeme d'e- 
quations differentielles ordinaires couplees que Ton peut aussi ecrire sous la 
forme 

X i = F i {n 1 ,...,n m ;X 1 ,...,X n ) i = l,2,...,n. (3.4) 

On suppose que F est continuement differentiable. Une solution XL(t) de ce 
systeme decrit une trajectoire dans l'espace des phases. Un point d'equilibre 
X e , formant a lui seul une trajectoire, est solution de l'equation implicite 
F(p; X e ) = 0. On dit aussi que c'est une solution stationnaire ou encore un 
point critique du systeme dynamique. 

On montre (theoreme de Cauchy) qu'une condition initiale X(0) donne nais- 
sance a une et une seule trajectoire X(t). On en deduit que, pour un systeme 
dynamique autonome, deux trajectoires l'espace des phases ne se coupent ja- 
mais. Deux courbes contenant des trajectoires se coupent en une trajectoire 
qui est un point d'equilibre. 

Les systemes dynamiques non autonomes a n degres de liberte s'ecrivent sous 
la forme 

X = F(^,X,t) (3.5) 

avec 2L G M n appartenant a l'espace des phases. On peut les ramener au cas 
des systemes autonomes a n + 1 degres de liberte en considerant le temps t 
comme une variable X n+ i = t supplementaire de l'espace des phases et en 
ajoutant l'equation X n+ \ = 1. Notons que ce nouveau systeme dynamique ne 
peut pas admettre d'equilibre. 

1.2 Bifurcation noeud-col pour a < 0. 

On suppose que a est negatif. La recherche des solutions stationnaires conduit 
a resoudre l'equation f(x) = que Ton ecrit 

F(fx, a; x) = fj, + a x 2 = . (3-6) 

L'ensemble des solutions stationnaire est Pensemble vide pour /j, < 0, le 
singleton {0} pour \i = et {—a, a} pour fi > avec a = \f— fi/a. 
On dit qu'un point d'equilibre x e , verifiant done f(x e ) = 0, est stable, si 
une petite perturbation n(0) autour de cet equilibre est amortie au cours du 
temps. Autrement dit, si la condition initiale x(0) = x e + u(0), avec u(0) 
"petit", donne naissance a une trajectoire x(t) = x e + u(t) convergeant vers 
x e , e'est-a-dire tel que u(t) tend vers zero lorsque t tend vers l'infini. 
Pour determiner si un equilibre x e d'un systeme dynamique x = f(x) est 
stable, on effectue le developpement limite suivant 



x(t) = u(t) = f[x e + u(t)} = f'(x e ) u(t) + 0(\u(t)\ 2 ) . (3.7) 
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On a utilise ici x e = et f(x e ) = 0. Tant que u(t) reste petit, la solution 
peut etre approchee en resolvant l'equation 

u = A u (3.8) 

avec A = f'(x e ). On dit que l'on a linearise le systeme dynamique autour de 
l'equilibre x e . 




Figure 3.2: Diagramme (fJ-,x) de la bifurcation noeud-col. a) Conditions 
initiales [fi,x(0)] repartie regulierement. b) Positions des points [/i, x(t)} au 
bout du temps t. 

En tenant compte de sa condition initiale, la solution de ce systeme s'ecrit 

u(t) = u(0) exp(A t) . (3.9) 

On voit alors que, pour A = f'(x e ) < 0, la perturbation u(t) tend vers zero. 
L'equilibre est alors stable. Lorsque f'(x e ) > 0, l'equilibre est instable et le 
systeme linearise cesse d'etre valide lorsque u(t), qui croit exponentiellement, 
devient trop grand. Pour le cas critique f'(x e ) = 0, on dit que l'equilibre est 
marginal. Pour cette valeur tres particuliere, l'etude de stabilite, lorsqu'on 
s'y interesse, necessite de pousser le developpement limite de f(x) autour de 
x e jusqu'au terme d'ordre 2, ce qui conduit a une equation non lineaire. 
Dans le cas general d'un systeme dynamique a n degres de libertes X = F(X), 
on montre que la stabilite d'un equilibre X_ e verifiant F(X e ) = depend du 
signe de la partie reelle de la matrice DF{X_ e ) jacobienne de l'application 
F et de composantes dFi/dXj. En effet, si on pose X = X_ e + U, on doit 
resoudre U = DF(X e ) U + 0(||{7|| 2 ). Dans le cas n = 1 qui nous interesse 
ici, la jacobienne n'a qu'une composante qui est /'. 

Pour l'exemple particulier qui nous interesse ici, on a f'(x) = 2 a x. La 
stabilite des equilibres est done donnee par les nombres A = f'(x_) = 2^/— a [i 
pour x- = —a et A = f'(x + ) = —2^/— a fi pour x+ = a. On peut done 
affirmer que l'equilibre x_ est instable et que l'equilibre x+ est stable. 
Un methode graphique pour la determination de la stabilite des equilibres d'un 
systeme dynamique a un seul degre de liberte consiste a tracer la fonction f(x) . 
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Ses zeros correspondent aux equilibres et les pentes en ces points indiquent 
leur stabilite : stable si la pente est negative, instable si la pente est positive. 

1.3 Cas a > et diagrammes de bifurcation 

L'examen du cas a > est identique. L'ensemble des solutions est {—a, a} 
pour fi < avec a = ^—fx/a, le singleton {0} pour fx = et l'ensemble vide 

pour fx > 0. On calcule A = f'(x-) = —2yJ—a fx pour l'equilibre x_ = —a 
et A = f'(x + ) = 2yJ—a [x pour x + = a. On peut done affirmer que l'equilibre 
X- est stable et que l'equilibre x+ est instable. 

On peut rassembler tous ces resultats en tracant deux diagrammes de bifur- 
cation. Pour a < 0, on trace d'abord les courbes des solutions x e (fx) dans un 
plan (/x,x) (parabole orientee vers les fx positifs). Par convention, on trace en 
trait plein la courbe des equilibres stables, ici x+(/x), et en traits pointilles la 
courbe des equilibre instables, ici x-(/x). Pour un fx fixe, il est alors facile de 
tracer les trajectoires. Celles-ci s'ecartent des points instables pour converger 
vers les points stables. On trace done des droites dirigees vers les x > pour 
ji < 0. Pour /i>0on trace des demi-droites allant de — oo a x_, des segments 
de droites allant de x_ a x+, puis des demi-droites allant de x+ a +oo. 
Pour a > 0, on trace un nouveau diagramme de bifurcation ayant Failure 
suivante. Pour [i < 0, on trace des demi-droites allant de — oo a x_, des 
segments de droites allant de x_ a x+, puis des demi-droites allant de x+ a 
oo. Pour /i > on trace des droites allant de — oo a +oo. 
On remarque que Ton passe du cas a < au cas a > en changeant a — > —a, 
(i — > —fj,, t — > — t et x — > x ou bien encore en changeant a — > —a, \i — > — /i, 

1 — > i et x — > —x. Ces symetries sont visibles en comparant les diagrammes de 
bifurcations et le sens des trajectoires, ou encore le trace des diverses fonctions 
x(t) en fonction du temps. 

Cette etude montre que la valeur \i = correspond a un changement du 
nombre d'equilibres. Pour cette valeur du parametre de controle, un equilibre 
stable et un equilibre instable se rejoignent. On parle alors de bifurcation 
"noeud-col" (saddle-node en anglais). 

2 Bifurcation fourche 

On considere maintenant l'equation differentielle ordinaire suivante : 

x = / ux + ax 3 . (3.10) 

On cherche a tracer le diagramme de bifurcation de ce systeme dynamique. 
On note des a present la symetrie fx — > — fx, t — > — t, x — ► — x. Une autre 
symetrie est obtenue en ne modifiant pas les parametres de controle ni le sens 
de l'axe des temps. Cette symetrie est la symetrie point x — > — x dans M 
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et Ton voit que si x{t) est une solution, la fonction y(t) = —x(t) est aussi 
une solution. Cette symetrie particuliere joue un role tres important pour la 
bifurcation fourche, comme nous allons le voir. 

2.1 Diagramme de bifurcation 

On considere tout d'abord le cas a < 0. L'ensemble des equilibres est obtenu 
en cherchant les racines de l'equation /(x e ) = fi x e + a Xg = 0. Cet ensemble 
est {0} pour fi < et {— a, 0,a} avec a = \fjjja pour \i > 0. La stabilite de 
depend du signe de /'(0) = fi. Pour \i < l'equilibre est stable, il est 
instable pour \i > 0. 

La stabilite de x_ = —a est la meme que celle de x+ = a dans la mesure 
ou f'(x-) = f'(x + ) = —2fi. Pour fjL > 0, c'est-a-dire lorsqu'ils existent, 
ces equilibres sont stables. On trace alors le diagramme de bifurcation cor- 
respondant a ce cas a < 0. Pour \i < 0, deux trajectoires (demi-droites) 
venant respectivement de — oo et +oo convergent le point d'equilibre x = 0. 
Pour fx > deux trajectoires (segments de droites) s'ecartent de x = pour 
converger respectivement vers x_ = —a et x+ = a, et deux trajectoires (demi- 
droites) s'ecartent respectivement de ces equilibres pour converger vers — oo 
et +oo. Le nombre d'equilibre et leur stabilite peuvent se voir graphiquement 
en tragant la fonction /(x) = F(fi, a; x) lorsque \x est negatif ou positif. 

4 O > o 1 o i o I o « o 6 J >F ! 1 ! | ! ! ! - 




-i - ? » » » » 9 ° » ° ° ♦ 



Figure 3.3: Diagramme (/i, x) de la bifurcation fourche. a) Conditions initiales 
[fj,, x(0)] repartie regulierement. b) Positions des points [fi, x(t)\ au bout du 
temps t. 

Dans le cas a > 0, l'ensemble des equilibres est {— a,0,a} avec a = y/—/j,/a 
pour fi < et {0} pour /i > 0. La stabilite de l'equilibre x = est la meme que 
pour le cas a < : stable pour \x < et instable pour \i > 0. La stabilite de 
x_ = —a et x+ = a depend du signe de f'(x-) = f'(x + ) = —2a// : ces deux 
equilibres sont instables sur leur domaine d'existence fi < 0. Pour fi < 0, le di- 
agramme de bifurcation fait apparaitre deux trajectoires (demi-droites) allant 
respectivement de —a a — oo et de +a a +oo ainsi qu'une trajectoire (segment 
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de droite) allant de —a a +a. Pour /i > 0, deux trajectories s'ecartent de 
l'equilibre x = pour aller respectivement vers — oo ou +oo. 

2.2 Brisure de symetrie 

On dit que la bifurcation est supercritique dans le cas a < et sous-critique 
(ou subcritique) dans le cas a > 0. Dans le cas supercritique, l'etat d'equilibre 
x = se destabilise lorsque \x devient positif pour donner naissance a deux 
etats d'equilibre symetriques dans l'amplitude a croit progressivement avec 
jjL comme ^/Jl. Dans le cas sous-critique, une perturbation d'amplitude finie 
avant le seuil de destabilisation de l'equilibre x = peut conduire a une 
trajectoire divergente si cette perturbation est sufnsamment grande. 
La symetrie point S : x — > — x joue un role important dans cette bifurcation. 
On rappelle que si x(t) est une solution la fonction y(t) = S x(t) = —x(t) 
est aussi solution. En ecrivant le systeme dynamique sous la forme x = 
f(x), l'invariance de ce systeme par la symetrie S s'ecrit S x(t) = f[Sx(t)] = 
S f[x(t)]. On remarque que l'equilibre x = est invariant par la symetrie 
(50 = 0) et que les equilibres bifurques sont deux a deux symetriques (S X- = 
x + et S x + = x ). 

Dans le cas general d'un systeme dynamique X = F_(X), on dit que l'operateur 
lineaire S_ est une symetrie du systeme si F(S X) = S F(X) pour tout vecteur 
d'etat 2L G Si n de l'espace des phases. Par definition S_ est un matrice n x n 
verifiant S S_ = I (identite) . Cette definition entraine que si X(i) est solution, 
alors Y_(t) = S_ X_(t) est aussi solution. Un equilibre X e est invariant par la 
symetrie si ^ X_ e = X_ e . En changeant l'origine de l'espace des phase, on peut 
toujours se ramener au cas X e = 0, ce que nous supposons ici. 
Lorsque l'equilibre symetrique X_ e = se destabilise, il peut donner naissance 
a de nouveaux equilibres qui sont eux-memes symetriques. Dans ce cas, il 
n'y a pas brisure de symetrie. Lorsqu'il y a brisure de symetrie, les equilibres 
bifurquees formes des paires (X_ , X + ) telles que S X__ = X_ + et S_ X_ + = X__ . 
Dans les cas generiques, une seule paire de telles solutions prend naissance a 
la bifurcation, et par symetrie, X_ e = reste un equilibre, devenu instable. 
On montre que la brisure de symetrie depend du vecteur propre l'operateur 
lineaire decrivant revolution des petites perturbations autour de X_ e = 0, 
associe a la valeur propre qui s'annule a la bifurcation. Lorsque ce vecteur 
propre est invariant par ^, la symetrie n'est pas brisee. Elle l'est dans le cas 
contraire. 

3 Bifurcation de Hopf 

On considere maintenant P equation differentielle ordinaire complexe suivante : 
i = (fi + iu) z + (a + i(3) \z\ 2 z (3.11) 
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oil z £ <T est une variable complexe. En posant z = x + i y, on se ramene a 
un systeme dynamique a deux degres de libertes. L'espace des phases est de 
dimension 2, tandis que l'espace de controle decrivant les parametres p, uj, a 
et (3 est de dimension 4. 

Tres souvent, on s'interesse a la partie reelle des solutions complexes z(t) de 
cette equation. Ces amplitudes complexes sont par exemple la projection d'un 
etat sur un vecteur propre associe a une valeur propre complexe (p + iu). On 
cherche a decrire le diagramme de bifurcation de ce systeme. 

3.1 Destabilisation de l'equilibre 

En posant z(t) = p(t) exp[i 9{t)\ avec p(t) > 0, le systeme s'ecrit 

p = fJiP + ap 3 et 6 = uj + (3p 2 . (3.12) 

On remarque que revolution du module p(t) est decouplee de celle de 9{t) 
et que son etude a deja ete realisee lors de l'etude de la bifurcation fourche. 
L 'evolution de la phase 6(t) depend alors de celle de p(t) de maniere simple. 
On peut alors tracer les diagrammes de bifurcation de la bifurcation de Hopf 
a la lumiere des diagrammes de bifurcation de la bifurcation fourche. Ces 
nouveaux diagrammes peuvent etre representee dans l'espace tridimensionnel 
(/x, x, y) oil x et y sont les parties reelles et imaginaires de z. Pour p fixe, les 
trajectoires parcourent done des courbes dans un espace de dimension 2 que 
Ton peut voir comme M 2 ou W. 

Pour p < 0, l'equilibre stable p = de la bifurcation fourche, correspond 
a un equilibre stable z = 0. La decomposition polaire z = pexp(i6) n'est 
pas valable pour cet equilibre meme (la phase est indeterminee) mais reste 
pertinente dans tout voisinage. Cet equilibre devient done instable pour p > 
0. On peut retrouver directement ce resultant en remarquant que le systeme 
s'ecrit 

z = (p + iu) z + 0{\z\ 3 ) (3.13) 

lorsque z est petit, e'est-a-dire proche de l'equilibre z = 0. L'integration du 
systeme linearise conduit alors a z(t) = z(0) exp(pt) exp(i uj t). On voit que 
\z(t)\ = exp(pt) tend vers zero pour p < et ne reste pas petit p > 0. Dans le 
voisinage de zero, les trajectoires forment alors des spirales logarithmiques de 
raison p et de pulsation to. Ces spirales convergent vers dans le cas stable, 
ou divergent a partir de dans le cas instable. 

3.2 Diagramme de bifurcation 

Considerons tout d'abord le cas supercritique a < 0. Pour p < 0, les modules 
p(t) des trajectoires vont de — oo vers 0. La vitesse angulaire 9 = uj + (5p 2 est 
tres rapide a l'infini et se stabilise a uj lorsque les trajectoires arrivent dans 
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le voisinage de zero. Si cette vitesse angulaire garde le meme signe de l'infini 
a lo, les trajectoires dans W ont Failure de spirales tournant toutes dans le 
meme sens. 

Pour p > 0, L'equilibre stable p e = a avec a = ^j—pja correspond a une 
trajectoire periodique stable z(t) = aexp(iuj e t) avec m e = uj + (3a 2 . Cette 
trajectoire decrit un cercle dans la mesure ou le module est constant. Si Ton 
perturbe cette solution en considerant une condition initiale dans le voisinage 
proche de la trajectoire periodique, la petite perturbation ainsi introduite 
convergera vers zero. On dit que cette trajectoire periode est un "cycle limite" , 
ce mot etant synonyme de stabilite. 




Figure 3.4: Diagramme (p,x) de la bifurcation de Hopf. a) Conditions ini- 
tiales [p, x(0)] repartie regulierement. b) Positions des points [p,x(t)] au bout 
du temps t. 

Le module p des autres trajectoires va de a a = yj—p/a ou de a a l'infini. 
Dans le premier cas, ces trajectoires decrivent des spirales divergente autour 
de zero et qui approchent exponentiellement le cycle limite par l'interieur. 
Dans le second cas, ces trajectoires decrivent des spirales venant de l'infini et 
convergeant vers le cycle limite par son exterieur. 

Le cas sous-critique a > conduit au diagramme de bifurcation suivant. 
Pour p < 0, il existe un trajectoire periode z(t) = aexp(iui e t), mais celle-ci 
est instable. Une perturbation vers l'interieur du cercle donne naissance a 
une spirale qui convergent vers zero. Une perturbation a l'exterieur du cercle 
conduit a une spirale tendant vers l'infini. Pour p > 0, les trajectoires sont 
des spirales s'ecartant de et convergeant vers l'infini. 

4 Formes normales et catastrophes 

Les systemes dynamiques x = p + ax 2 et x = px + ax 3 avec x £ M ou 
z = (p + iuj) z + (a + if3) \z\ 2 z avec z G (U, decrivent les trois bifurcations 
generiques de l'equilibre au sens qui a ete defini dans l'introduction. Les 
bifurcations noeud-col et fourche correspondent a des instabilites stationnaires 
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tandis que la bifurcation de Hopf correspond a une instability oscillatoire. 
Plusieurs remarques peuvent etre effectuees pour completer Petude de ces 
bifurcations. 



4.1 Bifurcations imparfaites 

Notons tout d'abord que la notion de genericite est relative a la classe de 
systemes dynamiques accessibles en parcourant l'espace de controle. Pour la 
bifurcation fourche, on a suppose que le systeme etait invariant par symetrie 
S pour toutes les valeurs de l'espace de controle. Cette situation se rencontre 
souvent dans la nature oil les symetries sont nombreuses. C'est justement 
la brisure de ces symetries (un equilibre bifurque n'est plus symetrique mais 
le systeme dynamique lui-meme reste symetrique) qui donne lieu a des bi- 
furcation fourches. Le mathematiciens pourrait cependant pretendre que la 
bifurcation fourche n'est pas generique dans un espace de controle ou Ton 
s'autorise a ce que le systeme dynamique lui-meme ne soit plus symetrique. 
C'est le cas par exemple du systeme dynamique x = \ix + ax 3 + e. Dans 
ce cas symetrie x — > — x n'est plus une symetrie du systeme des que e est 
different de zero. On peut alors tracer des diagrammes de bifurcations dans 
le plan (/x, x) pour e petit, et observer que la bifurcation fourche degenere en 
une bifurcation noeud-col, proche d'une branche stable (cas super-critique) 
ou instable (cas sous-critique). On parle alors de bifurcation imparfaite. 



4.2 Exposants critiques 

Pour les trois bifurcations, l'amplitude des etats bifurques croit comme la 
racine carre de Pecart au seuil. La theorie des formes normale, qui ramene les 
bifurcations generiques de P equilibre a Pun de ces trois modeles, montre que 
cet exposant critique 1/2 est universel, c'est-a-dire independant du systeme 
dynamique etudie. Ce resultat est important lorsque Pon veut determiner 
experiment ale me nt un seuil de bifurcation a partir de la mesure des ampli- 
tudes de Petat bifurque (cas supercritique) . 

Un autre exposant critique universel concerne les temps caracteristiques de 
convergence vers un equilibre dans le cas stable ou de croissance exponen- 
tielle a partir de l'equilibre dans le cas instable. Ces temps varient comme 
Pinverse de Pecart au seuil, faisant ainsi apparaitre l'exposant critique — 1. 
La determination de ce temps en fonction d'un parametre de controle per- 
met de mesure experiment alement un seuil de bifurcation, meme dans le cas 
sous-critique. 
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4.3 Bifurcations sous-critiques et hysteresis 

La theorie des formes normales indiquent que le comportcment des trois 
systemes dynamiques considered pour decrire les bifurcations noeud-col, four- 
che et de Hopf, est universel dans le voisinage des equilibres obtenus pour 
des systemes dynamiques quelconques et pour des valeurs des parametres 
proches du seuil. Les trajectoires qui tendent vers l'infini pour les formes nor- 
males peuvent converger vers d'autres equilibres ou regimes pour le systeme 
dynamiques que l'on souhaite etudier. Cette situation se rencontre en parti- 
culier pour les bifurcations fourche ou de Hopf sous-critiques, pour lesquelles 
il n'existe pas de nouvel equilibre une fois l'equilibre initial destabilise. 
On peut donner des exemples de telles situations a l'aide de systemes dy- 
namiques simples dont l'universalite n'est pas comparable a celle des trois 
formes normales que Ton a etudie. Considerons par exemple le systeme dy- 
namique reel 

x = fix + ax 3 + 7X 5 (3.14) 
avec x £ M et le systeme dynamique complexe 

z = (n + iuj) z + (a + i/3) \z\ 2 z + (7 + iS) \z\ 4 z (3.15) 

avec z G W. Dans le cas sous-critique a > on observe un phenomene 
d'hysteresis. Lorsque fx devient positif, une petite perturbation de l'equilibre 
devenu instable s'amplifie pour converger vers un nouveau regime (equilibre ou 
cycle limite) d'amplitude finie des le seuil (a condition que 7 soit negatif). Ce 
nouveau regime stable persiste lorsque Ton diminue le parametre de controle 
jjL, meme en-dega du seuil de la bifurcation. En continuant a diminuer fi, 
regime subit une bifurcation saddle-node (d'equilibres ou de cycles) et le 
systeme retombe vers l'equilibre stable 0. On a ainsi decrit un cycle d'hysteresis 
et il existe une gamme de valeurs du parametre fi pour laquelle deux regimes 
stables coexistent. Dans ce cas la position de la condition initiale dans l'un 
ou l'autre des bassins d'attraction des deux regimes conditionne l'etat asymp- 
totique du systeme. 

4.4 Theoreme des fonctions implicites et catastrophes 

La recherche des equilibres d'un systeme dynamique x = f(x) = F(fi; x) con- 
duit la resolution de Pequation implicite F(fj,; x e ) = 0. Le theoreme des fonc- 
tions implicite indique que l'on peut prolonger une courbe solution x e (fi) tant 
que ^(/i;x e ) = f'(x e ) 7^ 0. Les valeurs critiques n c telles que ^(/j,;x e ) = 
f'(x e ) = correspondent a des bifurcations ou le nombre d'equilibre change. 
L'analyse de stabilite de l'equilibre x e du systeme dynamique x = f(x) faisant 
intervenir le signe de f'(x e ), il est done normal qu'un changement de stabilite 
soit accompagne d'un changement dans le nombre d'equilibres. C'est le cas 
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de la bifurcation noeud-col et de la bifurcation fourche. Dans le cas de la 
bifurcation de Hopf, le theoreme des fonctions implicites n'est pas viole sur la 
recherche des equilibres z = x + iy. II l'est pour la recherche equilibres pour 
le module p en posant z = pexp(i6). 

L 'etude systematique des cas de violation du theoreme des fonctions implicites 
et de leur voisinage dans l'espace de controle fait l'objet de la theorie des 
catastrophes. La catastrophe la plus simple est la pli, represents par sa forme 
normale F(fi; x e ) = fi + ax 2 . En variant un seul parametre, on passe de zero 
a deux solutions. C'est ce qui se passe pour les equilibres de la bifurcation 
noeud-col. Par rapport a l'etude des catastrophes, qui ne s'interesse qu'aux 
equilibres, l'analyse dynamique apporte des informations sur la stabilite de 
ces equilibres. 

4.5 Catastrophes de codimension 2 : la fronce 

II est interessant ici d'etudier, a titre d'exemple, un systeme dynamique faisant 
apparaitre la catastrophe de dimension 2 qu'est la fronce. Une fronce est 
l'intersection de deux plis dans l'espace de controle. Nous considerons done 
le systeme dynamique 

x = F(p, q; x) = f(x) = — q + px — x 3 (3.16) 

dont l'espace de controle est l'ensemble des couples (p, q) £ M 2 . Considerons 
tout d'abord des chemins dans l'espace de controle ou p est maintenu constant 
et ou q decrit l'espace reel. Si p < 0, il n'existe qu'un seul equilibre comme le 
montre le trace de la courbe f(x) = —q+px — x 3 . Cet equilibre est alors stable 
car la derivee f'(x) = p — 3x 2 est negative. Si p > 0, le trace de la courbe 
f(x) fait apparaitre deux extrema egaux a — q ± 2(p/3) 3//2 . Si \q\ < 2(p/3) 3//2 , 
e'est-a-dire q 2 < 4p 3 /27, il existe trois equilibres, un qui est instable entoure 
de deux qui sont stables. Dans ce cas, le diagramme de bifurcation {p fixe, q 
variable) fait apparaitre deux bifurcations noeud-col pour les valeurs critiques 
q = 2(p/3) 3 / 2 et q = — 2(p/3) 3 / 2 . On note aussi un phenomene d'hysteresis : 
en suivant un equilibre lorsque q varie, le systeme bascule brutalement vers 
l'autre equilibre lorsque le seuil de la bifurcation noeud-col est depasse. En 
faisant alors varier le parametres dans l'autre sens, on retourne a l'equilibre 
stable initial mais en depassant l'autre seuil de bifurcation. Ainsi, pour une 
valeur de q situee entre les deux seuils, on peut atteindre deux equilibres 
stables selon la condition initiale ou le chemin suivi dans l'espace de controle. 
On resume cette analyse en tracant la courbe 27 q 2 — 4p 3 dans l'espace de 
controle (p, q) . Cette courbe fait apparaitre un point de rebroussement au 
voisinage de (p, q) = 0. On peut aussi tracer la courbe F(p, q, x e ) = dans 
l'espace (p,q,x e ). Cette visualisation, que Ton peut essayer de materialiser 
en deformant convenablement une feuille de papier, donne le nom de fronce a 
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cette catastrophe de codimension 2. Notant pas ailleurs que des chemins q fixe 
et p variable produisent une bifurcation fourche dans le cas particulier q = 
et des bifurcations fourches imparfaites dans le voisinage de cette valeur. 
La theorie des catastrophes combinee a la theorie des formes normales montre 
que cette description qualitative, ainsi que l'exposant critique 3/2 du point 
de rebroussement dans l'espace de controle, sont des caracteristiques uni- 
verselles des systemes pour lesquels deux bifurcations noeud-cols se coupent 
dans l'espace de controle. La portee de ce resultat explique certainement le 
succes de la theorie des catastrophes dans 1' etude des systemes physiques. 



Conclusion 



Nous avons etudie en detail trois equations modeles des bifurcations noeud- 
col (x = [i + a x 2 ), fourche (x = fix + a x 2 ) et de Hopf (i = (/x + iuS)z + (a + 
i0)\z\ 2 z). 

Les bifurcations noeud-col et de Hopf sont generiques dans l'ensemble des 
systemes dynamiques quelconques (sans symetrie) de M n des lors que Ton 
s'interesse a la stabilite d'un equilibre en faisant varier un parametre. Ceci 
signifie que, sur une surface de codimension un dans l'espace des parametres 
de controles, un equilibre stable se destabilisera en suivant le diagramme de 
bifurcation de l'une des ces deux bifurcations. Pour determiner laquelle, il faut 
lineariser le systeme dynamique autour de l'equilibre et analyser les valeurs 
propres de la matrice jacobienne ainsi obtenu. En faisant varier un parametre 
de controle, l'equilibre se destabilise par passage d'une valeur propre a travers 
l'axe des imaginaires purs. S'il s'agit d'une valeur propres reelles, on est en 
presence d'une bifurcation noeud-col. S'il s'agit d'une paire de valeurs propres 
complexes conjugees, on est en presence d'une bifurcation de Hopf. 
La bifurcation fourche est generique dans l'ensemble des systemes dynamiques 
invariant par une symetrie 5 de M n lorsque Ton suit un equilibre lui-meme 
invariant par cette symetrie. Cette bifurcation correspond alors a une brisure 
de symetrie et correspond au changement de signe d'une valeur propre reelle. 
Les trois equations modeles decrivant les trois bifurcations generique de l'equilibre 
sont en fait l'ordre dominant d'un developpement asymptotique au voisinage 
des seuils de changement de stabilite. On dit que ces equations sont les 
"formes normales" des ces trois bifurcations. On montre alors que la topolo- 
gie des trajectoires du systeme complet au voisinage des equilibres est la meme 
que celle des formes normales. 
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SYSTEMES DYNAMIQUES 

N degres de liberte 

X = F{p,X) 

avec XeM n et[i£ M m . 
Un degre de liberte 

x = f(x) = F(wx) 

Stabite 

Equilibre x e avec f(x e ) stable ssi \i = f'(x e ) < 0. 



BIFURCATIONS GENERIQUES 
Bifurcation noeud-col 



Bifurcation fourche 



x = /j, + a x 2 



x = fi x + a x 3 



Bifurcation de Hopf 

z = (n + ioj) z + (a + i j3) \z\ 2 z 



EXERCICES ET PROBLEMES 



EXERCICE 3.1 Bifurcation noeud-col analytique 



On peut donner l'expression analytique des solutions x{t) du systeme dy- 
namique i = /i + ax 2 decrivant la bifurcation noeud-col. Le trace prealable 
des diagrammes de bifurcation permet d'eclairer la discussion en fonction des 
valeurs des parametres lors du calcul analytique. 

1) Donner l'expression analytique des solutions dans le cas oi\x > 0. 

2) Donner l'expression analytique des solutions dans le cas afi < 0. 

3) Verifier que les diagrammes de bifurcation obtenus pour afi < ou a/j, > 
correspondent bien au comportement des solutions analytiques. 

Corrige page 100 



EXERCICE 3.2 Bifurcation fourche analytique 



Comme pour le cas de la bifurcation noeud-col, il est interessant de considerer 
les solutions analytiques du systeme dynamique x = \x x + ax 3 , decrivant la 
bifurcation fourche, a la lumiere des diagrammes de bifurcations qui donnent 
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le comportement asymptotique a priori des solutions. 

1) Donner l'expression analytique des solutions dans le cas a/j, > 0. 

2) Donner l'expression analytique des solutions dans le cas a/j, < 0. 

3) Verifier que les diagrammes de bifurcation obtenus pour afi < ou afi > 
correspondent bien au comportement des solutions analytiques. 

Corrige page 102 



EXERCICE 3.3 | Travelling or Standing waves 

On s'interesse a des solutions d'un modele ID qui s'ecrivent sous la forme : 

u(x,t) = Re[ Zl (t) e iklX + z 2 (t) e~ iklX ] (3.17) 

ou Re designe la partie reelle. On suppose que les amplitudes complexes z\ 
et z 2 verifient le systeme : 

z\ = ([i + iu>) z\ - (a + i j3)\zi\ 2 zi - + i5) \z 2 \ 2 zi 

z 2 = (ii + i uj) z 2 - (a + i (3) \z 2 \ 2 z 2 - (^ + i 5) \zi\ 2 z 2 . (3.18) 

On suppose que fi > 0, a > et [3 > 0. 

1) Montrer que l'etude de ce systeme peut se ramener, sans perte de generalite, 
a celle du systeme reel : 

R 1 = R 1 - (r\ + v R,£) Ri et R 2 = R 2 - (r 2 2 + v R^j R 2 , 

(3.19) 

avec i?i > 0, R 2 > et v > 0. 

2) On s'interesse aux solutions stationnaires telles que R 2 = 0. Calculer la 
valeur d'equilibre R\ = Rtw 

3) Dans le cas ou (R\,R 2 ) = (Rt w ,0), decrire revolution spatio-temporelle 
de la fonction u(x, t) correspondante. Illustrer cette description a l'aide 
de quelques traces schematiques. 

4) On s'interesse aux solutions stationnaires telles que R± = R 2 = R sw - 
Caculer la valeur d'equilibre R sw - 

5) Dans le cas ou (R\,R 2 ) = (R sw , R sw ), decrire revolution spatio-temporelle 
de la fonction u(x, t) correspondante. Illustrer cette description a l'aide 
de quelques traces schematiques. 

6) On s'interesse aux petites perturbations autour de l'un des trois equilibres 
(Rtw, 0), (0, Rtw) et (R sw , Rsw)- Ecrire le systeme linearise autour de l'un 
quelconque de ces equilibres. 

7) Calculer la valeur critique v = u c au-dela de laquelle les regimes "ondes 
propagatives" et "onde stationnaire" de la fonction u(x, t) echangent leur 
stabilite. 



CORRIGES DES EXERCICES ET PROBLEMES 
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CORRIGES DES EXERCICES ET PROBLEMES 



Corrige 3.1 Bifurcation noeud-col analytique 



l)Le premier cas est obtenu pour a fi > 0, c'est-a-dire a < avec \x < ou 
a > avec fx > 0. Dans ces deux cas, l'ensemble des equilibres est vide. On 
s'attend done a des solutions x(t) allant de — oo a +oo ou de +oo a — oo. Dans 
ce cas, on pose b = y/n/a. et le systeme s'ecrit 



x 



H + ax = n 



{x/bf + 1 



(3.20) 



On en deduit que arctg (x/b) = \x 6i+arctg (xo/t) ou x(0) = xq est la condition 
initiale. La solution s'ecrit alors 



x(t) 



tg 



arctg 



xq\ 
b J 



tg m 



■ML 
b 



(3.21) 



On peut voir l'ensemble des solutions comme etant la famille des courbes 



x(i) = b tg 



f(t-«d) 



(3.22) 



ou to est un reel quelconque, c'est-a-dire l'ensemble des translatees en temps de 
la courbe particuliere btg (xt). Par rapport a la condition initiale x(0) = xo, 
a t = on a to = — (6//x)arctg (xo/b). Ces trajectoires vont de — oo a +oo 
pour n > (et done a < 0) et de +oo a — oo pour \x < (et done a > 0). On 
note d'autre part que les trajectoires mettent un temps fini T = nb/fi pour 
relier les deux infinis. 

2)Le second cas est obtenu pour a [i < 0, c'est-a-dire a < avec ^ > ou 
a > avec /x < 0. Dans ces deux cas, l'ensemble des equilibres est {—a, a} 
avec a = yS—fi/a. On s'attend done a des solutions x(i) allant de — oo a a, 
de — a a a ou de a a +oo. Dans ce cas, le systeme s'ecrit 



x 



[i + ax 2 = a(x — a)(x + a) 



(3.23) 



Un decomposition en elements simples permet d'ecrire cette equation sous la 
forme 

— ( — —) = • (3.24) 

2aa \x — a x + a J 

L 'integration en temps conduit au logarithme d'une valeur absolue qui conduit 
a l'expression 

(3.25) 



x-a x -a 2aat xp-a _ 2 a t 
e = e a 



x + a xq + a 



x + a 
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en remarquant que era = Li/a. On en deduit la solution 

x(t) _ ^(e 2 ^ + l) + (e 2 ^-l) 
a ^( e 2^_i) + ( e 2^ + i) 

En effectuant quelques transformations, la solution s'ecrit finalement 

^-th(^) 



(3.26) 



x(t) = a 



1 - ^th (±t) 

a V a > 



(3.27) 



On peut classer ces solutions en trois families. Si xo £ [—a, a], la solution x(t) 
issue de cette condition initiale appartient a la famille des courbes 



x(t) 



-a th 



(3.28) 



oil to est un reel quelconque, e'est-a-dire l'ensemble des translatees en temps de 
la courbe particuliere —a th (^i) . Par rapport a la condition initiale x(0) = xo 
a t = on a to = (a//x)argth (xq/cl). Ces trajectoires vont de —a a +a pour 
/i > (et done a > 0) et de +a a —a pour fx < (et done a < 0). 
Si xo > a ou xo < — a , les solution x(t) issues de ces deux types de conditions 
initiales appartiennent a la famille des courbes 



x(t) = — ath 



/i 



(t - to) 



(3.29) 



avec t > to si x o > a e t t < to pour xo < —a. Comme to est un reel 
quelconque, ces courbes sont l'ensemble des translatees en temps de la courbe 
particuliere —ath (^t). Par rapport a la condition initiale x(0) = xo a t = 
on a to = — (o//x)argth (a/xo). Ces trajectoires vont de +oo a +a ou de —a 
a — oo pour \i > (et done a < 0) et de — oo a —a ou de +a a +oo pour 
(i < (et done a > 0). Pour ces trajectoires, to est aussi le temps fini qu'il 
faut pour atteindre Pinfini a partir de la condition initial xo- 
3) On verifie que les diagrammes de bifurcation obtenus pour a/j, < ou 
a/j, > correspondent bien au comportement des solutions analytiques qui 
vient d'etre decrit. 



Corrige 3.2 Bifurcation fourche analytique 



l)Dans le cas a fx > 0, on sait que est le seul equilibre et que les trajectoires 
joignent et l'infini. On pose alors 



x = fix + ax = n 



x x 

b \ ¥ 



+ i 



(3.30) 
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avec b = \J fi/a. Un decomposition en elements simples permet d'ecrire le 
systeme sous la forme 



x / 1 x/b 
b \xjb ~ l + x 2 /b 2 



li . (3.31) 



L 'integration de ce systeme conduit a des logarithmes de valeurs absolues puis 
a 

x(t) 1 1 

b ^/(l + 6 2 /x§) exp(-2 / ui) - 1 ^exp[-2/i (t - t )] - 1 

avec to = ^Ln(l + b 2 /xq). Cette famille de trajectoires est constitute des 
translatees en temps des courbes x(t) = 6/^/exp(— 2fit) — 1 definies pour t < 
lorsque /x > (et done a > 0) et pour t > lorsque lorsque // < (et done 
a < 0). Ces trajectoires relient l'infini et l'etat xo en un temps fini egal a to- 
2)Dans le cas a < 0, on sait que l'ensemble des equilibres est {— a, 0,a} 
avec a = ^—fi/a. On pose alors 

x = fix + ax 3 = a(x + a)(x — a)x . (3.33) 

Un decomposition en elements simples permet d'ecrire le systeme sous la 
forme 



+ 2- = 2a z a . (3.34) 

x + a x — a a 

En remarquant que a 2 a = —fi, l'integration de ce systeme conduit a des 
logarithmes de valeurs absolues puis a 

^ = ±^=^^^^ (3.35) 
a ^(a 2 /x 2 -l)exp(-2fit) + l 

On n'etudie alors que le cas x(t) > 0, le cas x(t) < s'en deduisant par 
symetrie. 

On peut classer ces solutions en deux families. Pour \xo\ < a les solutions 
s'ecrivent sous la forme 

x(t) = b/yfexp[-2[i(t-to)] + l (3.36) 

avec to = (l/2/x)Ln (a 2 /xg — 1). Ce sont les translatees en temps de la trajec- 
toire 6/Vexp(-2^it) + 1 qui est definie pour t G M. Ces trajectoires vont de 
a a pour /j, > (et done a < 0) et de a a pour fi < (et done a > 0). 
Pour IxqI > a les solutions s'ecrivent sous la forme 



x(t) = b/yjl-exp[-2fi(t-t )] 



(3.37) 
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QUESTIONNAIRES A CHOIX MULTIPLES 



avec to = (l/2/x)Ln (1 — a 2 /xq). Ce sont les translatees en temps de la trajec- 
toire bj \/l — exp(—2pt) qui est definie pour t > si p > (et done a < 0) et 
pour t < si p < (et done a > 0). Ces trajectoire relient l'etat xq et l'infini 
en un temps qui est egal a to- 

3) On verifie que les diagrammes de bifurcation obtenus pour ap < ou 
ap > correspondent bien au comportement des solutions analytiques qui 
vient d'etre decrit. 



Corrige 3.2 Travelling or Standing waves 



l)Dans un premier temps, on pose zi = pi(t)e ldl ^ pour n = 1,2. On 
obtient alors pi = p pi — a pf — 7 p 2 m pi et 9\ = u — (3 pf — 5 pf p m . On 
adimensionne les equations en posant t = t n T et p\ = p n Ri ou t n et p n 
sont des unites arbitraire. En choisissant t n = p et p n = \fpja et en 
posant v = 7/a, on obtient les equations indiquees qui decrivent revolution 
des modules p\ et p2- L'evolution des phases 9\ et 62 en decoule simple- 
ment. 2) On a Rt w = 1. 3)Ce cas correspond a une onde progative qui se 
deplace vers la droite (Travelling wave). 4) On a R sw = + v. 5)Ce 

cas correspond a une onde stationnaire (Standing wave). 6)La jacobienne de 
l'application non lineaire autour d'un point [Rx^R^) est donnee par la matrice 
\-ZR\-vR\ -2vR x R 2 \ ++ f-2 



-2v Rl R 2 \-ZR\-vR\y Cette matrice S ' &rit V 1 

autour de Pequilibre [Rtw, 0) = (1, 0) et ^ ^ ^ ) au t° ur de Pequilibre 

(R SW ,R SW ) = -^== (1, 1). Les valeurs propres des ces matrices sont respec- 
tivement (-2, 1 - v) et [-2(1 + 1/), -2(1 - 1/)]. 7)Pour < 1/ < 1, les ondes 
propagatives (tw) sont stables et les ondes stationnaires (sw) sont instables. 
C'est le contraire pour v > 1. On a done v c = 1. 



QUESTIONNAIRES A CHOIX MULTIPLES 

Bifurcation noeud-col 



QCM 3.1 



1) Le systeme dynamique x = p + a x 2 pour a < et /x > admet 
un equilibre stable et un equilibre instable 
aucun equilibre 



A 



B 



C I trois equilibres dont deux stables 
2) Le systeme dynamique x = p + a x 2 pour a < et p < admet 



A un equilibre stable et un equilibre instable 



APM-INPT thu-bifgen (2003), O. Thual November 14, 2009 



103 



B aucun equilibre 



C | trois equilibres dont deux stables 
3) Le systeme dynamique x = fi + a x 2 pour a > et /i > admet 
un equilibre stable et un equilibre instable 
aucun equilibre 



A 



B 



C trois equilibres dont deux stables 



QCM 3.2 



Bifurcation fourche 

1) Le systeme dynamique x = fj, x + a x 3 pour a < et fi > admet 
trois equilibres dont deux instables 
aucun equilibre 



A 



B 



C 



trois equilibres dont deux stables 
2) Le systeme dynamique x = [i x + a x 3 pour a < et fi < admet 
trois equilibres dont deux instables 
un seul equilibre stable 
trois equilibres dont deux stables 



A 



B 



3) Le systeme dynamique x = /x x + a x 3 pour a > et fi < admet 
trois equilibres dont deux instables 
un seul equilibre stable 



B 



C trois equilibres dont deux stables 



QCM 3.3 Bifurcation de Hopf 



1) Le systeme dynamique z = (fi + ioj) z + (a + i (3) \z\ 2 z pour a < et 
H > admet 



A 



B 



un equilibre instable et un cycle limite 

un seul equilibre stable 

un equilibre stable et un cycle instable 

2) Le systeme dynamique z = (/i + iw) z + (a + i [5) \z\ 2 z pour a < et 
H < admet 



C 



A un equilibre instable et un cycle limite 
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B un seul equilibre stable 

un equilibre stable et un cycle instable 



C 



3) Le systeme dynamique z = (fi + iuj) z + (a + i (5) \z\ 2 z pour a > et 
(i < admet 



B 



un equilibre instable et un cycle limite 
un seul equilibre stable 



C un equilibre stable et un cycle instable 



CORRIGES DES QCM 

Bifurcation noeud-col 



Reponses 3.1 



1)A 2)B 3)B 



Reponses 3.2 



Bifurcation fourche 



1)C 2)B 3)A 



Reponses 3.3 Bifurcation de Hopf 



1)A 2)B 3)C 
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OBJECTIFS PEDAGOGIQUES 

Cet article pedagogique est une introduction a l'etude des systemes dyna- 
miques simples, c'est-a-dire aux families de trajectoires des systemes d'equa- 
tions differentielles ordinaires couplees. Plusieurs objectifs de formations sont 
vises : 

• Assimiler la notion de stabilite d'un equilibre a partir de la linearisation 
des equations. 

• Comprendre la determination de la topologie des trajectoires au voisi- 
nage de l'equilibre a partir du systeme linearise. 
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NOTATIONS 



• Se familiariser avec Petude parametrique de la stabilite d'un equilibre 
en faisant varier un parametre. 

• Decouvrir ou consolider un certain nombre de concepts de bases comme 
par exemple la definition des systemes hamiltoniens. 

• Explorer quelques approches asymptotiques comme la derivation d'e- 
quations d'amplitudes a l'aide de la methode des echelles multiples. 

• Motiver l'exploration numeriques de la dynamique d'un certain nombre 
de modeles simples issus de la physique. 

Les competences a acquerir lors de Petude des cet article pedagogique sont 
les suivantes : 

• Etre capable de lineariser un syteme dynamique autour d'un equilibre, 
ce qui reveint a calculer la jacobienne d'un champ de vecteurs. 

• Etre capable de decrire la topologie des trajectoires autour d'un equilibre 
en determinant les valeurs propres et vecteurs propres de la jacobienne. 

• Etre capable de tracer le portrait de phase d'un systeme dynamique 
dependant d'un potentiel. 

• Connaitre la definition d'un systeme dynamique hamiltonien ainsi que 
ses proprietes de base (conservation de Penergie et des aires). 

• Etre capable de calculer une equation d'amplitude a l'aide de la methode 
des echelles multiples. 

• Etre capable d'effectuer une exploration numerique d'un systeme dy- 
namique quelconque a l'aide d'un schema numerique simple. 



Le niveau requis pour la lecture de cette article pedagogique se situe autour 
de celui d'une Licence scientifique. La lecture de "P Article Pedagogique Mul- 
timedia" intitule "Bifurcations generiques de l'equilibre" peut etre utile pour 
situer les bifurcations observees dans un cadre general. Elle n'est cependant 
pas indispensable. 



PRE-REQUIS 



NOTATIONS 



A(t) 



a 



Notation pour y 'fi/a ou y / — fx/a 
Amplitude complexe 
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(A, B, C) Composantes de U dans une base propre 

b Parmetre de controle du modele de Lorenz 

x(t) Derivee de la fonction x(t) 

DF(X) Jacobienne du champ de vecteur F(X_) 

div Operateur divergence d'un champ de vecteur 
Derivee partielle par rapport a t 
Derivee partielle par rapport a T 

^ Derivee partielle par rapport a q 

^ Derivee partielle par rapport a p 

E Energie d'un systeme hamiltonien 

Eq Valeur constante de l'energie 

/ Constante pour le terme de forcage 

F Constante pour le terme de forcage 

F(n; 20 Fonction vectorielle quelconque 

g Gravite (m s -2 ) 

G(X) Champ de vecteur 

h Parmetre de controle 

H(q,p) Hamiltonien 

k Parametre de controle 

L(t) Matrice 2x2 a coefficients periodiques 

M Constante 

0(e) Du meme ordre de grandeur que e 

p Variable de l'Hamiltonien H(q,p) 

p Variables pi de l'Hamiltonien H(q,p) 

q Variable de l'Hamiltonien H(q,p) 

q Variables % de l'Hamiltonien H(q,p) 

q e Valeur d'equilibre de q 

r Parmetre de controle du modele de Lorenz 

r c Valeur critique de r 

r* Valeur critique de r 

rd Valeur critique de r 

R Module de Z(0) 

R Matrice 2x2 

s Valeur propre 

(s a ,Sb,s c ) Valeurs propres 

si + \ s[ + \ si + ^ Valeurs propres pour 

S Symetrie point x — ► — x 

S Matrice 3x3 des composantes d'une symetrie 

t Temps (s) 

T e Periodicite (s) 

u(i) Variable d'un syteme dynamique 
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NOTATIONS 



u (t,T) 


Ordre zero de la methode des echelles multiples 


U!(t,T) 


Ordre un de la methode des echelles multiples 


u m 


Amplitude complexe 


u{t) 


Derivee seconde de u(t) 


U = (u, v, w) 


Vecteur dans le voisinage d'un equilibre 


U = (u,v) 


Vecteur dans le voisinage d'un equilibre 


u 2 


Norme au carre U ■ U_ 


V m 


Amplitude complexe 


V r 


Partie reelle du vecteur propre <f> 




Partie imaginaire du vecteur propre <\> 


V{u) 


Potentiel 


Vn{u) 


Derivee seconde de V{u) 


(x,y,z) 


Variables d'un systeme dynamique 


(x e ,y e ,z e ) 


Equilibres d'un systeme dynamiques 


X = (X 1 ,...) 


Vecteur de composantes X; 


xit) 


Trajectoire solution d'un systeme dynamique 


x e 


Equilibre d'un systeme dynamiques 


x + 


Equilibre particulier du modele de Lorenz 


x_ 


Equilibre particulier du modele de Lorenz 


x b 


Equilibre bas du pendule 


x h 


Equilibre haut du pendule 


X n 


Estimation numerique de X au temps nAt 


x',x",... 


Estimations numeriques de X 




Trajectoire solution d'un systeme dynamique 


z = x + iy 


Variable complexe d'un systeme dynamique 


Z(t) 


Amplitude complexe 


z* 


Complexe conjugue de Z 


\z\ 


Module de Z 


a 


Constante 


P 


Constante 


7 


Constante 


S 


Constante 


A 


Discriminant d'une equation du second degre 


A 


Constante oo 2 /w 


At 


Pas de temps d'un methode numerique 


e 


Petit parametre 


A 


Parametre constant 


A 


Parametre de frottement pour le pendule amorti (s^ 1 ^) 


a 


Parametres de controle 


p 


Module de z 


a 


Parmetre de controle du modele de Lorenz 


<P 


Argument de Z(0) 



Ill 



<p Vecteur propre reel ou complexe 

4>* Complexe conjuge de (f> 

((f) , 4> b , <f> c ) Vecteurs propres 

uj Parametre de controle des modeles de pendule 

uj Partie imaginaire de s (s^ 1 ) 

uj e Valeur particuliere de co 



COURS ECRIT 
Introduction 

L'objet de ce chapitre est de donner quelques elements utiles pour la compre- 
ension des instabilites observees en mecanique des fluides. 
L'exemple du modele de Lorenz est instructif a plusieurs titres. II permet tout 
d'abord de montrer que l'etude de stabilite d'un equilibre se ramene au calcul 
des valeurs propres de la matrice du systeme linearise. Ensuite, l'examen de 
la premiere bifurcation de ce modele est coherent avec le fait que la "bifur- 
cation fourche" est la bifurcation stationnaire generique lorsqu'une symetrie 
est brisee. Enfin, le comportement chaotique observe pour certaines valeurs 
des parametres permet de mettre en evidence la sensibilite aux conditions ini- 
tiales dans un systeme deterministe, ce que Ton appelle l'effet papillon dans 
le contexte de la meteorologie dont est d'ailleurs issu ce modele. 
L'exemple du pendule simple permet d'introduire la famille des systeme dy- 
namique dependant d'un potentiel qui est un cas particulier de la famille 
des systeme dynamiques hamiltoniens. Comme l'energie et le volume dans 
l'espace des phases est invariant, ces systemes sont appeles systeme conser- 
vatifs. lis presentent la particularite de ne pas posseder d'attracteurs, c'est- 
a-dire de sous-ensembles de l'espace des phases vers lesquels convergent les 
trajectoires de leurs voisinages. Au contraire, les systemes dissipatifs ont la 
propriete de pouvoir exhiber de tels attracteurs. 

La presentation suivante etudie les autres exemples classiques d'oscillateurs : 
pendule amorti, pendule resonant, oscillateur de Van der Pol et oscillateur 
parametrique. Ces exemples permettent de mettre en evidence un certain 
nombre de phenomenes que Ton retrouve dans de nombreuses systeme physi- 
ques sujets a des oscillations : portraits de phases en forme d'ellipses conver- 
gentes ou divergentes, phenomene de resonance, bifurcation de Hopf, pheno- 
mene de resonance parametrique, exposants de Floquet, etc. La determination 
de l'amplitude des oscillations bifurquees au voisinage du seuil d'instabilite 
est effectuee a l'aide d'une methode asymptotique aux echelles multiples pour 
l'oscillateur de Van der Pol et l'oscillateur parametrique. 
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1 Le modele de Lorenz et les systemes dissipatifs 

Le modele de Lorenz est l'un des systemes dynamique les plus simples per- 
mettant d'illustrer la competition entre un forgage exterieur et la dissipi- 
sation interne d'un systeme physique, tout en exhibant des comportements 
dynamiques non triviaux. 

En 1963, Lorenz a etudie numeriquement un systeme de trois equations 
differentielles cense representer tres grossierement la convection thermique 
dans P atmosphere [E. N. Lorenz, J. Atmos. Sci. 20 (1963) 130]. La demarche 
permettant de passer des equations fluides a ce systeme est exposee dans de 
nombreux ouvrages, en plus de l'article original de Lorenz et n'est pas abordee 
ici. Seul les regimes dynamiques de ce modele sont exposes ici. 

1.1 Les equations du modele de Lorenz 

On se contente ici d'enoncer que modele de Lorenz est le systeme dynamique 
qui s'ecrit 

x = —a x + a y 
y = —x z + r x — y 

z = x y — b z (4-1) 

avec X(t) = [x(t),y(t), z(t)} G M 3 , r > 0, a > et b > 0. 
L'espace des phases, constitue des vecteurs 2L £ -K 3 , est de dimension 3. On 
dit aussi que le systeme dynamique possede trois degres de liberte. 
L'espace de controle est constitue des vecteurs \x = (r, a, b) G M+. Le choix 
de Lorenz dans son article de 1963 est a = 10 et b = 8/3 avec r variable. 
Physiquement, r est proportionnel au gradient thermique vertical impose au 
fluide, a au nombre de Prandtl et 6 a l'elongation de la boite contenant le 
fluide. La solution d'equilibre trivial x = y = z = 0du systeme correspond 
physiquement a un regime ou le fluide est au repos et ou la chaleur se transmet 
uniquement par diffusion moleculaire. C'est Petat de base conductif. On va 
voir que pour r suffisamment grand, cet equilibre est instable et qu'il laisse la 
place a des regimes ou le transfert de chaleur est realise par diffusion et par 
convection. 

1.2 Equilibres du modele 

Le modele de Lorenz est done un systeme dynamique de la forme X = 
E.(fJ-',2Q- II est autonome dans la mesure ou la fonction F ne depend pas du 
temps. Les equilibres, points fixes, points critiques ou encore solutions sta- 
tionnaires de ce systeme sont les points X_ e = (x e , y e ,z e ) tels que F(X e ) = 0. 
La resolution du systeme F(X e ) = conduit au tableau suivant de solutions : 



1 . LE MODELE DE LORENZ ET LES SYSTEMES DISSIPATIFS 113 



< r < 1 : Une solution 


1 < r : Trois solutions 


2L e = 




x e e{x_,x + } 



avec X__ = (—a, —a, r — 1), X+ = (a, a, r — 1) et a = ^(r — 1). 
La stabilite d'un equilibre X_ e s'obtient en etudiant les valeurs propres de la 
jacobienne DF{X_ e ) obtenue en linearisant le systeme autour du point X_ e . Le 
calcul des composantes dFi/dXj de cette jacobienne conduit a 

/ -a a \ 
RF(X e ) =[r-z e -1 -x e . (4.2) 

V Ve X e -b J 



1.3 Stabilite de l'equilibre nul 



La stabilite de l'equilibre trivial X_ e = s'obtient en etudiant les valeurs 
propres s de la matrice 



(4.3) 



—a 


a 


r 


-1 





-b 



QF(0) = 

dont 1 'equation caracteristique est 

(s + b) \s 2 + (<r + l)a + a(l - r) 



= . 



(4.4) 



Une premiere valeur propre est s = —b. Le discriminant de Pequation du 
second degre permettant de determiner les deux autres est A = {a — l) 2 +4ar. 
Comme il est positif, il y a deux racines reelles de somme — (cr+1) et de produit 
<t(1 — r). On peut done resumer le signe des trois valeurs propres reelles de 
DF(0) sur le tableau suivant : 



< r < 1 : est stable 


1 < r : est instable 


signe de s 6 {— , — , — } 


signe de s G { — , — , +} 



Lorsque r < 0, les trois valeurs propres sont negatives : l'equilibre X_ e = est 
done stable. Lorsque r > 0, une des valeurs propres est positive : l'equilibre 
est done instable. Pour r = 1, la stabilite est marginale. 

Notons Sb = —b la valeur propre independante de r, s c (r) la valeur propre 
toujours negative et Sd(r) la valeur propre qui devient positive pour r > 1. 
Les vecteurs propres associes sont notes <p b , <f> (r) et <fi d (r). Ces vecteurs sont 
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sigma=10, b=2.67, r=0.3 sigma=10, b=2.67, r=0.3 




Figure 4.1: Trace de particule pour r = .3. a) t=0, b) t=.3, c) t=1.2. 
L'equilibre est stable. 

dermis a une constante pres et on peut par exemple choisir eft = (0,0,1) 
pour le vecteur propres associe a la valeur propre s& = —b. L'expression 
analytique des autres valeurs propres et vecteurs propres est triviale mais 
n'est pas explicitee ici par souci de concision. 

1.4 Portraits de phase au voisinage de l'equilibre nul 

Si nous notons U(t) les trajectoires X_(t) = + ILif) proches de 0, l'approxi- 
mation de linearisation indique qu'elles sont solutions du systeme dynamique 
lineaire U_(t) = DF(0)U_ ce que Ton peut expliciter de la maniere suivante : 

it = —o u + <7 v 
v = r u — v 

w = —bw, (4-5) 

oil l'on a note U = (u,v,w). Par rapport au systeme initial, on voit que l'on 
a neglige les termes non lineaires — x z et x y. On peut alors poser 

U(t) = Bit) ^ + C(t) c + D(t) l d (4.6) 

ou (B, C, D) sont les composantes de U dans la base propre de la matrice 
DF(0). En reportant cette decomposition dans le systeme lineaire ci-dessus 
gouverne par cette matrice on obtient trivialement les equations 

B = s b B , C = s c C , D = s b D , (4.7) 

ce qui se resout en 

B(t) = B{0)e Sbt , C(t) = C{0)e Sct , D[t) = D{0)e Sdt . (4.8) 

Ce resultat demontre le fait que est stable lorsque les trois valeurs propres 
sont negatives et qu'il suffit d'une valeur propre positive pour obtenir une 
divergence des trajectoires dans la direction instable. 
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sigma=10, b=2.67, r=1 




y 



Figure 4.2: Vecteurs propres de la statibilite de pour r = 1 : ck, <f> (1) et 

II permet aussi de tracer le portraits de phase, c'est-a-dire l'ensemble de tra- 
jectoires, dans le voisinage de 0. Le cas le plus interessant a examiner est 
obtenu pour r proche de la valeur critique r = 1. Pour cette valeur, les 
vecteurs propres associes aux valeurs propres (sb,s c ,Sd) = [— b, — (a + 1),0] 
sont 




Dans un plan affine engendre par les vecteurs [^,(1), <j) (1)], les trajectoires 
convergent toutes vers l'intersection de ce plan de coordonnees (B, C) avec 

tT + l 

l'axe 4> le long des courbes d'equations C/C(0) = [B/B(0)] b . Cette con- 
vergence est exponentielle. 

On voit done que toutes les trajectoires convergent vers l'axe engendre par 
au moins pour celles qui sont proches de 0. On montre que e'est en fait 
le cas pour toutes les trajectoires, meme eloignees de 0, mais il faut utiliser 
pour cela un raisonnement (fonctionnelle de Liapounov) qui depasse le cadre 
de cette presentation. 

Pour r legerement inferieur a un, l'orientation des vecteurs propres est legere- 
ment modifiee. En revanche, les trajectoires sur l'axe 4> d convergent vers zero, 
alors qu'elles etaient reduites a un continuum de points (tout l'axe) pour r 
exactement a un. Toutes les trajectoires du systeme lineaire convergent done 
vers qui est done stable. 

Pour r legerement superieur a un, les trajectoires divergent, lentement, sur 
l'axe 0,. II arrive alors un moment ou 1' approximation lineaire n'est plus 
valable, si bien que ces trajectoires quittent l'axe 0,. II est naturel de penser 
que ces trajectoires vont converger vers l'un des equilibres bifurques a partir 
de la valeur critique r = 1. II convient maintenant d'etudier la stabilite de 
ces nouveaux equilibres. 
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Figure 4.3: Valeurs propres en fonction de r £ [0, 2]. a) Stabilite de : Sd(r), 
s b {r) et s c (r) b) Stabilite de X W : s ( +\r), s[ +) (r) et si + \r) 
c) Superposition des deux graphes. 



1.5 Stabilite des equilibres bifurques 

La stabilite de l'equilibre X e = X + = (a, a, r— 1) avec a = \Jb(r — 1) s'obtient 
en etudiant les valeurs propres s de la matrice 



DF(X + ) =1 -1 -o (4.10) 




dont l'equation caracteristique est 

s 3 + s 2 (a + b + 1) + s b{a + r) + 2b(r - l)a = . (4.11) 



Suivant les valeurs des parametres de controles (r, a, b) , ce polynome de degre 
trois peut avoir trois racines reelles ou bien une racine reelle et deux racines 
complexes conjuguees. Notons s^~\r) , si + \r) , s^\r) ces trois racines. Le 
signe de la partie reelle des trois racines de X_ e determine la stabilite du point 
d'equilibre X e . Sans entrer dans une discussion parametrique complete sur 
cette stabilite, nous allons nous concentrer ici sur les cas ou r est dans le 
voisinage de la valeur critique r = 1. 

Comme X_ + = pour r = 1, on peut faire le lien avec les trois valeurs propres 
determinant la stabilite de et ordonner les notations de telle sorte que Ton 
ait : 

si + \l) = -b, 4+)(l) = -( CT + l) et 4 +) (l) = 0- (4.12) 

En posant r = 1 + e ou e > est un petit parametre, on peut done effectuer 
les developpements limites suivants : 

4 +) (l + e) = -b + l3e + 0{e 2 ) 

+ = -(tr + l)+7e + 0(e 2 ) 
4+)(l + e ) = 5e + 0(e 2 ). (4.13) 
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Figure 4.4: Valeurs propres en fonction de r £ [0,30]. Stabilite de X^' 

4 +) W> 4 + V) et si +) (r) 



igma=10, b=2.67, r=1 .5 



sigma=10, b=2.67, r=1 .5 



=10, b=2.67, r=1 .5 




Figure 4.5: Trace de particule pour r 
L'equilibre X_( + ^ est stable. 



1.5. a) t=0, b) t=3, c) t=6. 



Pour e suffisamment petit, les valeurs propres s[ + ^ et restent reelles et 
negatives. La stabilite de X e est depend done du signe de et done du signe 
de 5. En reportant le developpement de dans le polynome caracteristique 
et en identifiant l'ordre dominant en e, on voit facilement que 5 = —2a/{a+\). 
On peut done affirmer que X + est stable lorsque r est dans un voisinage 
de la valeur critique r = 1. L'etude de stabilite de X__ conduit au meme 
polynome caracteristique et done a la meme conclusion. En fait, la symetrie 
S qui a X = (x, y, z) associe S_X_ = (— x, —y, z) est une symetrie du systeme 
dynamique. En effet, on voit que F(S_X_) = si bien que si X_(t) 

est un solution, alors Y^(t) = S2L(t) en est une autre. On remarque que 
l'equilibre X_ e = verifie S_X_ e = X_ e , tandis que les deux equilibres bifurques 
a partir de la valeur r = 1 verifient S X + = X_ et S2L- = X + . On dit 
que la bifurcation obtenue pour r = 1 brise la symetrie S_. L 'ensemble de 
ces considerations montre que l'on est en presence d'une bifurcation fourche 
supercritique. 

Les trajectoires dans le voisinage de l'equilibre X, s'obtiennent en tragant 



,(+) 
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tout d'abord les vecteurs propres de la matrice DF{X_ Jr ). Pour r proche de 
1, ces vecteurs propres sont proches des vecteurs propres </>,, </> c (l) et ^ rf (l) 
que l'on a deja explicites pour etudier les trajectoires autour de 0. Dans 
un plan affine engendre par cp b et 4> c , les trajectoires sont identiques a celles 
obtenues dans le voisinage de 0. Sur l'axe <j> les trajectoires convergent 
vers l'equilibre X_ + . Le portrait de phase dans le voisinage de X- et pour 
r toujours proche de 1, s'obtient en tragant les trajectoires par symetrie par 
rapport aux trajectoires du voisinage de X_ + . 

Pour r proche de 1, on peut reunir les trajectoires du voisinage de et les 
trajectoires des voisinages de X + et X__ en tracant une trajectoire partant 
de dans la direction du vecteur propre de DF(0) proche de <p d et rejoignant 
X + le long de la direction propre de DF(X + ) proche de 4> d . La trajectoire 
symetrique relie et X_ + de la meme maniere. 



1.6 Autres regimes lorsque r devient grand 

Lorsque r augmente au-dela du voisinage de r = 1, on peut chercher s'il existe 
une valeur critique r = r* pour laquelle les equilibres X_ + et JC__ deviennent 
instables (simultanement a cause de la symetrie). On rappelle que la stabilite 
de ces equilibres depend du signe de la partie reelle des racines du polynome 
caracteristique : 

s 3 + s 2 (cr + b + 1) + s b(o + r) + 2b{r - l)o = . (4.14) 

La destabilisation de ces equilibres par changement de signe d'une valeur 
propre reelle est impossible car si s = on a forcement r = 1. On peut done 
rechercher pour quelles valeurs des parametres on peut obtenir deux racines 
iu et — iui de partie reelle nulle. En reportant la valeur s = iuj dans Pequation 
caracteristique, on obtient les deux conditions —to 2 (a + b+1) +26(r — l)a = 
et — iuo 3 + iujb(a+r) = 0. En eliminant lo 2 entre les deux equations, on obtient 
finalement 

_ ff (a + b + 3) 
r *~ ( a -b-l) ■ (4 ' 15j 

Pour les valeurs classiques a = 10 et b = 8/3, on obtient r* = 470/19 ~ 24, 73. 
On peut remarquer qu 'entre le voisinage de r = 1 ou les trois racines sont 
reelles, et le voisinage de r = r* ou il y a une racine reelle et deux racines 
complexes conjuguees, il existe une valeur r = telle que le polynome ad- 
mette un racine double negative. Pour r legerement superieur a r^, il y a deux 
valeurs propres complexes conjuguees de partie reelle negative et les portraits 
de phase dans le voisinages de 2L+ et 2L- sont des spirales elliptiques con- 
vergeant exponentiellement vers ces equilibres stables. 

La destabilisation de des equilibres X + et X_ correspond done a une insta- 
bility oscillatoire et done une bifurcation de Hopf. Deux valeurs propres com- 
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iigma=10, b=2.67, r=5 sigma=10, b=2.67, r=5 sigma=10, b=2.67, r=5 




Figure 4.6: Trace de particule pour r = 5. a) t=0, b) t=l, c) t=3. L'equilibre 
est stable. 



sigma=10. b=2.67, r=30 sigma=10, b=2.67. r=30 sign 




Figure 4.7: Trace de particule pour r = 30 > r*. a) t=0, b) t=.l, c) t=2. Le 
regime est chaotique. 

plexes conjuguees traversent l'axe des imaginaires lorsque le parametre r fran- 
chit la valeur critique r* . Dans le quart de plan (a, b) , il existe deux domaines 
connexes dans lesquels cette bifurcation est supercritique, le complementaire 
correspondant au cas oil la bifurcation est sous-critique. La determination 
analytique de ces domaines depasse le cadre de ce cours (calcul de formes 
normales). Indiquons simplement que les valeurs classiques a = 10 et b = 8/3 
correspondent au cas ou la bifurcation est sous-critique. 

Dans ce cas, la simulation numerique du systeme dynamique montre que l'on 
transite abruptement vers un regime chaotique ou la trajectoire tourne autour 
d'un des equilibres instables X_+ et 2L- comme si elles y convergeait avant 
de basculer aleatoirement vers l'autre equilibre pour y repeter le meme type 
de comportement. On observe que toutes les trajectoires convergent vers 
une trajectoire chaotique que l'on nomme "attracteur etrange". On montre 
que la dimension fractale de cette trajectoire est superieure a un et que la 
distance entre deux conditions initiales initialement tres proches s'amplifie 
tres rapidement. Cette sensibilite aux conditions initiales, ainsi que le chaos 
deterministe observe avec ce systeme dynamique simple, ont servi de base 
a ce que l'on appelle "l'effet papillon". Si revolution de 1' atmosphere est 
comparable a revolution d'un systeme dynamique comme celui de Lorenz, 



120 APM-INPT thu-sysdyn (2003), O. Thual November 14, 2009 



la prise en compte ou non de l'effet du battement d'une aile de papillon 
dans les conditions initiales peut changer radicalement revolution du temps 
(atmospherique) au bout d'un temps (en nombre de jours) estime a environ 
7 jours. 

2 Le pendule simple et les systemes conservatifs 

Le pendule simple est le systeme dynamique prototype des modeles conser- 
vatifs, qui decrivent des systemes pour lesquels le forgage exterieur et la dissi- 
pation interne sont nulles. C'est un cas particulier des systemes dynamiques 
dependant d'un potentiel qui sont eux-memes un cas particulier des systemes 
hamiltoniens. 

2.1 Equation du pendule simple 




S 



Figure 4.8: Pendule simple 

On considere un pendule simple forme d'une masse m assimilee a un point 
et suspendue a un fil de masse negligeable et de longueur I. En l'absence de 
frottement, Tangle u(t) que fait le pendule avec la verticale verifie Pequation 
differentielle ordinaire 

u + u? sinn = . (4.16) 

oil lo = \J gjl et g est la gravite. Cette equation ne fait pas intervenir la 
masse. 

Cette equation differentielle de degre deux peut se transformer en un systeme 
dynamique de degre unl= E_(X_) en posant X_ = (q,p) et en ecrivant 

q = p 

p = —J 1 sing. (4-17) 

L'espace des phases de ce systeme est l'ensemble des X_ 6 M 2 et l'espace de 
controle est constitue de l'ensemble des valeurs du parametre to 2 £ M + . 
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On remarque l'existence d'une invariance qui est la translation T^^X = X + 
(2vr,0). Si X(t) est une trajectoire, alors Y(t) = T 2jT X(t) = X(t) + (2vr,0) 
en est une autre. On dit aussi que l'espace des phases est Pensemble des 
X_ G T 1 x M oil T 1 est le cercle de perimetre 2ir. 

2.2 Stabilite de l'equilibre haut 

Les deux equilibres de ce systeme, obtenus en resolvant F(X e ) = 0, sont 
Xb = (0,0) et Xh = (^jO)- Ces equilibres correspondent a la position basse 
(Xb) e ^ a la position haute (Xh) du pendule. 
La jacobienne de l'application F est la matrice 

DF(X) = ( 2 ° v I) . (4.18) 

La stabilite de Xh es t done obtenue en calculant les valeurs propres de la 
matrice 

RF(K h ) = ^2 I) , (4-19) 

e'est-a-dire en resolvant l'equation s 2 - co 2 = 0. Les valeurs propres sont done 
s = u), a laquelle est associe le vecteur propre (p + = (1,lo) et s = —uo a laquelle 
est associe le vecteur propre §~ = (1, —uj). L'equilibre est done instable, la 
direction d'instabilite etant le long de l'axe 4> + . Le portrait de phase autour 
de X_ h est constitue d 'hyperboles dont les asymptotes sont donnees par les 
directions 4> + et <fi~ . 

2.3 Stabilite de l'equilibre bas 

La stabilite de Xb = Q est obtenue en calculant les valeurs propres de la 
matrice 

£*m)=(_° 2 I) (4-20) 

e'est-a-dire en resolvant l'equation s 2 +uj 2 = 0. Les valeurs propres sont done 
s = iuj, a laquelle est associe le vecteur propre (p = et s* = —ioj a 

laquelle est associe le vecteur propre (f>* = (1,—icu). 

L 'etude de la stabilite de Xb = Q revient a poser X(t) = Xb + U(t) e * a 
negliger les termes d'ordre U 2 . On se ramene done a l'etude de l'equation du 
pendule lineaire u + u 2 u = qui s'ecrit aussi sous la forme 




(4-21) 



On peut indiquer directement la solution de l'equation du second ordre en 
posant u(t) = a cos(ujt + tp) ou bien poser Z(t) = u + i^v pour combiner 
astucieusement les equations du systeme d'ordre un. 
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On obtient alors Z = iuZ qui se resout en Z(t) = Z(0)e iult = Re iu/t+i v. 
Comme u(t) = Z[t) + Z*(t) = 2Rcos{u>t + icp), ces deux approches se re- 
joignent en posant a = 2R. On a alors v(t) = u(t) = —a to sm(ut + i(p). Les 
trajectoires sont done des ellipses d'axes paralleles aux axes des coordonnees 
(u, v), engendrees par les vecteurs a(l,0) et a(0,u>). Le sens de rotation est 
contraire au sens trigonometrique. 

Une autre maniere de presenter la deuxieme approche consiste a poser 

EW = zW# + z- Wr = (.JW + fJ (() 

En reportant dans le systeme lineaire et en utilisant le fait que (j> et (ft* sont 
vecteurs propres, on obtient le systeme 

Z = iujZ et Z* = -iuZ*. (4.23) 

On remarque que les ellipses formees par les trajectoires sont engendrees 
par des couples de vecteurs proportionnels aux vecteurs V_ r = Re(0) et 
Y_i = Im((/)). On verifie facilement que cette propriete s'applique de maniere 
tres generale aux portraits de phase associes a des valeurs propres complexes 
conjuguees obtenues en etudiant la stabilite d'un equilibre. 






2.4 Systemes dependant d'un potentiel 

Le pendule simple est un cas particulier de la famille des equations differentielles 
s'ecrivant sous la forme u = —V'[u) ou V(u) est une fonction derivable quel- 
conque appelee potentiel. On peut poser par exemple V(u) = u 2 (l — cosu) 
pour le pendule non-lineaire et V(u) = \oj 2 u 2 pour le pendule lineaire appele 
aussi pendule harmonique. 
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On peut ramener ces equations a 
posant X_ = (q,p) et en ecrivant 



des systemes de degre unl = F(2Q en 



Q 
P 



P 



-V'(q) . 



(4.24) 



Les equilibres sont des points X e = (q e , 0) tels que X e est un extremum de la 
fonction V. La stabilite de ces equilibres est obtenue en calculant les valeurs 
propres de la matrice 



c'est-a-dire en recherchant les solutions de s 2 + V"(q e ) = 0. Si V"(q e ) < 0, 
la valeur q e correspond a un maximum du potentiel V. Les valeurs propres 
sont alors s = ±y / —V"(q e ) et l'equilibre est instable. Si V"{q e ) = uj 2 > 0, la 
valeur q e correspond a un minimum du potentiel V. Les valeurs propres sont 
alors s = ±iuo e et l'equilibre est marginal. 

En restant dans le cadre d'un potentiel V{u) quelconque, on considere la 
fonction E(t) = \p{t) 2 + V[g(t)] associee a une trajectoire [q(t),p(t)} solution 
du systeme. On montre facilement que E{t) = p(t)p(t) + q(t) V'[q(t)] = 
-p(t) V'[q(t)] +p(t) V'[q(t)] = 0. On appelle energie la quantite E(t) = E 
qui est done un invariant du mouvement. Cette energie s'ecrit aussi E{t) = 
|n 2 + V(u) et peut etre vue comme la somme de Penergie cinetique \ii 2 et 
de l'energie potentielle V{u). 




Figure 4.10: Syteme dependant d'un potentiel. a) Forme d'un potentiel V(q). 
b) Portrait de phase dans le plan (q,p). 

En utilisant l'invariance de l'energie on peut afnrmer que les trajectoires sont 
les courbes H(q,p) = \p 2 + V(q) = Eq. Grace a cette interpretation, on peut 
tracer le portrait de phase d'un systeme defini par une fonction V quelconque. 
II suffit pour cela d'imaginer une bille qui roule le long d'une colline dont 
l'altitude est le potentiel V. 




(4.25) 
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2.5 Systemes dynamiques hamiltoniens 

Les systemes dependant d'un potentiel sont des cas particuliers de la famille 
des systemes hamiltoniens qui s'ecrivent 

dH 

q = 

d FT 

P = -^-(9,P) (4-26) 
ou H(q,p) est une fonction derivable quelconque. 

En posant H(q,p) = \p 2 + V(q) on retrouve le cas des systemes dependant 
d'un potentiel. Mais on peut considerer d'autres formes de PHamiltonien 
H(q,p) comme par exemple H = q p. 

Etant donnee une solution [q(t),p(t)], on montre que la quantite E(t) = 
H[q(t),p(t)] est encore un invariant du mouvement. En effet, 

E_.dH .dH _dHdH dH dH () 
^ dq ^ dp dp dq dq dp 

On montre de plus que les aires sont conservees dans l'espace des phases. En 
effet, si on note 2L = F(2Q ay ec X_ = (q, p) le systeme dynamique hamiltonien, 
on peut calculer la divergence 

Cette propriete entraine qu'un volume (ici une aire) de conditions initiales 
reste constant au cours du temps, comme dans le cas d'un ecoulement isochore. 
On parle alors de systemes conservatifs. De tels systemes n'admettent pas 
d'attracteurs, c'est-a-dire d'ensembles vers lesquelles toutes les trajectoires 
d'un bassin d'attraction pourraient converger. 

L'invariance de l'hamiltonien H et des aires dans l'espace des phases se 
generalise au cas des systemes hamiltoniens H(q, p) dans M 2n dont l'expression 
generale est 

dH 

dpi ~ ~ 
d FT 

A = -^(1.2) (4.29) 
pour i = 1, n avec q G M n et p G M n . 



3 Autres oscillateurs et phenomenes non lineaires 

On considere maintenant des exemples d'oscillateurs qui prennent en compte 
les phenomenes de frottement ainsi que des forgages externes. Contrairement 
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aux systemes conservatifs, ces systemes sont dissipatifs. Dans certains cas, 
il s'etablit un equilibre (au sens large) entre le forcage et la dissipation. Cet 
equilibres sont parfois des attracteurs dans la mesure ou un ensemble de con- 
ditions initiales peuvent y converger. La particularite des exemples ci-dessous 
est le fait que ces equilibres se traduisent par des regimes oscillant en temps. 

3.1 Le pendule amorti 

Lors de l'oscillation du pendule, le frottement de Pair peut se modeliser par 
une force proportionnelle a la vitesse. On obtient done l'expression 

il + 2A ii + u? sin u = . (4.30) 

Le parametre A > est done un facteur d'amortissement du mouvement du 
pendule. 

Les equilibres sont les memes que pour le pendule non amorti. Leur stabilite 
et le portrait de phase de leur voisinage s'obtient en linearisant les equations 
autour de la position haute (instable) ou de la position basse (stable). 
Une maniere rapide d'etudier la stabilite de la position basse u = consiste 
a chercher des solutions complexes de la forme u{t) = u m e st avec u m G W et 
s G (J* dans l'equation linearisee autour de zero qui s'ecrit u + 2\ii + uj 2 u = 0. 
Comme l'equation est lineaire, la partie reelle d'une solution complexe est une 
solution. On doit alors chercher les racines de l'equation s 2 + 2\s + u 2 = 0. 
Lorsque A < to, il existe deux racines complexes conjuguees 

s = -\±iuy/l - X 2 /lu 2 . (4.31) 

Le retour a l'equilibre s'effectue par des oscillations exponentiellement amor- 
ties. Lorsque A est petit, les oscillations sont faiblement amorties et la pulsa- 
tion est proche de to. 

Lorsque A > u, il existe deux racines reelles 

s = -X (li^l-c^/A 2 ) . (4.32) 

Le retour a l'equilibre exponentiellement sans oscillation. On parle alors de 
"sur-relaxation" . 

L 'etude de stabilite autour de la position haute fait apparaitre uniquement des 
valeurs propres reelles dont l'une est positive. En connectant les portraits de 
phase des voisinages des deux equilibres, on peut tracer toutes les trajectoires 
de l'espace des phases M 2 de coordonnees (u, ii). 

3.2 Le pendule resonant 

On peut maintenant considerer le cas ou les oscillations du pendule sont en- 
tretenues en exergant sur la masse ponctuelle une force oscillant a la pulsa- 
tion uj e (en imposant par exemple un champ magnetique lorsque la masse est 
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Figure 4.11: portrait de phase du pendule amorti pour T = 2-k/uj = 1 et 
A = 1. 



aimantee). Les equations du mouvement peuvent alors s'ecrire 

u + 2Au + uj 2 sinu = f cos(u> e t) . (4.33) 

Cette equation differentielle ordinaire de degre deux peut se mettre sous la 
forme d'un systeme dynamique de degre unl = E_(X_,t) avec X_ G M 2 . Ce 
systeme n'est pas autonome dans la mesure ou F_ depend du temps. 
La position basse u = n'est plus un equilibre du systeme dynamique si 
/ 7^ mais on peut quand meme s'interesser aux trajectoires situees dans le 
voisinage de X_ = 0, surtout si le forgage / est un petit parametre. Dans ce 
cas, ces trajectoires sont solutions de l'equation linearisee 

u + 2Xu + J 2 u = fcos(uj e t) . (4.34) 

Sauf cas particulier discute ci-dessous, les solutions de cette equation lineaire 
(non homogene) sont de la forme u(t) = u m e lWe 1 + u* n e~ lUe 1 ou u m £ W et 
s e (F sont des coefficients complexes. En reportant dans l'equation lineaire 
on obtient 

(lo 2 -oj 2 e + 2i\oj e )u m = f '/2 . (4.35) 

En notant A = uj 2 /oo 2 et A = X 2 /uj 2 , l'amplitude de la reponse du pendule 
au forgage est done 

f 1 

(4.36) 



1 m| 2ul 7(1-A)H4AA • 

Cette reponse est maximale pour A = 1— 2A et s'ecrit alors \u m \ = 2^ 4A(i_a) • 
Lorsque A est tres petit (faible dissipation) la reponse est tres grande pour 
pour pour A ~ 1 e'est-a-dire oj e ~ uj. On est en presence du phenomene de 
resonance. 

Lorsque A = et uj = uj e = uj, e'est-a-dire pour le cas de la resonnance, on 
cherche une solution sous la forme u(t) = t {y m e tu)t + v^ n e~ lu ' t ) avec v m £ (E. 
On trouve alors v m = — i-^j ce qui conduit a la solution u(t) = ^ t s'm(ujt). 
La reponse du pendule croit done lineairement avec le temps a la resonance. 
Lorsque les trajectoires s'eloignent du voisinage de u = 0, 1' approximation 
lineaire n'est plus valable et il faut tenir compte de la nonlinearite contenue 
dans le terme sin u. La reponse du pendule lorsque A varie et pour differentes 
valeurs de / peut etre obtenu numeriquement ou par un developpement 
asymptotique lorsque / est petit. On peut alors observer le phenomene 
d'hysteresis correspondant au fait que deux reponses sont possibles pour une 
meme valeur des parametres de controle. 
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3.3 L'oscillateur de Van der Pol 

Une autre fagon de forcer le pendule peut etre obtenue en imposant une 
force proportionnelle a la vitesse et dans le meme sens, au moins pour les 
petites oscillations. On peut imaginer une telle force sur un pendule en util- 
isant l'induction magnetique que provoque la masse ponctuelle lorsqu'elle est 
electriquement chargee. Mais les equations du pendule avec une dissipation 
negative pour les petites oscillations apparaissent naturellement dans un mon- 
tage electronique portant le nom d'oscillateur de Van der Pol. Ces equations 
s'ecrivent 



u + (k u 2 — 2 /i) it + u> 2 sin u = 



(4.37) 



La stabilite de l'equilibre u = s'obtient en etudiant la partie reelle des 
valeurs propres du systeme linearise qui s'ecrivent s = /i ± iuj. La position 
basse u = est done stable pour fi < (on retrouve la notion de dissipation 
en posant [i = —A) et instable pour [i > 0. La valeur critique de la bifurcation 
est /j, = 0. On s'attend alors a etre en presence d'une bifurcation de Hopf. 
II est possible de calculer l'amplitude de saturation de l'instabilite par les 
termes nonlineaires au voisinage de la bifurcation /j, = 0, grace a des deve- 
loppements asymptotiques ou fi est le petit parametre. On peut par exemple 
effectuer un developpement en echelles multiples en posant /i = M e et 



u(t) = yfe Uo(t, e t) + e u\(t, e t) + o(e 2 ) 



(4.38) 



ou les fonctions Ui(t,T) sont considerees comme des fonctions de deux vari- 
ables independantes et doivent rester d'ordre 1. Ce choix de developpement 
est conforme a ce que Ton attend d'une bifurcation de Hopf ou une oscil- 
lation sur une echelle de temps d'ordre 1 coexiste avec une amplification 
sur une echelle de temps inversement proportionnelle a l'ecart au seuil et ou 
l'amplitude de Poscillation bifurquee croit comme la racine carree de l'ecart 
au seuil. 

II convient alors de calculer les expressions suivantes : 



u 



u 



smu 
u 2 
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dt \dT at 
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u 3 + o{u A ) = sfe 



uq + e tii 



-Ur 



+ 



4 2 ) 



e Uq + o(e) 



(4.39) 



En reportant dans l'equation, l'ordre dominant 



ce qui conduit a la solution u (t,T) = A{T)e lu}t + A*{T)e~ lu)t ou A(T) est un 
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fonction complexe arbitraire de la variable T. A l'ordre suivant, en e-^/e, on 
obtient P equation 

d 2 ui 2 d 2 u 2 <9u <9u 1 3 

n=3 

= c ™( r ) einu)t ( 4 - 4 °) 

n=— 3 

oil les fonctions C n (T) s'expriment en fonction de A(T) et de ses derivees. 
Pour que u\(t,T) reste d'ordre 1, comme suppose dans le developpement, 
il est necessaire d'imposer C\(T) = C_i(T) = 0. En effet, en l'absence de 
cette condition dite de compatibility, le terme u\ serait proportionnel a t et 
done non borne. On parle aussi de condition de resonance. En calculant 
l'expression de C\(T), cette condition s'ecrit done 

dA , , . ( k 3 .\ , 2 



iT MA -[S + ^ i )^ A - (4 - 41) 
En posant Z(t) = yfeA(et) e tujt , on obtient finalement la solution 

u(t) = Z(t) + Z*(t) + 0(\Z\ 3 ) (4.42) 

oil Pamplitude complexe Z(t) est solution de l'equation 

k 3 

Z = (fi + iu) Z - a\Z\ 2 Z avec «=- + —«. (4.43) 

3 2 oj 

On voit done que la bifurcation de Hopf est supercritique pour k > et 
sous-critique pour k < 0. 



3.4 L'oscillateur parametrique 

Une autre fagon d'exciter un pendule consiste a faire osciller sa longueur 
ou bien le champ de gravite apparent en faisant osciller verticalement son 
point d'attache. On note uj e la pulsation de ces oscillations de forgage dit 
parametrique dans la mesure oil il modifie les parametre de controle du pend- 
ule simple. Lorsque ce forcage est petit, l'equation du pendule parametrique 
s'ecrit : 

u + 2A it + u? [1 + h cos(w e t)\ sin u = . (4.44) 

L'espace de controle de ce systeme est constitue des parametres (A, h) € 
M\. En effet, un changement d'unite de temps permet de voir que seul le 
rapport Lo/uj e est pertinent pour l'etude des solutions. 

Cette equation differentielle ordinaire de degre deux se ramene a un systeme 
dynamique d'ordre un de la forme X = F(2L, t) avec X G M 2 qui est 
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dependant du temps. La stabilite de l'equilibre X_ = s'obtient en linearisant 
les equations ce qui conduit ici au systeme lineaire 

U = Ut) U avec Ut) = ( _^ „ ^ f) JJ. (445) 

L 'etude des systemes dynamiques lineaires J7 = L{t)U_ avec J7 € M n et oil 
L(i) est une matrice reelle a coefficients periodiques de periode T e est appelee 
probleme de Floquet. On montre que l'espace vectoriel des solutions com- 
plexes est engendre par n solutions de la forme 4>.{t) e Sj 1 pour j = 1, n ou 
les (j>j(t) sont des fonctions complexes periodiques de periode T e et les Sj des 
nombres complexes appeles exposants de Floquet. 

La demonstration de ce resultat est simple en considerant la matrice R telle 
que U[T e ) = R IZ(0). On peut construire numeriquement cette matrice en 
integrant le systeme dynamique lineaire de partir d'un ensemble de 

conditions initiales formant une base de lR n . Si on note r,j pour j = 1, ...,n 
les valeurs propres complexes de R, les coefficients de Floquet sont obtenus 
en posant Sj = ^Lnr^ pour j = 1, ...,n. La propriete U_(t + T e ) = RU_(t) 
permet de conclure en construisant les fonctions 4>-{t) a partir des solutions 
integrees entre et T e en prenant successivement les n vecteurs propres de R 
comme conditions initiales. 

Meme pour le cas simple du pendule parametrique (on parle d'equations de 
Mathieu), la determination analytique des exposants de Floquet fait l'objet 
de developpements mathematiques complexes (determinants infinis, ...). Dans 
l'espace de controle (A, —,h), on montre que l'equilibre se destabilise prin- 
cipalement pour to/uJ e ~ 1/2, sur un intervalle d'autant plus large que h est 
grand, a condition de depasser une valeur une valeur critique d'autant plus 
grande que la dissipation A est grande. En plus de la resonnance parametrique 
uj/uj e ~ 1/2, on observe les resonances uj/uj e ~ p avec p G IV*, mais sur des 
domaines de parametres d'autant plus restreints que p est grand. 
On peut effectuer un developpement asymptotique lorsque Pintensite h du 
forcage, le coefficient A de la dissipation et le detuning 5 = uj 2 /uj^ — 1/4 sont 
petits. On suppose que h = eF, A = eAet5 = Ae 

u(t) = y/e u (t,et) + e ui{t,et) + o(e 2 ) (4.46) 

comme precedemment. A l'ordre dominant, on obtient uo(t, T) = A(T)e tuJot + 
A* (T)e~ tuJot avec ujq = co e /2. La reponse rapide du pendule s'effectue done a 
la pulsation moitie de celle du forcage, e'est-a-dire avec une periode double. 
L'ordre suivant necessite d'imposer une relation de compatibilite qui conduit 
a l'equation d'amplitude 

A A 

— = (- A + iiy ) A + fi A-i[3\A\ 2 A (4.47) 
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avec v = 2A wq, u = loq F/A et j3 = luq/4. L'etude de cette equation, en posant 
par exemple A = Rexp(i6), permet de tracer le diagramme de stabilite dans 
le plan iy,F) pour differentes valeurs de A. L'instabilite est supercritique 
pour v > et sous-critique pour v < 0. 



F 






instable / 


sous- 




critique 


\ / 


stable 

i * 



v 
Figure 1: stabilite de I'equilibre du pendule parametrique 

4 Exploration numerique de systeme dynamiques 

L'objectif de ce paragraphe est de motiver la programmation de systemes 
dynamiques pour observer numeriquement la richesse et l'universalite des 
comportements identifies par Papproche theorique. Dans cet esprit, nous 
rappelons quelques methodes numeriques simples pour 1' approximation des 
equations d'evolution. Nous enumerons ensuite une liste de systemes dy- 
namiques qui offre un large champ d' experimentations numeriques. Cette 
liste est completee par les exercices et problemes presentes a la fin. 

4.1 Methodes numeriques simples 

On considere un systeme d'equations differentielles ordinaires couplees de la 
forme 

X = F{X) oil X{t) G M N . (4.48) 
En notations indicees ce systeme s'ecrit 

X i = F i (X 1 ,X 2 ,...,X N ) pourz = l,2,..,iV. (4.49) 

On suppose que F est differentiable. Le theoreme de Cauchy indique alors 
qu'une condition initiale X(0) donne naissance a une et une seule solution 

On presente ici quelques schemas d' approximation des solutions de ces systemes 
en discretisant le temps en pas de temps egaux a At. On considere alors une 
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succession d'etats X°, X 1 , X 2 , 2L n , ••• qui doivent approximer une solution 
X(t) aux temps t n = n At. 

Le schema d'Euler consiste a effectuer l'iteration 

vn+1 vn 

~ At ~ = F(2L n ) (4.50) 

a partir de la condition initiale X® = X(0). 
Le schema Euler arriere s'ecrit 

vn+1 vn 

~ At~ =^ n+1 ) (4-51) 

et necessite done de resoudre une equation implicite pour X a+l a chaque pas 
de temps. 

Le schema Leap-Frog (saute-mouton) est defini par Pequation explicite 

vn+1 _ vn—\ 

- - £(i "> (4 - 52) 

mais il fait intervenir les deux pas de temps precedant le nouveau pas de temps 
calcule. II doit done etre initialise avec un schema a un pas de temps (par 
exemple Euler). II est bon de savoir qu'il a tendance a decoupler les pas de 
temps pairs et impairs, ce qui se traduit par des oscillations numeriques que 
l'on peut eliminer en melangeant, a intervalles reguliers ou systematiquement, 
les pas de temps pairs ou impairs. 

Le schema Adams-Bashforth fait lui aussi intervenir deux pas de temps, 
sans presenter le phenomene d'oscillation numeriques du schema precedent. 
II s'ecrit 

yn+l v-n q i 

~ At ~ = 2 l_{X n ) - 2 F(2C^) . (4.53) 

Le schema "predicteur-correcteur" fait intervenir une etape intermediaire 
issue du schema d'Euler. II s'ecrit 

vl vn vn+1 vn 

^^=E(X n ) et ~ Af - =F(X I ). (4.54) 

Un schema bien connu pour sa precision est le schema de Runge-Kutta 
quatrieme ordre qui s'ecrit 

vl vn vll vn 

At - K — ' At/2 - y — ' 

vIII vn yIV vn 

At/2 ' At - v — ' 

X n+1 = X n + ^- [F{X!) + 2F{X n ) + 2F{X HI ) + F(X IV )\ (4.55) 
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Notons que Ton peut combiner plusieurs de ces schemas pour une meme 
equation. Si on decompose par exemple F(X) = A 2L + G(X) oil A est 
un operateur lineaire, le schema "Leap-Frog / Crank-Nicolson" s'ecrit 

yn+l vn—1 vn+1 , vn—\ 

Tous ces schemas peuvent se mettre sous la forme X n+1 = C(2L n , X n ~ l ; At), 
apres avoir resolu les equations implicites lorsque c'est le cas. On peut alors 
rappeler sommairement les notions de consistance, stabilite et convergence 
d'un schema. 

Un schema est consistant si 

Lim \\X(t + At)-C[X(t),X(t-At);At]\\ = 
At -> (4.57) 

pour une solution X[t). II est precis a l'ordre p si cette quantite tend vers 
zero comme (At) p . 

Un schema est stable si pour tout temps T il existe une borne superieure M 
finie telle que \\X n \\ < M pour tous les At et tous les n tels que nAt < T. 
Enfin, un schema est convergeant si pour toutes les conditions initiales 
l'approximation X™ converge vers la solution X(nAt) quand At tend vers 
zero et n vers l'infini avec t = nAt fixe. 

Le theoreme de Lax indique que pour un probleme bien pose et pour un 
schema consistant, la stabilite est equivalente a la convergence. 
On montre que les schemas d'Euler et Euler arriere sont precis au premier 
ordre, que les schemas de Crank-Nicolson, Leap-Frog et predicteur correcteur 
sont precis au second ordre, et que le schema de Runge-Kutta est precis au 
quatrieme ordre. 

L'analyse de la stabilite des schemas est du ressort de l'etude des systemes 
dynamiques discrets, encore appeles "mappings". 



4.2 Le zoo des systemes dynamiques continus 

On presente ici plusieurs specimens de systemes dynamiques c'est-a-dire de 
systemes d'equations differentielles ordinaires couplees. Ces exemples peu- 
vent servir de champ experimental pour des simulations numeriques ainsi que 
d'illustration des concepts qui seront developpes dans ce cours. 
Cage des formes normales 

On presente ici quatre systemes dynamiques a un degre de liberte, c'est-a- 
dire ne comprenant qu'une seule equation differentielle ordinaire d'ordre 1. 
Ces equations sont etudiees en detail dans le chapitre intitule "Bifurcations 
generiques de l'equilibre" . Ce sont des formes normales dans la mesure ou 
elles sont representatives de bifurcations de l'equilibre observees dans des 
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systemes dynamiques plus complexes pouvant avoir un grand nombre de 
degres de liberte. La theorie des formes normales montrent qu'en effectu- 
ant un developpement asymptotique au voisinage de la bifurcation on ramene 
la dynamique du systeme a celle de l'une ces formes normales. 

1. Bifurcation noeud-col : x(t) G M 

x = fi + a x 2 . (4.58) 

2. Bifurcation fourche : x(t) £ M 

x = fi x + a x 3 . (4.59) 

3. Bifurcation de Hopf : z(t) G W 

z = (fi + iuj)z + {a + i[5)\z\ 2 z . (4.60) 

4. Catastrophe en fronce : x(t) G M 

x = —q + px — x 3 . (4.61) 

Cage des vedettes du chaos 

On presente maintenant deux systemes dynamiques a trois degres de liberte 
souvent utilises pour illustrer le concept de chaos deterministe. En simu- 
lant numeriquement ces systemes pour certaines valeurs des parametres de 
controle, on montre que le comportement des solutions est erratique, sensible 
aux conditions initiales et que les trajectoires simulees decrivent un ensemble 
que l'on appelle souvent "attracteur etrange". L 'etude du premier exemple 
fera l'objet du chapitre "Modele de Lorenz" . 

1. Modele de Lorenz : X(t) = [x(t),y(t), z(t)\ G M 3 

x = —a x + a y avec a > 
y = —x z + r x — y avec r > 
z = x y — b z avec b > . (4.62) 

Le choix de Lorenz dans son article de 1963 est a = 10 et f3 = 8/3. Pour 
r suffisamment grand, ce systeme a un comportement chaotique. 

2. Modele de Rossler : X(t) = [x(t),y(t), z(t)} G M 3 

x = —y — x 
y = x + ay 

z= b + xz — cz (4.63) 

Un choix de parametres conduisant a un comportement chaotique est 
a = b = 0.2 et c = 5. 
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Cage des oscillateurs 

On presente un certain nombre de systeme dynamique exprimes sous la forme 
d'une equation differentielle d'ordre 2. En posant v = u, on se ramene a un 
systeme dynamique d'ordre 1 a deux degres de libertes. Ces systemes ont 
souvent un comportement oscillant, et decrivent done des oscillateurs. Leur 
etude fait l'objet du chapitre intitule "Oscillateurs" . 

1. Pendule simple : u(t) G M 

il + uj 2 s'mu = . (4.64) 

2. Oscillateur de Duffing : u(t) G M 

il + uJ 2 u = en 3 . (4.65) 

3. Pendule lineaire amorti : u(t) G M 

u + 2\u + lo 2 u = 0. (4.66) 

4. Pendule tournant : u(t) G M 

u + 2A u + {uj 2 n — uj 2 cos u) sin u = . (4.67) 

5. Pendule resonant : u(t) G M 

u + 2A u + uj 2 sin u = f cos(u; e t) . (4.68) 

6. Oscillateur de Van der Pol : u(t) G M 

u + (ku 2 -2/j)u + uj 2 1 u = . (4.69) 

7. Oscillateur force : u(t) G M 

u + 2\ii + uj 2 n s\nu = f cos(uj e t) . (4-70) 

8. Oscillateur parametrique : u(t) G M 

u + 2Xu + ujI [1 + h cos(uj e t)} sin u = . (4.71) 
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Conclusion 

Plusieurs elements de base dans l'etude des systemes dynamiques ont ete 
aborde a travers l'examen du modele de Lorenz, du pendule simple et autres 
oscillateurs. La linearisation de ces systemes autour de leurs equilibre a ete ef- 
fectuee systematiquement. L'etude du systeme lineaire ainsi obtenu s'effectue 
a l'aide des vecteurs propres et valeurs propres de la matrice jacobienne. On 
obtient ainsi la stabilite de P equilibre etudie ainsi que la topologie des trajec- 
toires autour de cet equilibre. 

Les changements de stabilite et de topologie des trajectoires obtenus en faisant 
varier un parametre de controle sont appeles des bifurcations. Nous avons 
ainsi observe une "bifurcation fourche" lors de la destabilisation de Pequilibre 
trivial du modele de Lorenz et une "bifurcation de Hopf" pour la destabilisation 
du repos dans l'oscillateur de Van der Pol. Ces bifurcations sont propres aux 
systemes dissipatifs qui modelisent la competition entre un forgage exterieur 
et une dissipation interne. 

Les systemes conservatifs modelisent des systemes physiques isoles de tout 
forgage et pour lesquels on neglige la dissipation d'energie. Les systemes 
hamiltoniens forme une large classe de modeles conservatifs qui incluent les 
sytemes dependant d'un potentiel et le pendule simple. 

A travers les quelques exemples presentes, plusieurs concepts importants ont 
ete rapi dement abordes : sensibilite aux conditions initiales, phenomene de 
resonance, instabilite parametrique, etc. Le calcul d'equations d'amplitude 
par la methode des echelles a ete rapi dement esquisse sur deux exemples. 
Enfin, une initiation a l'exploration numerique des systemes dynamiques a 
ete effectuee. 

Cette introduction est une invitation a approfondir l'etude des systemes dy- 
namique dont la comprehension est essentielles pour comprendre la dynamique 
de nombreux systemes physiques, en particulier dans le cas des la mecanique 
des fluides. 



FORMULAIRE 

EXERCICES ET PROBLEMES 



EXERCICE 4.1 Systemes conservatifs 



Tracer les portraits de phase des systemes dynamiques suivants 

1) Systeme 1 : u = —u(u 2 — 1). 

2) Systeme 2 : u = -(u - l)(u - 2)(u - 3) 

3) Systeme 3 : x = —x + 2y et y = —x + y. 
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4) Systeme 4 : q = 2pq/r — p s /r 2 et p = Ap — p 2 /r 



EXERCICE 4.2 Modele de Rossler 



On considere le systeme dynamique suivant : 



± = -y-z (4.72) 
V = x-Zy (4.73) 
z = -fi + xz . (4.74) 

Ce systeme est un cas particulier du modele de Rossler (a = —3, b = —fi et 
c = 0) qui decrit la cinematique d'une reaction chimique entre des constituants 
dont les concentrations sont les quantites x, y et z. 

1) Indiquer les ensembles que Ton appelle espace des phases et espace de 
controle pour ce systeme dynamique. 

2) Calculer les points d'equilibre en fonction des valeurs des parametres de 
controle. 

3) Effectuer le calcul de stabilite pour le point d'equilibre correspondant au 
cas fi = 0. 

4) Tracer le portrait de phase (trajectoires) restreint au plan Oxy au voisi- 
nage de ce point d'equilibre pour fi = 0. 

5) Dans le cas general /i / 0, montrer que la stabilite des points fixes est 
determinee par les racines d'un polynome de degre trois que l'on calculera. 

6) Calculer les racines de ce polynome lorsque l'une d'entre elles est nulle. 
Relier ce resultat avec l'une des questions precedentes. 

7) Calculer la stabilite des equilibres dans un voisinage dans le voisinage de 
la valeur critique fi = 0. 

8) Indiquer la nature de la bifurcation observee et tracer son diagramme. 

9) Le systeme dynamique est-il invariant par la symetrie point par rapport 
a O ? Relier cette reponse a la nature de la bifurcation. 



PROBLEME 4.3 



Exemples de bifurcations 



Bifurcations des systemes dynamiques conservatifs 

1) Discuter, en fonction du parametre de controle fj,, le comportement du 
systeme dynamique reel 

x = fix + y (4-75) 
y = -x - fiy . (4.76) 
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Stabilite du point fixe (solution stationnaire), portraits de phases (traces 
des trajectoires dans l'espace des phases), solutions analytiques, dia- 
gramme de bifurcation (illustrant les changements de dynamique en fonc- 
tion du parametre de controle), exposants critiques, etc. ... 
2) Discuter le trace de la famille de courbes definie par l'equation 

\(p 2 + q 2 )+m = E (4.77) 
dans le plan reel (p, q) (relier cette question a la precedente) . 



Bifurcations des systemes dynamiques dissipatifs 

Etudier les systemes dynamiques suivants (stabilite des points fixes, portraits 
de phases, solutions analytiques, diagrammes de bifurcations, exposants cri- 
tiques, brisure de symetrie, etc. ...) pour les valeurs du parametre de controle 
/U indiquees : 

3) Pour (x, y) G M 2 et ^ G M voisin de \i c = 1 : 

x= -2x + y (4.78) 
y = -2-xy + 2fxy. (4.79) 

4) Pour Ae Met fi£ M voisin de fi c = : 

A + 3A + 2A-2fiA + A 3 = . (4.80) 



5) Pour z e (F et n £ M : 

z = (n + i) z + (2 + i)\z\ 2 z - (1 + 2i)\z\ A z . (4.81) 



PROBLEME 4.4 



Transitions de phases : une analogie 



On considere un milieu continu dont les variables d'etat (pression, tempera- 
ture, etc.) sont proches d'une transition de phase liquide-vapeur. Le but de ce 
probleme est d'etudier un modele unidimensionnel simple dont le comporte- 
ment permet d'effectuer une analogie avec une telle transition de phase. 



O 



liquids 



\ gaz 



Figure 4.12: Transition de phase dans un milieu continu unidimensionnel. 

Dans cette analogie, le milieu est contenu dans un tube fin d'axe Ox et le 
parametre d'ordre reel A(x, t) £ 1R de la transition obeit a l'equation aux 
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derivees partielles suivante : 



8A 

~dt 



I^A-a A 3 - A 5 + 



d 2 A 



(4.82) 



dx 2 



oil a £ M depend des variables d'etat ainsi que ji G M qui est est proportionnel 
a l'ecart a la temperature de liquefaction. Dans ce modele simple, on suppose 
que A = represente la phase liquide et que les solutions A / sont associees 
a la phase gazeuse. 

Transition du second ordre 

On suppose ici que a = 1 et on s'interesse uniquement aux solutions ho- 
mogenes en espace de la forme A(x, t) = R(t). 

1) Montrer que R(t) est regi par un systeme dynamique 



et donner l'expression de la fonction F + qui depend du parametres \x. 

2) Tracer le diagramme de bifurcation de ce systeme dynamique en fonction 
du parametre de controle /x. 

3) Commenter l'analogie avec une transition de phase du second ordre liquide- 
vapeur. 

Transition du premier ordre 

On suppose maintenant que a = — 1 et on note 



le systeme dynamique regissant revolution des solutions homogenes d'espace 
A(x,t) = R(t). 

4) Tracer le diagramme de bifurcation de ce systeme dynamique en fonction 
du parametre de controle /x. 

5) On reserve desormais le qualitatif de phase gazeuse aux valeurs de A 
proches d'une solution stationnaire R(t) = R g ^ stable pour le systeme 
dynamique (4.84). Calculer la valeur /x* du parametre [x delimitant 
Pintervalle [//*, 0] pour lequel la phase liquide R\ = et une phase gazeuse 
R g / coexistent. 

6) On s'interesse maintenant aux solutions stationnaires A(x, t) = Q(x). 
Montrer que le profil Q(x) est solution d'un systeme dynamique de la 
forme 



R = F+(R) 



(4.83) 



R = F-(R) 



(4.84) 




(4.85) 



et donner l'expression du potentiel V en fonction de fi. 
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7) Tracer les allures de la fonction V(Q) lorsque fi varie, en particulier dans 
l'intervalle [/x*,0]. On s'interessera surtout a la comparaison des valeurs 
des extrema du potentiel V. 

8) En deduire qu'il existe une valeur fi c £ [/i*, 0] du parametre fj, pour laque- 
lle la phase liquide et la phase gazeuse ont le meme potentiel, et calculer 

fJ-c- 

9) Tracer le portrait de phase du systeme dynamique (4.85) pour fi = fi c . 
Commenter l'existence d'une trajectoire heterocrine Qh{x) (i.e. reliant 
deux points fixes) entre la phase liquide (Q = Ri = 0) et la phase gazeuse 
(Q = R g / 0) dans le cadre de l'analogie avec une transition de phase. 

10) Dans le cadre de cette analogie on appelle "densite de potentiel thermo- 
dynamique" associee a A(x,t), la fonction W^[A(x,t)] = — y[yl(x,t)] + 
Kfx) 2- ^ e P ms ' on dit qu'une phase homogene stationnaire, c'est-a-dire 
A(x,t) = Ri ou A(x,t) = R g , est "metastable" si elle est stable au regard 
des perturbations homogenes mais si son "potentiel thermodynamique" 
W(R) = —V(R) est plus eleve que le "potentiel thermodynamique" de 
l'autre phase qui est alors qualifiee de "phase stable" . Reformuler le di- 
agramme de bifurcation du systeme dynamique (4.84) en precisant cette 
fois la "stabilite" (trait plein), "metastabilite" (trait tirete) ou "insta- 
bilite" (trait pointille) des differentes phases. 



Construction de Maxwell 

On suppose toujours que a = —1. On s'interesse maintenant aux solutions 
de type front s'exprimant sous la forme A(x, t) = B(x — ct) ou la constante c 
est la vitesse de propagation du front et ou B(— oo) = Ri = 0. 

11) Montrer que B(y) est solution du systeme dynamique 

B dB , , , 
+ c— + V (B) = (4.86) 



dy dy 

et interpreter le parametre c de ce systeme en terme de coefficient de 
frottement. 

12) On suppose que /x G [//*, fj, c ] est fixe et Ton fait varier c a partir de la valeur 
c = 0. Montrer qualitativement qu'il existe une valeur c = C(/i) pour 
laquelle le systeme dynamique (4.86) admet une trajectoires heteroclines 
Bh(y) reliant la phase liquide Bh^—oo) = Ri = a la phase gazeuse en 
l?h(+oo) = R g . Indiquer le signe de C(fj.) pour /i G [^*,/i c ]. Merries 
questions pour /j, £ [/U c ,0]. 

13) En multipliant l'equation (4.86) par -^-B^y) et en integrant de — oo a 
+oo, calculer C(/u) en fonction de la difference de potentiel thermody- 
namique entre les phases liquide et gazeuse et du terme J<^°°[^{y)] 2 dy. 
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14) Deduire de l'une des deux questions precedentes (ou des deux) qu'un front 
se propage toujours de la phase metastable vers la phase stable. 

15) Commenter l'analogie entre la transition stable-metastable a p = p c et 
une transition de phase liquide-vapeur du premier ordre. 



Point critique 

On considere maintenant le cas general a GR. 

16) Montrer que suivant le signe de a on peut se ramener, par un changement 
de variable simple, au cas a = 1 (transition du second ordre) ou a = — 1 
(transition du premier ordre). 

17) Montrer que le point (p, a) = dans l'espace de controle de l'equation 
(4.82) peut etre qualifie de point critique pour l'analogie avec la transition 
de phase liquide vapeur. 

18) Completer eventuellement ce probleme : exposants critiques des transi- 
tions, brisure de symetrie, formulation variationnelle de (4.82), fonction- 
nelle de Liapounov, ... 



PROBLEME 4.5 



Experimentations numeriques 



Les systemes dynamiques suivants complement la liste des exemples men- 
tionnes mentionnes dans le corps du texte. On pourra en faire l'etude de 
maniere analytique (equilibres, stabilite, developpements asymptotiques, ...) 
ainsi que numerique. 

1) Oscillateurs de Lienard : u(t) € M 

u + (-e + r u 2 ) u + (cjq u + p u 3 ) = (4.87) 

2) Oscillateur force : u(t) G M 

u + -fii + Lolu = f + ku 2 (4.88) 

3) La boussole : u(t) G M 

u + aii + Msinu- Psm(u-t) = . (4.89) 

4) Equations de Volterra : X(t) = [x(t),y(t)\ G R 2 

x = A x — a x y 

y = -py + a x y . (4.90) 

5) Equations de Volterra complexes : Z(t) = [A(t),B(t)\ £ d 2 et (p, a, (3) £ 

A= p A + (a \A\ 2 + \B\ 2 ) A 

B= p B + (a \B\ 2 + j3 \B\ 2 ) B . (4.91) 
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6) Calcul de forme normale : X(t) = [A(t),B(t)} G M 2 et (n,\,a,/3) G M 4 

A= n A + a A B 

B= -XB + I3A 2 . (4.92) 

7) Systeme hamiltonien "Todda Lattice" : (q,p) G JR 2n 

, n n 

H (g,p) = ~T,p 2 j + E ex p(^-i - n) ( 4 - 93 ) 

3=1 3=1 

8) Systeme hamiltonien "Henon-Heiles" : (q, p) G JR 4 

^(£,P) = ^(Pi +Pl) + \ [fi + 92 + 2 gfe - \ql) (4.94) 

9) Etudier l'equation d'amplitude 

dA _ 

= (-A + iz/)A + ^4-i/3|A| 2 .4 (4.95) 



et comparer avec l'equation du pendule parametrique 

u + Att + LjQ [1 + h cos(u; e t)] u = . (4.96) 

dont elle est issue avec f = 2Aujo, M = ujqH/4, (3 = wo/4, wo = w e /2, 
5 = lo 2 Juj 2 - 1/4 = eA, A = eA, h = eH et 

u(t) = yfe [A(et)e iu}ot + Z(et) e -^°* + O(e)] . (4.97) 

CORRIGES DES EXERCICES ET PROBLEMES 
QUESTIONNAIRES A CHOIX MULTIPLES 

En chantier 

CORRIGES DES QCM 

En chantier 
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